Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,159	0	G	A	42.1%	89.3 / 46.0	76	intergenic (+139/‑75)	thrC/yaaX	L‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (38/6);  new base A (1/31);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 3.95e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5013‑5308                                                                                                                                                       |                                                                                                                                                      atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑‑AAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATaa                                                                                                                                                  <  2:140266/151‑1 (MQ=255)     gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                               <  2:556915/151‑1 (MQ=255)         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                           >  2:1091084/1‑151 (MQ=255)         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                           <  2:398540/151‑1 (MQ=255)                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑AGAAAAGG‑GAAAAGGAGAAATTTTTAATTAATGCGGGAA‑CTT‑TGGGATTTAGGTTGCGGGAGATAACCAacgccggt                                                                                                                                       >  2:1040974/1‑145 (MQ=255)                    tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                                >  1:225031/1‑151 (MQ=255)                          aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                          >  2:428248/1‑151 (MQ=255)                           aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                         >  2:517737/1‑151 (MQ=255)                           aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                         >  1:97678/1‑151 (MQ=255)                            tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                        >  2:906386/1‑151 (MQ=255)                               ggcggttttttttttaATGCAGCCCGGTTTTTTTAAGAAAAAATTAGGAAAAAAAA‑TAACAGGA‑AAAAAGGAAAAATTCAAAAAAAAGGGGTTAT‑CTG‑GGAGATTGGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTaa                                                                                                                           <  2:519519/133‑1 (MQ=255)                                gCCGGTTTTTTTTTTAGGCAGCCCGGTTTTTTTAAGAAAAAATTAGGAAAAAAAA‑TGACAGGA‑AAAAAGGAAAAATTCAAAAAAAAGGCGTTAT‑CTT‑AGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                    <  1:1007974/151‑1 (MQ=255)                                      tttGTTTTTAGGCACCCCGTTTTTTTTAAAAAAAAATTAGGGAAAAAAA‑TGAAAGGA‑AAAAAGGAAAAATTCTAAATAAATGCGTTAA‑CTT‑AGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                               <  2:64739/150‑1 (MQ=255)                                      tttGTTTTTAGCCAGCCCGTTTTTTTTAGGAAAAATTTAGGAAAAAAAA‑GAACAGGA‑AAAAAGGAAAAATTCTAAATAATTGGGTTAA‑TTT‑AGAGATTAGTATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                  <  2:44161/147‑1 (MQ=255)                                       ttGTTTTTAGGCGGCCCGGTTTTTTTAAAAAAAAATTAGGAAAAAAAA‑TGACAGGA‑AAAAAGAAAAAATTCTAAAAAAAGGCGGAAA‑CTA‑AGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                              <  2:63334/150‑1 (MQ=255)                                          tttttAGGCAGCCCGGTTTTTTAATAAAAAAATTAGGAAAAAAAA‑GAAAAGGA‑AAAAAGAAAAAATTCAAAAAAAAGGCGTTAT‑CTT‑AGAGATTAGGATGGCGAAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                           <  1:466326/149‑1 (MQ=255)                                            tttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        >  2:744648/1‑151 (MQ=255)                                             ttATGCAGCCCGTTTTTTTTAGAAAAAAATAAGGGAAAAAAA‑AGAAAGGA‑AAAAAGAAAAAATTCTAAAAAAAGGCGTTAA‑CTT‑AAAAATGGGGATTGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                         <  1:791585/149‑1 (MQ=255)                                              tATGCAGCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAA‑TGAGAGGA‑AAAAAGAAAAAATTCTAAAAAAATGCGGTAA‑TTT‑AGAGATGGGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                       <  1:1178234/150‑1 (MQ=255)                                               acgCACCCCGTATTTTTAAAAAAAAAATTAGAGAAAAAAA‑TAAGAGGA‑AAAAAGAAAAAATTCTCAAAACATGCGATAC‑TTT‑AAAAATTGGGATTGGGAAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        <  1:654405/146‑1 (MQ=255)                                               aTGCAGCCCGGTTTTTTAAAGAAAAAATTAGGAAAAAAAA‑GGACAGGA‑AAAAAGGAAAAATTCAAAAAAAAGGGGTAAA‑CTA‑AGAGATAAGTATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                      <  1:1101121/150‑1 (MQ=255)                                                tGCACCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAA‑TGACAGGA‑AAAAAGGAAAAATTCTAAAAAAATG‑TGTAA‑CTTAAGAGATTAGTATGGCGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATaa                                                                                                      <  2:216331/149‑1 (MQ=255)                                                  cAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                  >  2:382857/1‑151 (MQ=255)                                                   aGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                 >  2:554341/1‑151 (MQ=255)                                                     ccggttttttttATAAAAAAATAAGGGAAAAAAA‑GAAGAGGA‑AAAAAGGAAAATTTCTAAAAAAAGGGGATAA‑TTC‑AGAGATTAGGGTGGGGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                               <  2:507795/145‑1 (MQ=255)                                                     ccggtttttttaATGAAGAAATT‑TGGAAAAAAAGTGAAAGGA‑AAAATGGAAAAATTTTAAAAAAAGGCGGTAT‑TTT‑AAAGATAGGGTTGGCGGAGAAAAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                               <  2:137815/145‑1 (MQ=255)                                                        gttttttttATAAAAAAATTAGGAAAAAAAA‑TGACAGGA‑AAAAAGGAAAAATTCTAAAAAAAGGCGGTAA‑TTT‑AGATATTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                             <  1:1010278/147‑1 (MQ=255)                                                        ggtttttttAGAAAAAAATTAGGAAAAAAAA‑GGAAGGGA‑AAAAAGGAAAAATTCTAAAAAAAGGCGGTAA‑TTA‑AGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                             <  2:716880/147‑1 (MQ=255)                                                          ttttttaATAAAAAAATTAGGAAAAAAAA‑TGACAGGA‑AAAAAGGAAAAATTCTAAATAAAGGGGTTAA‑TTT‑AGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                             <  2:369226/147‑1 (MQ=255)                                                             ttttATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        >  2:931457/1‑150 (MQ=255)                                                               taaaaaaaaaaTTAGGAAAAAAAA‑AGAAAGGA‑AAAAAGAAAAATTTCAAAAAAATGGGGGTAA‑TTA‑AGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                      <  1:807580/142‑1 (MQ=255)                                                                  tGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                  >  1:666371/1‑151 (MQ=255)                                                                   aaaaaaaTTAGGGAAAAAAA‑TGAAAGAA‑AAAAAGGAAAAATTCTAAAAAATTGCGTTAA‑TTT‑AGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                 <  1:424091/147‑1 (MQ=255)                                                                    aaaaaaTTAGGAAAAAAAA‑GGACAGGA‑AAAAAGGAAAATTTCTAAATAAATGCGTTAT‑CTT‑AGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                   <  2:1195904/145‑1 (MQ=255)                                                                        aataagggAAAAAAA‑GGACGGGA‑AAAAAGAAAAAATTCAAAATAAAGGGGTTAA‑CTT‑AGAATTAAGGATGGGGAAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            <  1:479367/145‑1 (MQ=255)                                                                        aaTTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            >  2:944874/1‑151 (MQ=255)                                                                        aaTTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            >  2:881382/1‑151 (MQ=255)                                                                            agggAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                         >  1:500653/3‑150 (MQ=255)                                                                             tGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                        >  1:643900/1‑150 (MQ=255)                                                                               ggaaaaaa‑aGAAAGGGAAAAAAGGAAAAATTCTAAAAAAAGGGGTAAA‑ATT‑AGAGATTAGGTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                          <  2:456053/144‑1 (MQ=255)                                                                                agaAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                      >  2:668438/1‑149 (MQ=255)                                                                                aaaaaaa‑TGAAAGGA‑AAAAAGGAAAAATTCTAAATAAAGGCGATAA‑TTA‑AAAAATTAGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                    <  1:950849/149‑1 (MQ=255)                                                                                aaaaaaa‑GGACAGGA‑AAAAAGGAAAAATTCAAAATAAATGCGGTAA‑TTA‑AAAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa                                                                     <  2:500534/148‑1 (MQ=255)                                                                                    aaa‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                                >  1:587449/1‑151 (MQ=255)                                                                                     aa‑TGAAGGGA‑AAAAAGGAGAAATTCAAAATAAATGCGAAAT‑CTA‑AGAGTTTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                               <  2:1105438/151‑1 (MQ=255)                                                                                        cgccAGGG‑AAAAAGGAGCAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                               >  1:290177/4‑149 (MQ=255)                                                                                          atagga‑aaaaaggaaaaaTCATAAATAATTG‑GGTAATGTT‑AAATATTAGATGTGGGAAAAATAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                           <  2:1025536/149‑1 (MQ=255)                                                                                           cAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                           >  2:735469/1‑150 (MQ=255)                                                                                                aaaaaaGGAAAAATCCTAAAAAAATGCGTTAA‑CTA‑AGAGATGAGAGTGGGTGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                             <  1:234537/143‑1 (MQ=255)                                                                                                 aaaaaGGAAAAATTCTAAAAAAGGGGGTTAA‑TTT‑AGAGTTAAGTATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                        <  1:891778/148‑1 (MQ=255)                                                                                                 aaaaaGGAAAAATTCTAAAAAAATGCGTTAA‑CTA‑AAAAATTAGATTTGAGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                           <  1:185763/145‑1 (MQ=255)                                                                                                 aaaaaGAAAAAATTCAAAAAAAAGGCGTTAA‑TTT‑AGAGATTAGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                      <  1:642636/150‑1 (MQ=255)                                                                                                  aaaaGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                    >  1:14442/1‑151 (MQ=255)                                                                                                  aaaaGGAAAAATTCAAAATAAATGCTAAAA‑TTA‑AGAATTTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                    <  2:997077/151‑1 (MQ=255)                                                                                                   aaaGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                   >  1:1184524/1‑151 (MQ=255)                                                                                                      ggAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                                >  1:182881/1‑151 (MQ=255)                                                                                                          aaaTTCTAAAAAAAGGCGGTAA‑CTT‑AGAGATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             <  1:193463/150‑1 (MQ=255)                                                                                                          aaaTTCTAAAAAAAGGCGGTAA‑CTA‑AGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             <  2:494524/150‑1 (MQ=255)                                                                                                           aaTTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                           >  1:501379/1‑151 (MQ=255)                                                                                                           aaTTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                           >  1:560230/1‑151 (MQ=255)                                                                                                             ttCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                         >  2:136731/1‑151 (MQ=255)                                                                                                                 cAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                     >  2:1092572/1‑151 (MQ=255)                                                                                                                  aaTAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                    >  1:82271/1‑151 (MQ=255)                                                                                                                  aaTAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                    >  2:1187227/1‑151 (MQ=255)                                                                                                                    tAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                  >  2:486548/1‑151 (MQ=255)                                                                                                                    tAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                  >  2:241374/1‑151 (MQ=255)                                                                                                                     aaaTGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                 >  1:938791/1‑151 (MQ=255)                                                                                                                        tGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                              >  2:286232/1‑151 (MQ=255)                                                                                                                              aa‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                        >  1:5411/1‑151 (MQ=255)                                                                                                                                     agagATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                   >  1:13479/1‑151 (MQ=255)                                                                                                                                       agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAACGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                  >  2:976693/1‑150 (MQ=255)                                                                                                                                         aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               >  2:894709/1‑151 (MQ=255)                                                                                                                                          ttgggttgggggAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               <  2:817238/140‑1 (MQ=255)                                                                                                                                          ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               >  1:137757/1‑150 (MQ=255)                                                                                                                                          tagggATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                     <  2:107060/141‑1 (MQ=255)                                                                                                                                                   cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg     >  1:845406/1‑151 (MQ=255)                                                                                                                                                       gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaa  >  1:5686/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAA‑TGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAA‑CTT‑AGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5013‑5308

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑