Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,486	0	C	A	45.7%	103.5 / 37.9	70	intergenic (+291/‑165)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (31/7);  new base A (0/32);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 2.26e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  NC_000913/29341‑29636                                                                                                                                                   |                                                                                                                                                       gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAaat                                                                                                                                                   <  2:323313/151‑1 (MQ=255)  ttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAaat                                                                                                                                                   >  2:581604/1‑150 (MQ=255)        tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                             >  1:741852/1‑150 (MQ=255)          aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                           >  1:126874/1‑150 (MQ=255)           tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                          >  2:849293/1‑150 (MQ=255)             tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                        <  1:793304/150‑1 (MQ=255)              aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                      <  2:614006/151‑1 (MQ=255)               aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                     <  2:252498/151‑1 (MQ=255)                    tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                 <  2:571825/150‑1 (MQ=255)                    tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                 >  1:766313/1‑150 (MQ=255)                     gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAt                                                                                                                               >  2:920595/1‑151 (MQ=255)                     gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAt                                                                                                                               >  2:307109/1‑151 (MQ=255)                       ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                              >  1:561327/1‑150 (MQ=255)                        cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                             >  2:518354/1‑150 (MQ=255)                         aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                             >  1:897930/1‑149 (MQ=255)                             ggcccttttttttCTGCTCTTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTCCCAAAACAGTTTAGATAATCTTTTTGTCGCTAAATGCCTGAAAAACTTGCATGGGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGtt                                                                                                                         <  1:705439/142‑1 (MQ=255)                                ccACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtt                                                                                                                    >  1:967969/1‑151 (MQ=255)                                    ttttttCTGCTCTTTTTTTTTTCATGCAT‑TTTTCTTCCTGCGAAATCTTTTCCAAAAACAGTTTAAATAATTTTTTTGTCGTTAAATCCCTAAAAAACTTGAATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                 <  1:566492/150‑1 (MQ=255)                                    ttttttCTGCTCGTTTTTTTTTCA‑GCAATTTTTTTTGCGGCGCAAGCGTTTCCAAAAACAGTTTAAATGATTTTTTTGCCGTTTAAGCCCTTAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                 <  1:40183/150‑1 (MQ=255)                                    ttttttCTGCTCGTTTTTATTTCATGAAT‑TCTTCTTGCTGCGAAAGCTTTTCCCAAAACAGGTTAAATGATTTTTTTGTCGTTTAATCCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                 <  1:819661/150‑1 (MQ=255)                                      ttttCTGCTCTTTTTTATTTCATGCAT‑TCTTCTTGCTGCGCAAGCTTTTTCCAGAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                               <  2:505460/150‑1 (MQ=255)                                            gCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                         >  2:728283/1‑150 (MQ=255)                                             cTCGTTTTTATTTCATGCAT‑TCTTCTTGCGCCGAAAGCGTTTCCAAAACCGGGTTAAATAATTTTTTTGCCGCTTAATGCGTGTAAACCATAAGTAAGCCAAAAAAAAATATAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                         <  2:658178/149‑1 (MQ=255)                                              tttttttttttTCATGCAT‑TCTTTTTGCTGCGCAAGCTTTTCCCAAAACAGGTTAGATGATTTTTTTGTCGTTTAATCCCTGTAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                         <  1:977798/145‑1 (MQ=255)                                              tCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                      >  1:922934/1‑151 (MQ=255)                                              tCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                      >  1:833923/1‑151 (MQ=255)                                               cGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                     >  2:106580/1‑151 (MQ=255)                                                 tttttATTTGAGGAAT‑TCTTCTTGCGGCGGAGGTGTTTTCAAGAAGGGGTTAGATGATTTTTTGGTCGCGTAGTGCGTGAAAAATATGCGTGGGCCAAAAAATAATAAAAAAAGTCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                    <  1:610952/150‑1 (MQ=255)                                                  ttttATTCCAAGCAC‑TCTTCTTCCCTCGCAGGTGTTTTCGATAACAGGTTAGATTTTCTTTTTGTCCGTTAATGCCTGAAAAACTGTTGTGGCCCAAAAAAAAATATAAAAAGTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                   <  2:335125/150‑1 (MQ=255)                                                     tttttCATTCAT‑TCTTTTGGCTGCGCAAGCTTTTTCAAAAACAGTTAAGATGATTTTTTTGTCACTAAATCCCTAAAAAACATGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                <  1:891701/148‑1 (MQ=255)                                                     tATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                >  1:734220/1‑150 (MQ=255)                                                         tCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                           >  2:467693/1‑151 (MQ=255)                                                          cATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                          >  2:672524/1‑151 (MQ=255)                                                              cAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGt                                                                                      >  1:842491/1‑151 (MQ=255)                                                                   cttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                  >  1:278764/1‑151 (MQ=255)                                                                     tcttGCTGCGCAAGCTTTTCCCAAAACAGGTTAGATGTTTTTTTTGTCGCTTAATGCTTTAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                 <  2:938973/150‑1 (MQ=255)                                                                       ttGCTGCGCAAGCTTTTTCCAGAACGGGTTAGATGATTTTTTTGTCGCTTAATGCCTGTAAAACATGAATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                <  2:443958/149‑1 (MQ=255)                                                                       ttGCTGCGAAAGCGTTTCCCAGAACTGGTTATATGATTTTTTTTTCTCTAAATCCCTGAAAAACATGCATGCCCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                               <  2:486638/150‑1 (MQ=255)                                                                         gctgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                             >  2:208554/1‑150 (MQ=255)                                                                             cgcAAGCTTTTCCCAAAAAAGGTAAGATGATTTTTTTGTCGTTTAAGCCCTGAAAAACATGAATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                           <  2:456342/148‑1 (MQ=255)                                                                             cgaaggcttttaCCAAAACAGGTTAAATAATTTTTTTGTCGATTATTGCCTGAAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                         <  1:669413/142‑1 (MQ=255)                                                                                         ccAGAACAGTTTAGATGATTTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                              <  2:220406/149‑1 (MQ=255)                                                                                           aaaaCAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGGGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                             <  2:760342/146‑1 (MQ=255)                                                                                             aaCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                        >  2:702080/1‑151 (MQ=255)                                                                                              aCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                       >  1:217504/1‑151 (MQ=255)                                                                                              aCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                       >  2:387127/1‑151 (MQ=255)                                                                                              aCAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                       >  1:826243/1‑151 (MQ=255)                                                                                              aCAGGTAAGATGATTTTTTTGTCGCTTAATGCCTAAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                       <  2:228256/151‑1 (MQ=255)                                                                                                 gtttAGATGATTTTTTTGTCGCTTAATGCTTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                      <  2:213888/147‑1 (MQ=255)                                                                                                   ttAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                   >  1:325580/1‑150 (MQ=255)                                                                                                       atgatTTTTTTGTCATTAAACCCTTGTAAAACGTGAATGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                 <  2:747379/148‑1 (MQ=255)                                                                                                       atgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                               >  1:560373/1‑150 (MQ=255)                                                                                                         gatCTTTTTGTCCCTTAATGCCTAAAAAACATTCATGCGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                              <  1:849168/149‑1 (MQ=255)                                                                                                           tttttttGTCGTTTAATCCCTATAAAACATGCATGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                <  1:242315/143‑1 (MQ=255)                                                                                                               tttGTTGCTAAATGCTTAAAAAACATGCTTGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                       <  2:734331/150‑1 (MQ=255)                                                                                                                       ttAATGCTTAAAAAACATGAAGGACCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                 <  2:923070/148‑1 (MQ=255)                                                                                                                        tAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGTGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                             >  1:298069/1‑151 (MQ=255)                                                                                                                         aatctctaaaaaaCATCCATGAGCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                <  2:498556/138‑1 (MQ=255)                                                                                                                         aaTGCCTGAAAACCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                            <  2:898062/151‑1 (MQ=255)                                                                                                                          aTGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                            >  1:308293/1‑150 (MQ=255)                                                                                                                          aTGCCTAAAAAACAGGAATGAGCCACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                           <  1:796847/151‑1 (MQ=255)                                                                                                                              cTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                       >  2:241988/1‑151 (MQ=255)                                                                                                                                  aaaaCATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                       <  2:362017/147‑1 (MQ=255)                                                                                                                                  aaaaCATGCATGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                     <  2:697634/149‑1 (MQ=255)                                                                                                                                  aaaaCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                   <  2:315532/151‑1 (MQ=255)                                                                                                                                    accatgcatgGGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                 <  2:219892/149‑1 (MQ=255)                                                                                                                                    aaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                  >  1:219836/1‑150 (MQ=255)                                                                                                                                       atgcatgACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag              <  1:307045/151‑1 (MQ=255)                                                                                                                                        tgcatgACCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATGTCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa               <  2:750571/149‑1 (MQ=255)                                                                                                                                                ccAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa      <  2:888579/150‑1 (MQ=255)                                                                                                                                                 cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat     >  2:173287/1‑150 (MQ=255)                                                                                                                                                   caAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt  >  1:858707/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                       GTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAA‑TCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATT  >  NC_000913/29341‑29636

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