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

ATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29350‑29646                                                                                                                                                      |                                                                                                                                                       aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                     >  1:607076/1‑150 (MQ=255) aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                     >  2:568536/1‑150 (MQ=255)  tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                                    >  2:23220/1‑150 (MQ=255)    tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                                  <  2:805173/150‑1 (MQ=255)    tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                                 >  1:56576/1‑151 (MQ=255)     aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                                 <  1:312474/150‑1 (MQ=255)     aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                                >  1:136879/1‑151 (MQ=255)          ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           >  1:704319/1‑151 (MQ=255)           tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           >  1:373771/1‑150 (MQ=255)           tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                           >  1:469746/1‑150 (MQ=255)               cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                       >  2:110/1‑150 (MQ=255)                aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                       >  1:302015/1‑149 (MQ=255)                 tttGTTCCTTTTTTTTCTGCTCGTTTTTATTCAATCCATTCTTCTTGCTGCGCAGGCTTTTCCCAAAACAGGTT‑AGATAATTTTTTTTTCTTTTAATCCC‑TAA‑AAAACATGCATGCGCCAAAAAAAAATATAAAAAACCCCGCCATTAAg                                                                                                                                     <  2:739451/150‑1 (MQ=255)                  ttGGTCCTCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCTTTTTCAAAAACAGTTT‑AGATAATTTTTTTGTCGCTAAATGCC‑TGA‑AAAACATGCATGAGCCACAAAATAATATAAAAAACCCCGCCATTAAGt                                                                                                                                    <  2:224001/150‑1 (MQ=255)                   tGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                                  >  1:745994/1‑151 (MQ=255)                      tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                               >  1:526123/1‑151 (MQ=255)                        cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                             >  1:672491/1‑151 (MQ=255)                             ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                                        >  2:828440/1‑151 (MQ=255)                                      cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                               >  2:394327/1‑151 (MQ=255)                                       gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                              >  2:503641/1‑151 (MQ=255)                                        ttttttttTCATGCAATTTTCTTCTTGCGAAGGTTTTTTCAAAACCAGTTA‑AAATGTTTTTTTTGTCCTTTAATGCT‑TTT‑AAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                                 <  1:640971/147‑1 (MQ=255)                                        tttttATTTCATGAAATCTTCTGGCTGCGCAAGCTTTTCCCAGAACGGTTT‑AGATGATTTTTTTGTCGCTAAATGCC‑TGA‑AAAACTTGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                                               <  1:720134/149‑1 (MQ=255)                                          tttATTTCATGAATATTTTTTGCTCGGCAATGTTTTTCAAGACCAGTTT‑AAATGAT‑TTTTTGTCCGTTAAATCC‑TTTAAAAAAATGTATGCCCCAAAAAATAATATAAAACACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                            <  2:685976/150‑1 (MQ=255)                                            tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                          >  1:22882/1‑150 (MQ=255)                                             aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAATTTTAGCGCCCATATCTCCAGAAt                                                                                                        >  2:566989/1‑151 (MQ=255)                                              tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       >  2:85110/1‑151 (MQ=255)                                              tttCAGGCATTCTTCTTGCGGCGCAGGCTTTTCCCAGAAAAGGTA‑AGATGATTTTTTTGTCTTTAAAGCCC‑GGA‑AAAACAGGCAGGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                        <  1:764687/150‑1 (MQ=255)                                                  aTGCAATCTTTTTGCGGCGCAAGCTTTTTCAAGACAAGGTT‑AAATAATTTTTTTTTCGTTTAATCCC‑TAA‑AAAACATGCTTGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                                    <  2:601395/150‑1 (MQ=255)                                                    gCATTCTTCTTGCTGCACAAGCTTTTTCAAGAACAGGTT‑AAATGATTTTTTTGTCTCTTAATGCC‑TGT‑AAAACATTCATGAGCAACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       <  2:89663/145‑1 (MQ=255)                                                      attttccttGCTGCGAAAGCTTTTCCAAAAACAGTTA‑AGATGATTTTTTTGTCGCTAAATCCC‑GGA‑AAAAAATGCATGACCCACAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                                     <  2:239023/139‑1 (MQ=255)                                                      attcttctTGCTGCGCAAGCGTTTACAAAAACAGTTT‑AGATAATTTTTTGGTCGTTTAACGCC‑TGA‑AAAAAATTAAGGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                                  <  2:49172/146‑1 (MQ=255)                                                             ttGCTGCGCAAGCGTTTTCAAAAAAAGGTT‑AAATAATCTTTTTGTCTCTTAATGCC‑TTA‑AAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                                         <  2:83029/150‑1 (MQ=255)                                                             ttGCTGCACATGTTTTTTCCAGAGCAGTTA‑AGATTATTTTTTTTTCGCTAAATCCC‑TAT‑AAAACAGGCATGCGCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                          <  2:68209/149‑1 (MQ=255)                                                                  gcaaaggctttttCAAGAACGGTTA‑AGATGATTTTTTTGTCGTTTAATGCC‑TAA‑AAAACATGCATGACCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                       <  2:577322/138‑1 (MQ=255)                                                                   cgaaacctttttCAAAAACGGTTT‑AGATGATTTTTTTGCCTTTAAATCCT‑TAA‑AAAACATTCATGAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                     <  2:314747/140‑1 (MQ=255)                                                                      aggctttttCAAAAACGGGTT‑AGATGTTTTTTTTGCCG‑TTAATTCCTTGA‑AAAACAGGTGTGACAAACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                                <  1:126702/145‑1 (MQ=255)                                                                      aagttttttCCAAAACGGGTT‑AAATGATCTTTTTTTCGTTAATTGCA‑TTA‑AAAACATGCATGTCCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                                <  2:221437/145‑1 (MQ=255)                                                                           ttttCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            >  1:652513/1‑149 (MQ=255)                                                                            tttCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                          >  2:48550/1‑150 (MQ=255)                                                                               caaGAACAGGTT‑AAATGATTTTTTTGTCTTTTAATGCT‑TGA‑AAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                      <  1:474660/149‑1 (MQ=255)                                                                                 aGAACAGTTT‑AGATGATTTTTTTGTCGTTTAATGCC‑TGA‑AAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    <  1:220789/151‑1 (MQ=255)                                                                                 aGAACAGTTT‑AGATGATTTTTTTGTCGCTTAATGCT‑TGA‑AAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                                     <  1:820699/150‑1 (MQ=255)                                                                                  aaaaaGGTTAAGATGATTTTTTTTTCTTTTATTCCC‑TGA‑AAAAATTGCTTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                      <  2:566108/145‑1 (MQ=255)                                                                                    aCAGTTT‑AGATCATCTTTTTTTCGATTAATGCC‑TGA‑AAAACATGCATGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                        <  2:136842/144‑1 (MQ=255)                                                                                            aGATTATTTTTTTTTCTCTTACTGCC‑TAA‑AAAACATGCATGAACCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                              <  2:114054/147‑1 (MQ=255)                                                                                            aGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGCTTTGCATTGAtt                                                          >  2:798088/1‑151 (MQ=255)                                                                                               tgttttttttGTCTTTTAATGCC‑TAA‑AAAACATGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                        <  1:742625/145‑1 (MQ=255)                                                                                               tgatTTTTTTGTCGTTTAATGCC‑TGA‑AAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                        <  1:188717/150‑1 (MQ=255)                                                                                               tgatCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                       >  1:120667/1‑151 (MQ=255)                                                                                                 atCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                     >  1:786668/1‑151 (MQ=255)                                                                                                  tCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                    >  2:457018/1‑151 (MQ=255)                                                                                                    ttttttTCTCTTAATCCC‑TAA‑AAAAAATTAATGACCAAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                         <  2:591716/144‑1 (MQ=255)                                                                                                       ttGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                               >  2:736392/1‑151 (MQ=255)                                                                                                       ttGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                               >  1:438929/1‑151 (MQ=255)                                                                                                        tGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                              >  1:87830/1‑151 (MQ=255)                                                                                                          tCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            >  1:257190/1‑151 (MQ=255)                                                                                                               tAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                       >  2:434008/1‑151 (MQ=255)                                                                                                                       ctaaaaaaCATTGATGACCCAAAAAAAAATAAAAAAAACCCCCGCATTTAATTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                    <  2:704096/144‑1 (MQ=255)                                                                                                                           aaacCATGCATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                                   <  1:538588/145‑1 (MQ=255)                                                                                                                           aaaaCTTGCATGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                                 <  1:85139/147‑1 (MQ=255)                                                                                                                           aaaaCATGAATGGCCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                             <  1:200710/151‑1 (MQ=255)                                                                                                                            aaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                            >  1:650173/1‑151 (MQ=255)                                                                                                                             aaCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                           <  1:799452/151‑1 (MQ=255)                                                                                                                                atgcatgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                         <  2:541862/150‑1 (MQ=255)                                                                                                                                  gaatgAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  1:182175/148‑1 (MQ=255)                                                                                                                                   cttgAGCAAAAAAAAAATAAAAAAAACGCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                       <  2:501314/147‑1 (MQ=255)                                                                                                                                   catgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                         <  1:303876/147‑1 (MQ=255)                                                                                                                                    attggCCAAAAAAAAAAAAAAAAAACCCCGCCCTTAAGTTGGCTTTTAGCGGCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                         <  1:791009/142‑1 (MQ=255)                                                                                                                                    atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  1:641531/1‑149 (MQ=255)                                                                                                                                         ccAAAAAATAAAATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  <  2:25364/148‑1 (MQ=255)                                                                                                                                             aaaaaaATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt            <  1:617579/150‑1 (MQ=255)                                                                                                                                             aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc           >  1:341358/1‑151 (MQ=255)                                                                                                                                             aaaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc           <  2:56566/151‑1 (MQ=255)                                                                                                                                                ataaattaAAACAGCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAGTGATTTACGTCATCATTGTGAATTAATATGCTAATAATGTGAGTGAATATtc           <  2:655955/142‑1 (MQ=255)                                                                                                                                                aaaataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct        <  1:659504/149‑1 (MQ=255)                                                                                                                                                    tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  1:446210/1‑149 (MQ=255)                                                                                                                                                      taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  1:721276/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       ATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCC‑TGT‑AAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29350‑29646

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