Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,496	0	T	A	51.6%	90.6 / 159.4	93	intergenic (+301/‑155)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (37/8);  new base A (0/48);  total (37/56)
Fisher's exact test for biased strand distribution p-value = 1.82e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-16

TTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29346‑29646                                                                                                                                                         |                                                                                                                                                       tttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                        >  2:374374/1‑151 (MQ=255)   tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                       <  2:1486933/150‑1 (MQ=255)   tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                       <  2:45588/150‑1 (MQ=255)    aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATaa                                                                                                                                                      <  2:133622/150‑1 (MQ=255)       gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                                   <  2:425129/150‑1 (MQ=255)       gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAACGCCTG‑TAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                                  <  2:70262/151‑1 (MQ=255)         aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                                 <  1:85743/150‑1 (MQ=255)            ttttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                             >  2:186463/1‑151 (MQ=255)             tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCg                                                                                                                                             <  2:1174480/150‑1 (MQ=255)               tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           <  2:64292/150‑1 (MQ=255)                  ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                        >  1:63489/1‑150 (MQ=255)                         gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCAAGAAAAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAACCATGCATGAGCAACAAATTAATATAAAAAATCCCGCCATTAAGTTGat                                                                                                                                >  1:675922/1‑150 (MQ=255)                             aCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                            >  2:669630/1‑151 (MQ=255)                              actttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                            >  1:1155175/3‑150 (MQ=255)                               ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                            >  2:1151163/1‑149 (MQ=255)                                 ttttCTGCTCTTTTTTTTTTCATACATTCTTTTTGCTGCGAAAGCTTTTTCAAAAACAGTTAAAATATTTTTTTTG‑TCTTTAAATCCCTA‑AAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAgc                                                                                                                         <  1:535723/150‑1 (MQ=255)                                 ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                        >  1:976037/1‑151 (MQ=255)                                   ttCTGCTCTTTTTTATTTCATGCAATCTTTTTGCTGCGCAAGCTTTTCCCAGAACAGGTTAGATGATTTTTTTG‑TCTTTAAAGGCCTG‑AAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                                       <  2:1093796/150‑1 (MQ=255)                                      tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                   >  1:556589/1‑151 (MQ=255)                                       gcttttttttATTTCATGCAATCTTCTTGCTGCGAAAGCTTTTCCCAGAAGAGGTTAGATGTTCTTTTTG‑TCGTTAAATGCCTG‑AAAAACATGCAGGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                   <  2:704727/145‑1 (MQ=255)                                       cctctTTTTTATTTCATGCATTCTTCTTGCTGCGCAAGCTTTTTCCAAAACAGGTTAGATGTTTTTTTTG‑TCGTTAAATGCCTT‑TAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                   <  2:467667/149‑1 (MQ=255)                                                tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                          >  1:1412305/1‑150 (MQ=255)                                                tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                          >  2:153442/1‑150 (MQ=255)                                                   ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                      >  1:1163428/1‑151 (MQ=255)                                                     cATGCAATCTTCTTCCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                                    >  1:1459870/1‑151 (MQ=255)                                                      aTGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                                   >  1:1396997/1‑151 (MQ=255)                                                          attttttttGCTGCG‑AAGGCTTTTCCAAGAAAGGTTAAAATATTTTTTTGCTCTTTTAATCCTTG‑AAAAACAGTAATGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                               <  2:837291/144‑1 (MQ=255)                                                          aaTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                               >  2:1007451/1‑151 (MQ=255)                                                              ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  2:880538/1‑151 (MQ=255)                                                               tcttGCCGCACAAGCTTTTTCAAAAACAGTTTAGATTTTTTTTTTG‑TCTCTTAATCCCTA‑AAAAACATGCAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                             <  1:468657/148‑1 (MQ=255)                                                                tttGCTGCGCAAGCTTTTCCCAGAACAGTTTAGATGATTTTTTTG‑TCGTTAAATGCCTG‑AAAAACTTGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                          <  2:143753/149‑1 (MQ=255)                                                                          aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                               >  1:1596259/1‑151 (MQ=255)                                                                           aGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                              >  1:927496/1‑151 (MQ=255)                                                                             cGTTTTCCAAAACAGGTTAGATGATTTTTTTG‑TCGTTTAATGCCTG‑AAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            <  2:1200038/151‑1 (MQ=255)                                                                             cGTTTTCAAAAACAGTTTAGATGATTTTTTTG‑TCGTTTAATCCCTG‑AAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            <  2:285332/151‑1 (MQ=255)                                                                              tttttCAAGAAAAGTTAAGATATTTTTTTTG‑TCGCTTAATGCCTA‑AAAAACATGCATGACCCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            <  1:763584/149‑1 (MQ=255)                                                                              ttttcCCAGAACGGTTTAGATAATTTTTTTG‑TCTTTTAATCCCTA‑AAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            <  1:425153/149‑1 (MQ=255)                                                                                 tcccAGAACAGTTTAGATTATTTTTTTG‑TCTCTTAAT‑CCTGATAAAACATGCATGACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                         <  1:289321/148‑1 (MQ=255)                                                                                  cccAGAACAGTTTAGATGTTTTTTTTG‑TCGCTTAATGCCTG‑AAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                        <  2:1205455/149‑1 (MQ=255)                                                                                  ccaaaaagggtttAGATGATTTTTTTG‑TCGCTAAATGCGTA‑AAAAACATGCATGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                       <  2:313460/140‑1 (MQ=255)                                                                                  ccaaGAACGGGTTAGATGATTTTTTTG‑TCGTTTAAGCCTTG‑AAAAACATGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                        <  2:721338/147‑1 (MQ=255)                                                                                   ccAGAACAGTTTAGATGATTTTTTTG‑TCGTTTAATGCCTG‑AAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                       <  2:1520201/150‑1 (MQ=255)                                                                                     aGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    >  1:1075950/1‑151 (MQ=255)                                                                                       aaaagtttAGATGTTTTTTTTG‑TCGCTTAATGCCTG‑AAAAACATGCATGACCAAAAAAATAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                                  <  1:79434/145‑1 (MQ=255)                                                                                        aCAGGTTAGATGTTTTTTTTT‑TCGTTTAATGCCTG‑AAAACCATGCATGAGCCACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                                 <  2:1418533/151‑1 (MQ=255)                                                                                         agggTTAGATGTTTTTTTTG‑TGGTTAAAGGCCTG‑AAAAACATGAATGAGCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                                <  1:629901/149‑1 (MQ=255)                                                                                          ggtttAGATGATTTTTTTG‑TCGTTTAATGCTTA‑AAAAACAGGCATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                               <  2:955941/148‑1 (MQ=255)                                                                                          ggtttAGATGATTTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGAATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                               <  2:1398317/148‑1 (MQ=255)                                                                                          aGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAt                                                               >  2:623047/1‑151 (MQ=255)                                                                                                  tgatCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                       >  2:292594/1‑151 (MQ=255)                                                                                                   gatCTTTTTG‑TCGCTTAATGCCTG‑AAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                       <  2:1534313/150‑1 (MQ=255)                                                                                                   attttttttt‑TCGTTTAATGCTTA‑AAAAACATGCATGACCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                      <  2:1179939/147‑1 (MQ=255)                                                                                                           tG‑TCGCTAAATGCCTG‑AAAAACATGAATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                               <  2:185070/150‑1 (MQ=255)                                                                                                            g‑TCGTTTAATGCCTG‑TAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                             <  2:775695/151‑1 (MQ=255)                                                                                                              tcctttAATGCCTA‑AAAAACAGGCATGACCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                             <  1:1111741/146‑1 (MQ=255)                                                                                                               cgtttAATCCCTG‑AAAAACATACATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                           <  1:1149086/148‑1 (MQ=255)                                                                                                                cgtttattgcttaaaaaaCATGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            <  1:926528/136‑1 (MQ=255)                                                                                                                 cTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                          >  1:833642/1‑150 (MQ=255)                                                                                                                 cTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                          >  2:1567845/1‑150 (MQ=255)                                                                                                                  ttAATGCCTA‑AAAAAAATGCATGACCCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGaa                                         <  1:490188/150‑1 (MQ=255)                                                                                                                   tAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                       >  1:907200/1‑151 (MQ=255)                                                                                                                   tAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                       >  2:1321957/1‑151 (MQ=255)                                                                                                                   tAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                       >  2:448794/1‑151 (MQ=255)                                                                                                                    aaTGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                       >  2:1060827/1‑150 (MQ=255)                                                                                                                      tGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                    >  1:1103409/1‑151 (MQ=255)                                                                                                                       gCCTG‑TAAAACATGCATGAGCCACAAAGTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                    >  2:1023830/1‑150 (MQ=255)                                                                                                                        gcctaaaaaaCATGCATGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                    <  1:1456882/144‑1 (MQ=255)                                                                                                                         cctaaaaaaCATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                  <  1:473525/146‑1 (MQ=255)                                                                                                                         cctaaaaaaCATGCAGGCCCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                                  <  1:1243882/146‑1 (MQ=255)                                                                                                                            aaaaaaCATGCATGACCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                                <  2:1037549/148‑1 (MQ=255)                                                                                                                            aaaaaaCATGCATGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                               <  1:1430260/149‑1 (MQ=255)                                                                                                                             tAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                               >  1:483315/1‑150 (MQ=255)                                                                                                                                 aCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                          <  2:639102/151‑1 (MQ=255)                                                                                                                                  aaggcatgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                         <  1:313172/148‑1 (MQ=255)                                                                                                                                    ggcatgACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  1:1469813/150‑1 (MQ=255)                                                                                                                                     gcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  2:1144319/1‑151 (MQ=255)                                                                                                                                     gcatgACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  1:374356/150‑1 (MQ=255)                                                                                                                                     gcatgAACCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  1:1073949/150‑1 (MQ=255)                                                                                                                                     gcaggaccaaaaaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  1:72173/139‑1 (MQ=255)                                                                                                                                      catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  2:663362/1‑150 (MQ=255)                                                                                                                                      aatgaccaaaaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtga                     <  1:33125/141‑1 (MQ=255)                                                                                                                                       atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  1:445800/1‑149 (MQ=255)                                                                                                                                       atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  1:788469/1‑149 (MQ=255)                                                                                                                                         gAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  <  1:292578/151‑1 (MQ=255)                                                                                                                                          agcaaaaaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga                 <  1:381081/145‑1 (MQ=255)                                                                                                                                          aGCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga                 <  2:1479274/151‑1 (MQ=255)                                                                                                                                            ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               >  1:1499315/1‑151 (MQ=255)                                                                                                                                            ccAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               <  1:1419651/151‑1 (MQ=255)                                                                                                                                             cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               >  2:894362/1‑150 (MQ=255)                                                                                                                                             cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat               >  2:299938/1‑150 (MQ=255)                                                                                                                                             caaaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata              <  2:1044153/148‑1 (MQ=255)                                                                                                                                              acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat             >  2:1228230/1‑151 (MQ=255)                                                                                                                                               aaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt            <  1:473074/150‑1 (MQ=255)                                                                                                                                                  aataatAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc         <  1:405071/151‑1 (MQ=255)                                                                                                                                                   ataataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct        <  2:715868/151‑1 (MQ=255)                                                                                                                                                     aataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      <  2:1183932/151‑1 (MQ=255)                                                                                                                                                       tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  1:1454214/1‑149 (MQ=255)                                                                                                                                                       tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  2:874603/1‑149 (MQ=255)                                                                                                                                                         taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  2:1577988/1‑151 (MQ=255)                                                                                                                                                         taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  1:1435841/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       TTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTG‑TCGCTTAATGCCTG‑TAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29346‑29646

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