Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,496	0	T	A	36.9%	143.7 / 36.5	84	intergenic (+301/‑155)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (38/15);  new base A (0/31);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.24e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646                                                                                                                                                        |                                                                                                                                                       tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑tata                                                                                                                                                       <  1:1039308/150‑1 (MQ=255)  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATaa                                                                                                                                                      >  1:493860/1‑150 (MQ=255)   aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATaaa                                                                                                                                                     >  1:1046917/1‑150 (MQ=255)    tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATaaaa                                                                                                                                                    >  2:1097388/1‑150 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATc                                                                                                                                                <  2:162816/151‑1 (MQ=255)         aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATccc                                                                                                                                              <  2:472432/151‑1 (MQ=255)           tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTTTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCg                                                                                                                                             <  1:1229016/150‑1 (MQ=255)           tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGc                                                                                                                                            >  2:581668/1‑151 (MQ=255)            ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGcc                                                                                                                                           >  1:652973/1‑151 (MQ=255)            ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGc                                                                                                                                            <  2:1147311/150‑1 (MQ=255)            ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGc                                                                                                                                            <  1:246681/150‑1 (MQ=255)             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGcc                                                                                                                                           >  2:411506/1‑150 (MQ=255)             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGcc                                                                                                                                           <  1:846288/150‑1 (MQ=255)              gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCa                                                                                                                                          <  1:1372378/150‑1 (MQ=255)              gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCa                                                                                                                                          <  1:486125/150‑1 (MQ=255)                    ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GGGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGtt                                                                                                                                   >  2:1018636/1‑151 (MQ=255)                      gttccttttttttCTGCTCTTTTTTTTTCAATGCAATTTTCTTGCGGCGCAAGCGTTTTCAAG‑AACAGTTTAGATGTTCTTTTTGTCGTTAAATGCTTGAAAAACATGCAG‑GAGCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGtt                                                                                                                                   <  2:737052/142‑1 (MQ=255)                      ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTg                                                                                                                                  >  2:1336415/1‑150 (MQ=255)                        tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                               >  2:1300466/1‑151 (MQ=255)                          cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                             >  2:1298767/1‑151 (MQ=255)                           aCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                            >  2:744742/1‑151 (MQ=255)                             ttttttCTGCTCTTTTTTTTTTCATGCATTTTTCTTGCTGCGCAACCTTTTCCCAA‑AACAGTTTAGATGATTTTTTTGTCTTTAAATGCTTGAAAAACATGCAG‑GACCCAAAAAAAAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                           <  2:29027/150‑1 (MQ=255)                              tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                         >  2:172502/1‑151 (MQ=255)                                  tCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                     >  1:513692/1‑151 (MQ=255)                                     gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                   >  2:192880/1‑150 (MQ=255)                                          tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                             >  2:367793/1‑151 (MQ=255)                                            tttatttaATGCAATCTTCTGTCTGAGCAACCTTTTTAAAC‑AACATTTTAGATGATTTTTTTGTCCTTTAATGCCGGAAAAATTTGTAT‑GGCCCAAAAAATAA‑TAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                            <  1:1305601/150‑1 (MQ=255)                                            tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                             >  1:245145/1‑149 (MQ=255)                                            tttATTTCATGCAATCTTCTTGCGGCGCAAGCTTTTCCAAA‑AACAGGTAAGATGATTTTTTTGTCTCTTAATGCCTGAAAAACATGCAT‑GACCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                            <  2:1248115/150‑1 (MQ=255)                                                tttCATGCATTCTTCTTGCTGCGCAAGCTTTTTCCAG‑AACAGGTTAAATAATTTTTTTGTCGCTTAATCCCTGTAAAACTTGCAT‑GACCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                          <  1:3351/148‑1 (MQ=255)                                                tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       >  2:1374073/1‑151 (MQ=255)                                                tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       >  1:994299/1‑151 (MQ=255)                                                            ttcttGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  1:1110610/1‑151 (MQ=255)                                                            ttcttGCTGCGCAAGCGTTTTCCAA‑AACAGGTTAGATGATTTTTTTTTCGCTTAATCCCTAAAAAACATGCAT‑GACCCACAAAATAA‑TAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                              <  1:1122296/148‑1 (MQ=255)                                                             tcttGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  2:672761/1‑150 (MQ=255)                                                               ttgcgccgcAAGCTTTTACCAA‑AAAAGTTTAGATAAGTTTTTTGTCG‑TTAAAGCCTGGAAAAACTTCTTGGACCCAAAAAATAA‑AAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                                 <  1:857387/120‑1 (MQ=255)                                                                 gctgcGCAAGCTTTTTCAAG‑AACAGGTTAAATGATTTTTTTGTCGTTTAATGCTTAAAAAACATGCAT‑GACCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                                         <  2:362151/148‑1 (MQ=255)                                                                   tgcGAAAGCGTTTCCAAA‑AACAGGTTAGATGATTTTTTTGTCGCTTAATGCTTGAAAAACATGCAT‑GACCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                       <  2:579518/148‑1 (MQ=255)                                                                    gcgcAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                                   >  1:455290/1‑151 (MQ=255)                                                                    gcgcAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                                   >  1:1012407/1‑151 (MQ=255)                                                                      gaaggCGTTTCCAAG‑AACAGTTTAGATGATTTTTTTGTCTCTTATTGCCTTAAAAACATGAAG‑GACCCAAAAAATAA‑TAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                    <  1:647827/144‑1 (MQ=255)                                                                         aggtttttCCAAAAACAGTTAAGTTGTTTTTTTTGCCGCTAAGCGCTTAAAAAACATCAAT‑AACCCA‑AAAATAATTATAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAGtc                                                                                   <  1:979314/143‑1 (MQ=255)                                                                         aGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                              >  2:911288/1‑151 (MQ=255)                                                                          gttttttCAAA‑AACAGGTTGTATTATTTTTTTGTCACTTAATGCCTTAAAAACATGAAT‑GAGCCACAAAAAAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                             <  1:70929/148‑1 (MQ=255)                                                                            gTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            >  1:1093575/1‑150 (MQ=255)                                                                                tCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                         >  2:1361550/1‑149 (MQ=255)                                                                                 ccAG‑AACGGGTTAGATGATCTTTTTTTCGTTTATTCCCTGTAAAAAATGCAT‑GCGCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                       <  1:729373/150‑1 (MQ=255)                                                                                   aG‑AACAGTTTAGATGATCTTTTTGTCGTTAAATGCCTGAAAAACAGGCAT‑GACCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    <  2:1191053/151‑1 (MQ=255)                                                                                    g‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                                   >  1:1380676/1‑151 (MQ=255)                                                                                    g‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                                   >  1:82822/1‑151 (MQ=255)                                                                                    aaaaCAGTTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGCAT‑GAGCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                       <  2:949043/146‑1 (MQ=255)                                                                                      aaCAGATTAAATCTTTTTTTTGTCGTTTATTGCATAAAAAACATGCAG‑GAGCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    <  2:642742/149‑1 (MQ=255)                                                                                        gaGGTAATATGATTTTTTTGTCCGTAATTGCCTGAAAAACGGGGCT‑GATCAAAAAAATAT‑TAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                                  <  2:737874/148‑1 (MQ=255)                                                                                            ttAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATGCAT‑GAGCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                             <  2:872063/150‑1 (MQ=255)                                                                                            ttAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                             >  2:1222631/1‑150 (MQ=255)                                                                                            ttAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                             >  1:999484/1‑150 (MQ=255)                                                                                                 tgttttttttGTCGTTTAATGCCTGAAAACCATGCAT‑GAGCCAAAAAAAAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                           <  2:203746/142‑1 (MQ=255)                                                                                                 tgatTTTTTTGTCGCTTATGCCCTGAAAACCATACAT‑GTGCCAAAAAATAA‑TAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                             <  1:919979/145‑1 (MQ=255)                                                                                                  gatCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                      >  2:414105/1‑151 (MQ=255)                                                                                                    tttttttGTCGTTTAATCCCTGTAAAACATGCAT‑GAACCAAAAAATAA‑TATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                        <  1:1291386/145‑1 (MQ=255)                                                                                                        tttGTCGTTAAATGCCTGAAAAACATGCAG‑GAGCCAAAAAAAAA‑TAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                                 <  2:1283176/150‑1 (MQ=255)                                                                                                            tCGTTTAATGCCTGAAAAACATGAAT‑GGCCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                               <  2:637531/148‑1 (MQ=255)                                                                                                            tCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            >  2:973294/1‑151 (MQ=255)                                                                                                             cGCTTAATGCCTGAAAAACATGCAT‑GAGCCAAAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            <  2:48030/150‑1 (MQ=255)                                                                                                               cTTATTGCCTTTAAAGCATGCAT‑GAGCCAAAAAATAA‑TAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            <  1:358347/148‑1 (MQ=255)                                                                                                                ttATTGCCTTAAAAACTTGCAT‑GACCCACAAAAAAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                              <  1:809344/145‑1 (MQ=255)                                                                                                                 taactctcaaaaaaCATGTAG‑GACCAAAAAAATAA‑AAAAAAAACGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                             <  1:940652/135‑1 (MQ=255)                                                                                                                  attGCCTAAAAAACATGCAT‑GAGCCACAAAATAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                          <  2:980711/145‑1 (MQ=255)                                                                                                                       cTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                                 >  1:670830/1‑151 (MQ=255)                                                                                                                         gaaaaaCAGCCAT‑GAGCCAAAAAATAA‑TAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                                <  2:494008/148‑1 (MQ=255)                                                                                                                         gTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATTTGc                                >  1:87678/1‑150 (MQ=255)                                                                                                                         gTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                                >  1:1108946/1‑150 (MQ=255)                                                                                                                          tAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                               >  1:571191/1‑150 (MQ=255)                                                                                                                           aaaaCTTGAAT‑GACCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                              <  1:79998/150‑1 (MQ=255)                                                                                                                           aaaaCGTGAAT‑AGGCCAAAAAATAA‑TAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                        <  1:399561/140‑1 (MQ=255)                                                                                                                           aaaaCATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                              >  2:1115227/1‑150 (MQ=255)                                                                                                                            aactttatat‑gAGCAACAAAATAA‑TATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                            <  1:453000/143‑1 (MQ=255)                                                                                                                             aaCATGCGT‑GACCCAAAAAAAAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                           <  2:644744/151‑1 (MQ=255)                                                                                                                                  gcat‑gAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  2:988475/1‑151 (MQ=255)                                                                                                                                     t‑gAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt                   >  1:1201221/1‑151 (MQ=255)                                                                                                                                          ccAAAAAATAA‑TATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  <  2:1144760/148‑1 (MQ=255)                                                                                                                                          ccAAAAAAAAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      <  1:1262803/144‑1 (MQ=255)                                                                                                                                            acaAAATAA‑TAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCTATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat             <  1:673927/151‑1 (MQ=255)                                                                                                                                            aaaaaaTAA‑TAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat             <  2:1047690/149‑1 (MQ=255)                                                                                                                                               aaaaaa‑TGTAAAAGACTCCGCCATTAAGTTGACTTTTAGCGCCCAGATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTGCGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               <  1:515741/146‑1 (MQ=255)                                                                                                                                                       ataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAgg   >  2:684449/1‑151 (MQ=255)                                                                                                                                                       ataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAgg   >  2:1183812/1‑151 (MQ=255)                                                                                                                                                        taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  1:1223145/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAG‑AACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCAT‑GAGCCACAAAATAA‑TATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646

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