Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,997,022	G→C	8.2%	noncoding (52/93 nt)	USA300HOU_t0037 ←	tRNA‑Ser

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,997,022	0	G	C	8.2%	246.1 / 11.0	82	noncoding (52/93 nt)	USA300HOU_t0037	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base G (25/50);  new base C (4/3);  total (29/53)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CTATAAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTCCTGCGTGACAGG  >  CP000730/1996883‑1997157                                                                                                                                            |                                                                                                                                        cTATAAAATATTGCGGGAGGCGGATTTGAACCCCCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGc                                                                                                                                        <  2:788675/141‑1 (MQ=255) cTATAAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGc                                                                                                                                        <  2:607808/141‑1 (MQ=255) cTATAAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGc                                                                                                                                        <  1:247355/141‑1 (MQ=255)  tataAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGcc                                                                                                                                       >  1:743008/1‑141 (MQ=255)  tataAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGcc                                                                                                                                       >  2:836417/1‑141 (MQ=255)    taAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTCAAGGCCCt                                                                                                                                     >  2:1032454/1‑141 (MQ=255)     aaaaTATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTg                                                                                                                                    >  1:611747/1‑141 (MQ=255)     aaaaTATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTg                                                                                                                                    <  1:326345/141‑1 (MQ=255)           ttGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGtt                                                                                                                              <  1:72584/141‑1 (MQ=255)            tGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGttt                                                                                                                             >  1:60018/1‑141 (MQ=255)            tGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGACGGctt                                                                                                                             >  1:80145/1‑141 (MQ=255)             gCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGtttt                                                                                                                            >  1:103548/1‑141 (MQ=255)             gCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGtttt                                                                                                                            >  2:435653/1‑141 (MQ=255)               gggAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCa                                                                                                                          <  2:931510/141‑1 (MQ=255)                 gAGGCGGAGTTGAACCACCGGCCTTCGGGTTATGAGCCCGACGAGCCACCGGGCTGCTCCATCCCGCTATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTGAGGCCCTGTCGGTTGTCAAg                                                                                                                        <  2:752298/141‑1 (MQ=37)                   ggCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGAc                                                                                                                      >  2:447394/1‑141 (MQ=255)                   ggCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCAGCGATACTCCAAAATAATGGCGGAGGCAGAGGGATTCGAACCCCAGCGGCAAGTTAAGGCCCTGTCGGTTTTCAAGAc                                                                                                                      <  1:432395/141‑1 (MQ=25)                         tttGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATcc                                                                                                                <  1:810456/141‑1 (MQ=255)                         tttGAACCACCGACCTTCATCGGGTGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATcc                                                                                                                <  1:826289/141‑1 (MQ=37)                         ttggAACCACCCACCTTCGGGTTAGGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATcc                                                                                                                <  2:650327/138‑1 (MQ=255)                           tGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCt                                                                                                              <  1:654525/141‑1 (MQ=255)                            gAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGCCCATGTCGGTTTTCAAGACCGATCCCtt                                                                                                             >  2:1018390/1‑141 (MQ=255)                                caccGCCCTTCGGTTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCTCGATAATAAAAAATAATGGCGGAGGAAGTGGTGTTCGAACCCCCGCTGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGc                                                                                                         <  2:38021/141‑1 (MQ=37)                                caccGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATcc                                                                                                                >  2:116370/1‑134 (MQ=255)                                caccGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATcc                                                                                                                <  1:116370/134‑1 (MQ=255)                                caccGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGc                                                                                                         <  2:878613/141‑1 (MQ=255)                                caccGACCCCCGGGTTACGAGCCCGACGGGCTACCGAACTGCCCCATCCCGCGATAGTACAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGc                                                                                                         <  2:60018/141‑1 (MQ=18)                                   cGACCGTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTACAAACAGCAATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCgg                                                                                                      <  2:493725/141‑1 (MQ=25)                                     aCCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGAc                                                                                                    <  1:250867/141‑1 (MQ=255)                                     aCCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGAc                                                                                                    >  2:558100/1‑141 (MQ=255)                                      ccTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACt                                                                                                   <  1:327367/141‑1 (MQ=255)                                        ttCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCTATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTg                                                                                                 <  1:556841/141‑1 (MQ=255)                                        ttCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATACAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTg                                                                                                 >  2:826282/1‑141 (MQ=255)                                        ttCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTg                                                                                                 <  2:288749/141‑1 (MQ=255)                                          cGGGGTATGAGCCCGACGAGCTACCGAACTGCGCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCACGACCGATCCCTTCAGCCGGACTTggg                                                                                               <  2:159204/141‑1 (MQ=255)                                            ggTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTa                                                                                             <  2:178279/141‑1 (MQ=255)                                               tATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtc                                                                                          <  2:819647/141‑1 (MQ=255)                                                aTGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCCGCCGGACTTGGGTATtcc                                                                                         >  1:1045828/1‑141 (MQ=255)                                                aTGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtcc                                                                                         <  2:1030110/141‑1 (MQ=255)                                                 tGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATtcct                                                                                        >  2:933038/1‑141 (MQ=255)                                                        gacgaGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCattat                                                                                 <  1:755577/141‑1 (MQ=255)                                                           gaGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGAAAtccgacgagactcg                                                                              >  2:741056/1‑130 (MQ=21)                                                                 ccGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAAAATtatatggac                                                                        <  1:558100/141‑6 (MQ=14)                                                                   gAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAAt                                                                      <  1:829734/141‑1 (MQ=255)                                                                        gCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTc                                                                                                                                  <  2:605526/76‑1 (MQ=17)                                                                        gCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGAGCCGTTAAGCCCCTGTc                                                                                                                                  >  1:605526/1‑76 (MQ=17)                                                                         cTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATATGTACCTCGcagg                                                                >  2:647338/1‑138 (MQ=255)                                                                         cTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTat                                                                <  1:737741/141‑1 (MQ=255)                                                                           ccttCCCGCGATAAGAAAAAGTGTTGGCGGAGCGTGAGGGGTGGGAGGCCGTGGTGCCGGTGAATGGCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTatta                                                              <  1:99129/138‑1 (MQ=12)                                                                           ccATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTatca                                                              <  2:763342/141‑3 (MQ=255)                                                                             aTCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTattaat                                                            <  1:771383/141‑1 (MQ=255)                                                                             aTCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTattaat                                                            >  1:152218/1‑141 (MQ=255)                                                                                  gcgATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATaa                                                       >  1:568700/1‑141 (MQ=255)                                                                                  gcgATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATaa                                                       >  1:568578/1‑141 (MQ=255)                                                                                  gcgATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATAATTATATGTAACTCGCTattcactaacc                                                       >  2:488436/1‑133 (MQ=18)                                                                                   cgATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATaaa                                                      <  2:72924/141‑1 (MQ=255)                                                                                   cgATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATaaa                                                      >  1:258015/1‑141 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  2:960054/141‑1 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  2:139947/141‑1 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  2:638053/141‑1 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  1:720807/141‑1 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  2:743008/141‑1 (MQ=255)                                                                                     ataataAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAAt                                                    <  2:856823/141‑1 (MQ=255)                                                                                         taaAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTaaa                                                <  1:622072/141‑1 (MQ=255)                                                                                              aactaaggCGGTGGAAGAGCGATTCGAGCCACGGCGGCCCGGTAAGGCCCTGTCGGGTTGCAAGACCGATCCCATCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAGATTAAATGGCg                                           <  2:40795/135‑1 (MQ=21)                                                                                               ataatGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCgg                                          <  1:892095/141‑1 (MQ=255)                                                                                               ataatGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCgg                                          <  1:547303/141‑1 (MQ=255)                                                                                               ataatGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGCTTTTTTTTACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTGTAGGTAAATCGCTATTAATTATAAAATTAAATGGCgg                                          <  1:782614/141‑1 (MQ=25)                                                                                                   tGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGtctc                                      <  1:230136/141‑1 (MQ=255)                                                                                                   ggtcGGAGGAAGAGGGATTCGACCGCCCGCGGCCCGTTAAGGCCCTGGCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGtctc                                      <  1:274238/138‑1 (MQ=255)                                                                                                    ggCGGAGGAAGAGGGATTCGAGCCCCCGCGGCCCGGTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCg                                     <  1:524359/141‑1 (MQ=255)                                                                                                        gaggaAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACgg                                 <  1:382673/141‑1 (MQ=255)                                                                                                                   aTTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTAATATAGGTAAATAGACCCTAAGAATAAAATTAAATGGCGGTCTCGACGGGAATCGAAccc                      <  1:520611/141‑1 (MQ=21)                                                                                                                    ttCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGTCGGGAATCGAACCcg                     >  1:899916/1‑141 (MQ=255)                                                                                                                        aaCCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGAt                 >  1:622405/1‑141 (MQ=255)                                                                                                                          cccccGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCAATATTATATGTAAATCGCTATTCTCGATCCAGCTACATGACGGTCTCCACGGGAATCGCACCCACGAtca               >  2:962950/1‑140 (MQ=17)                                                                                                                            cccGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTcc             <  2:611747/141‑1 (MQ=255)                                                                                                                            cccGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTcc             >  1:425512/1‑141 (MQ=255)                                                                                                                               gcgGCCCGTTAAGGCCCTATCGTGTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTCCTGc          <  1:826282/141‑1 (MQ=255)                                                                                                                                cgGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTCCTGCg         >  1:590995/1‑141 (MQ=255)                                                                                                                                       ttcaGGCCCTGTCGCTTTTCCAGACCGGCCCCTTCAGCCGCACCAAGGTATTCCTCCACTATTATAGGTAAAGCGCTATTAATTATAAGATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTCCTGCGTGAcagg  <  1:1082373/138‑1 (MQ=25)                                                                                                                                       ttAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATAATTATAGGTAAATCGCTATTATTTATACAATTAAATGGCGGTCTCGACGGGAATCGAACCCCCGATCTCCTGCGTGAccgg  >  2:79405/1‑141 (MQ=37)                                                                                                                                            |                                                                                                                                        CTATAAAATATTGCGGGAGGCGGATTTGAACCACCGACCTTCGGGTTATGAGCCCGACGAGCTACCGAACTGCTCCATCCCGCGATAATAAAAAATAATGGCGGAGGAAGAGGGATTCGAACCCCCGCGGCCCGTTAAGGCCCTGTCGGTTTTCAAGACCGATCCCTTCAGCCGGACTTGGGTATTCCTCCATTATTATAGGTAAATCGCTATTAATTATAAAATTAAATGGCGGTCTCGACGGGAATCGAACCCGCGATCTCCTGCGTGACAGG  >  CP000730/1996883‑1997157

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑