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

TTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29337‑29635                                                                                                                                                        |                                                                                                                                                      ttttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCcac                                                                                                                                                       >  1:392108/1‑150 (MQ=255)     gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACaaaa                                                                                                                                                   <  2:1108416/150‑1 (MQ=255)     gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACaaaa                                                                                                                                                   <  1:549516/150‑1 (MQ=255)     gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAaat                                                                                                                                                  <  1:366335/151‑1 (MQ=255)      ttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAaat                                                                                                                                                  >  1:568081/1‑150 (MQ=255)       tgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAaata                                                                                                                                                 >  2:875373/1‑150 (MQ=255)        gttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAaataa                                                                                                                                                <  2:807758/150‑1 (MQ=255)           ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                             >  1:64269/1‑150 (MQ=255)                gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCCCAAAATAATATaaaaa                                                                                                                                        >  2:828547/1‑150 (MQ=255)                gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                        >  2:213217/1‑150 (MQ=255)                   aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                    >  1:363868/1‑151 (MQ=255)                         gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                               <  1:1002325/150‑1 (MQ=255)                            cATTTGGTTCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                            >  2:852717/1‑150 (MQ=255)                            cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                            >  1:803893/1‑150 (MQ=255)                            cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                            >  1:1156628/1‑150 (MQ=255)                              tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGt                                                                                                                         >  1:263878/1‑151 (MQ=255)                                     cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                  >  1:142293/1‑151 (MQ=255)                                       cTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                 >  2:785833/1‑150 (MQ=255)                                        ttttttCTGCTCTTTTTTTTTTCATGAATTTTTTTTGCGGCGCAAGCTTTTCCA‑AGAACAGGTTAGATGATTTTTTTGTCGTTAAATGCTT‑GAAAAACATGAATGGCCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                  <  1:672324/148‑1 (MQ=255)                                         tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                              >  2:809175/1‑151 (MQ=255)                                           tttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                            >  1:435799/1‑151 (MQ=255)                                                 ttCGTTTTTATGTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATCATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCGCata                                                                                                       >  2:1152340/2‑150 (MQ=255)                                                   cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                     >  2:360158/1‑150 (MQ=255)                                                     tttttATTTCATGAAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                  >  2:532862/1‑151 (MQ=255)                                                      ttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                  >  2:901726/1‑150 (MQ=255)                                                        ttttttCATGCATTTTTCTTGCTGCGCAAGTTTTTTCA‑AAAACAGTTTAGATAATTTTTTTGTCGCTTAATCCTT‑AAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                <  2:1029917/147‑1 (MQ=255)                                                           ttccaagaactcttgtTGCTCCGAACGTTTTTCCC‑AGAAAAGTTTAGATATTCTTTTTTCCGATTAGTGCCT‑GAAAACCTTGTATGAACCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                             <  1:21182/140‑1 (MQ=255)                                                               aTGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                        >  2:1146962/1‑151 (MQ=255)                                                                 gCAATCTTTTTGCTGCACAACCGTTTTCC‑AGAACAGGTTCGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCCTGAGCCACAAAATAATATAAACAATCCCGACATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                       >  2:747342/1‑150 (MQ=255)                                                                    attcttcttaccgacAAGTGTTTTCCAAAAACGGTTTAGATTATTTTTTTGACGTATAATCCTT‑GTAAAAAATGATTGACCCAAAAAAAAATATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                        <  2:131810/133‑1 (MQ=255)                                                                    aTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                   >  2:477648/1‑151 (MQ=255)                                                                    aTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                   >  1:390581/1‑151 (MQ=255)                                                                      cttcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                 >  2:656035/1‑151 (MQ=255)                                                                              tgcACAAGCTTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGCTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                         >  2:823243/1‑151 (MQ=255)                                                                                   agggtttttCA‑AGAACAGTTTAGATGATTTTTTTGTCGTTAAATGCTT‑GAAAACCATGCAGGGGCCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                        <  1:1052160/142‑1 (MQ=255)                                                                                    aggcttttCCAAGAACAGGTTAGATTTTTTTTTTGTCGTTTAATGCCT‑GAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                              <  1:1143823/137‑1 (MQ=255)                                                                                      ctttttCA‑AGAAAAGTATAGTTAATTTTTTTGCCTC‑TAATTCCTAGTAAAAAAAGCATGACCCAGAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                    <  1:466838/146‑1 (MQ=255)                                                                                       ttttcCC‑AAAACAGTTTATATGTTTTTTTTGTCTTTTACTGCCT‑AAAAAACATGAATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                   <  2:1114669/147‑1 (MQ=255)                                                                                        ttttCCAAGAACAGGTTAGATATTTTTTTTTCCGTTAAATCCTT‑GAAAAACATGCAGGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                     <  1:981453/146‑1 (MQ=255)                                                                                         tttccaaaaaaaGGTTAGATGTTTTTTTTGTCGTTTAATCCCT‑GAAAAACATGCAGGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                   <  1:146033/136‑1 (MQ=255)                                                                                         tttCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                              >  1:471824/1‑151 (MQ=255)                                                                                               agactaGGTTAGTTTTTTTTTTTGTGGCTAAATGCAA‑ATAAAACATGCAGAAACAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                            <  2:843094/143‑1 (MQ=255)                                                                                                 aaCAGGTTAGATGTTTTTTTTGCCATTTAATGCTT‑GAAAAACTTGCAGGCGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                            <  1:815653/146‑1 (MQ=255)                                                                                                 aaCAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                       >  1:1199016/1‑151 (MQ=255)                                                                                                 aaCAGGTAAGATGTTTTTTTTGTCGTTTAATCCTT‑GAAAAACATGAATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                        <  1:971245/150‑1 (MQ=255)                                                                                                       ttaattgatTTTTTTGTCGTTAAAGCCCT‑GAAAAACTTCAATGACCCAAAAAAAAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                   <  1:654300/144‑1 (MQ=255)                                                                                                               tttttttGTCTTTTAATGCCT‑AAAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                            <  2:690311/146‑1 (MQ=255)                                                                                                               tttttttGTCGTTTAATGCCT‑GAAAAACAGGAATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                            <  1:109437/146‑1 (MQ=255)                                                                                                                  ttttGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                      >  2:1133170/1‑151 (MQ=255)                                                                                                                   tttGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                      >  1:134814/1‑150 (MQ=255)                                                                                                                        cctttAATCCCT‑AAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                   <  1:760459/145‑1 (MQ=255)                                                                                                                         gCTTATTGCCT‑GTAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                               <  2:171351/151‑1 (MQ=255)                                                                                                                             aaTGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                            >  1:896591/1‑150 (MQ=255)                                                                                                                                       aaaaCATGAATCACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGA‑TTGCATTGATTTACGTCATCATTGTGAAt                             <  2:1093292/139‑1 (MQ=255)                                                                                                                                       aaaaCAGGCATGACCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                   <  1:136852/150‑1 (MQ=255)                                                                                                                                        agacttctttgAGCCAAAAAAAAATATAAAACGTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                 <  2:93434/142‑1 (MQ=255)                                                                                                                                            atgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa               >  1:484216/1‑149 (MQ=255)                                                                                                                                             tgcttaacccAAAAAATAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa              <  2:195302/141‑1 (MQ=255)                                                                                                                                             tgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt            >  2:987967/1‑151 (MQ=255)                                                                                                                                              gcaggAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag             <  1:1045186/149‑1 (MQ=255)                                                                                                                                                atgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt            <  1:61647/148‑1 (MQ=255)                                                                                                                                                   aGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg       >  2:19723/1‑150 (MQ=255)                                                                                                                                                      cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  2:1073559/1‑150 (MQ=255)                                                                                                                                                      cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat    >  1:1142303/1‑150 (MQ=255)                                                                                                                                                       acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat  >  2:1141963/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      TTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAT  >  NC_000913/29337‑29635

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