Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,297	0	G	C	49.0%	122.6 / 60.3	99	intergenic (+103/+126)	yhdZ/rrfF	putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base G (9/41);  new base C (48/0);  total (58/41)
Fisher's exact test for biased strand distribution p-value = 3.71e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  NC_000913/3423150‑3423445                                                                                                                                                       |                                                                                                                                                         aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTGC‑TT‑CCCCCCCCCCTCCgcg                                                                                                                                                        >  1:402892/1‑148 (MQ=255)  aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGcc                                                                                                                                                      <  1:494083/151‑1 (MQ=255)       gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTTC‑TG‑CCCGCCACCCCCCCGGCCGtt                                                                                                                                                   >  2:369347/1‑149 (MQ=255)       gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGc                                                                                                                                                 <  2:1047331/151‑1 (MQ=255)         gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCCCCCGGGCGTTGCtt                                                                                                                                               >  1:312751/1‑151 (MQ=255)           gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCG‑GGCC‑TG‑CCCCCCCCCCCCCGCGcccc                                                                                                                                                    >  1:238122/1‑142 (MQ=255)            tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTG‑CCCCCCCCCCCCCGGCCGGTTCCt                                                                                                                                                >  2:941429/1‑148 (MQ=37)            tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCGCCCCCCTCCCGGGCGTTGCtt                                                                                                                                               >  2:669291/1‑148 (MQ=255)             aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TGCCCCCCCCCCCCCCGGGCCGt                                                                                                                                                    >  2:1111338/1‑143 (MQ=255)              cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCa                                                                                                                                           <  1:806998/150‑1 (MQ=255)              cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CT‑CCCCCCCCCCCCCCGGCCGCTGCtt                                                                                                                                               >  1:374328/1‑146 (MQ=37)                aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCCCCCTCCCGGCCGGTGCttttcaa                                                                                                                                          >  1:312522/1‑144 (MQ=25)                   gCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCC‑CG‑CCCCCCCCCCCCCCGGCCGCTGCTTTGCAA‑CG‑tt                                                                                                                                    >  2:838693/1‑150 (MQ=255)                    cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCACCCCCCCGGCCGGTTCTTCGCAA‑CG‑tt                                                                                                                                    >  2:446566/1‑149 (MQ=255)                    cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCCCCCCCCCCCCGGGCGGTTCTTCTCAA‑CA‑tt                                                                                                                                    >  1:193806/1‑149 (MQ=255)                     aTTTTTATCGCAGGTAATCCATTTATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCaa                                                                                                                                 <  1:218189/151‑1 (MQ=255)                     aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCCCCCGGCCGTTGCTTTGCAA‑CG‑tt                                                                                                                                    >  2:470479/1‑148 (MQ=255)                       ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCC‑TC‑CCCCCCCCCCCCCCGCCCGTT‑CTTCCCAA‑AG‑CTCaaa                                                                                                                                >  1:842570/1‑150 (MQ=255)                        tttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCAGGCC‑CC‑CCCCCCCCCCCCCGGGCCTTTCTTTGGCAA‑AG‑TTCa                                                                                                                                  >  1:544822/1‑148 (MQ=255)                          tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATcc                                                                                                                             <  2:135579/150‑1 (MQ=255)                          tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCCCCCGGCCGTTTCTTCGCAACCG‑TCCAAATc                                                                                                                              >  1:1105866/1‑150 (MQ=25)                           aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCC‑TG‑CCCCCCACCCCCCCGGCCGGTGCTTTGCAA‑AC‑TTCAAATc                                                                                                                              >  2:974906/1‑148 (MQ=255)                           aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCg                                                                                                                            <  2:797733/150‑1 (MQ=255)                           aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCCCCCCCCCGGCCGTTGCTTTTCAA‑CC‑TTCAAAAcccc                                                                                                                           >  2:753417/1‑149 (MQ=21)                           aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCCCCACCCCCCCCGCCGGTGCGTTGCAA‑CA‑ttata                                                                                                                                 >  1:510321/1‑142 (MQ=255)                            tCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TT‑CCCCCCACCCCCCCGGCCGCTTCCTCTCCA‑CC‑TTCAAAtac                                                                                                                             >  1:861853/1‑146 (MQ=255)                               cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTccc                                                                                                                       <  2:32085/151‑1 (MQ=255)                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCCCCCCCCCGGCCGTTGCTTTTCAA‑CG‑TTCCAATCCGCccccc                                                                                                                      >  1:94954/1‑150 (MQ=21)                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCCCCCGGCCGGTGCTTCGCAA‑CG‑TTCAAAACCGCTccc                                                                                                                       >  2:68504/1‑150 (MQ=255)                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCCCCCGGCCCCTGCTTCGCAA‑CC‑TTCAAAAccccccccc                                                                                                                      >  2:723289/1‑150 (MQ=255)                                 ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCcgg                                                                                                                     <  1:431736/151‑1 (MQ=255)                                    aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCGCCACCCCCCCGGCCGTTGCTTTGCAA‑AG‑TTTAAATCCCCCCCcggc‑gg                                                                                                                 >  2:825053/1‑151 (MQ=255)                                    aaTCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑CTCC‑CG‑CCCCCCACCCCCCCGGCCGCTTCTTCGCAA‑AG‑CTCCAATCCCCTCCccgc‑gg                                                                                                                 >  2:317661/1‑151 (MQ=255)                                       ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGAtt                                                                                                              <  2:933652/151‑1 (MQ=255)                                             aTTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCt                                                                                                        <  2:196289/151‑1 (MQ=255)                                              ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCt                                                                                                        <  2:663636/150‑1 (MQ=255)                                                gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TC‑CCCCCCCCCCCCCCCGGCGCTGCTTCTCAA‑CG‑TTCAAAACCCCTCCCCGC‑GGGCTTTTCCttc                                                                                                      >  1:945508/1‑148 (MQ=255)                                                 aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTc                                                                                                    <  2:718235/151‑1 (MQ=255)                                                 aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCCCCCTCCCCGCCGGTTGTTCTCAA‑AC‑TTTAAATACCCCCCCCGC‑GGGTTTGt                                                                                                           >  2:652242/1‑144 (MQ=255)                                                 aaTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCACCCTCCCGGCCGCTGCTTCGCAA‑AG‑TTCAAATTCGCCCCCCGC‑GGAttttt                                                                                                           >  2:325136/1‑142 (MQ=255)                                                  aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCGCCCCCCCCCCGGCCGTTTCTTCGCAA‑CG‑TTTAAAACCCCTCCCCGC‑GGGTTTTTCCTACTc                                                                                                    >  2:217791/1‑150 (MQ=17)                                                   tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CT‑CCCGCCCCCCCCCCGGCCGGTGCTTTGCAA‑CC‑TTTAAATACCCCCCCCCC‑GGGTTTTTCCTAc                                                                                                      >  1:272861/1‑147 (MQ=255)                                                   tGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCCCCCCCCCCCCGGCCGGTGCTTCGCAA‑CA‑TTCCAAACCACCCCCCGC‑CGGttttttcttc                                                                                                      >  2:1098105/1‑140 (MQ=255)                                                    gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCgg                                                                                                  <  2:646897/150‑1 (MQ=255)                                                     ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCCCCCCCCCGGCCGGTTCTTTGCAA‑CG‑TTTAAATCCGCTCCCCGC‑GGGGTTTTCCTACt                                                                                                     >  2:703122/1‑146 (MQ=17)                                                     ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCCACCTCCCGGCCGTTGCTTTGCAA‑CG‑ATCAAAACCCCCCCCCGC‑GGATTTGTCCTACt                                                                                                     >  2:469170/1‑146 (MQ=17)                                                     ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCACCCCCCCGGCCCTTGCTTCGCAA‑AC‑TTTAAATCCCCCCCCCGC‑GGGTTTGTCCTACTCtgg                                                                                                 >  2:723173/1‑150 (MQ=17)                                                       aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGgaga                                                                                              <  2:510324/151‑1 (MQ=255)                                                       aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TC‑CCCCCCACCCTCCCGGCCGTTGATTTGCAA‑CC‑TTCAAATCCCCCCCCGGC‑GGGGTTTTCCTTCTCcgg                                                                                                 >  2:737114/1‑148 (MQ=17)                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTTC‑TT‑TCCCCCCCCCTCCCGGCGCTTTTTTCGCAA‑CG‑TTCAAAATCCCTCCCCGC‑GGGTTTTTCCttccctgg                                                                                                 >  2:597584/1‑137 (MQ=17)                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAgtgt                                                                                          <  2:1006466/151‑1 (MQ=255)                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAgtg                                                                                           <  2:613725/150‑1 (MQ=255)                                                             aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTCCTCGGGAGAGTGTTc                                                                                        <  2:488387/151‑1 (MQ=255)                                                              aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCa                                                                                       <  2:505922/151‑1 (MQ=255)                                                                aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCAcc                                                                                     <  1:1038710/151‑1 (MQ=255)                                                                 gCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCCCCCCCCCCGCCGCTTCTTTCCCA‑CG‑TTTAAATTCCCCCCCCGC‑GGGTTTTTTCCACCCCCGAGAGTGTTCAcc                                                                                     >  2:648388/1‑150 (MQ=12)                                                                  cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGa                                                                                   <  1:1107368/151‑1 (MQ=255)                                                                      aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGa‑caaa                                                                              <  1:372305/151‑1 (MQ=255)                                                                       aGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCCCCCCCCCGGGCGTTGCTTCGCAA‑AG‑TTCAAATCCGCTCCCGGC‑GGGTTTTTCCTCCTCCGGGGAGGGGTCACCCa‑caaa                                                                              >  1:14706/1‑150 (MQ=14)                                                                         gCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCCCCCCCCCCCCGGCCGGTGCTTCCCAA‑CG‑TTCCAATCCCCCCCCCGC‑GGGTTTTTCCTACTCCGGAGAGTGTTCACCCa‑c                                                                                 >  2:10700/1‑145 (MQ=255)                                                                          ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑TT‑CCCGCCCCCCCCCCGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCCGC‑GGGTTTTTCCTACTCGGGAGAGGGTTCACCGa‑caaacaa                                                                           >  2:487199/1‑150 (MQ=17)                                                                           cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCCGGCCGTTGCTTCCCAA‑CG‑TTCAAATCCGCTCCCCGC‑GGGTTTTTCCTACTCGGGGGAGGGTTTACCGa‑caaacaa                                                                           >  2:490280/1‑149 (MQ=17)                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGa                                                                       <  1:626811/150‑1 (MQ=255)                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGAt                                                                      <  2:389342/151‑1 (MQ=255)                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGAt                                                                      <  2:438660/151‑1 (MQ=255)                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CG‑CCCCCCCCCCCCCCGGCCGTTGCTTTGCCA‑CG‑TTTCAATTCGCCCCCCGC‑CGGTTTTTCCTTCTCCGGGGAGGGGTTACCCC‑CAAaaaacaaa                                                                       >  1:389083/1‑148 (MQ=9)                                                                              acTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGCTTATGGCGGTC‑CGCC‑CT‑CCCCCTTCCCACCAGGCCGTTTTTTCGCAA‑AA‑ATGAAATCACCACCCGGG‑GGGTTTTTCCCTCTCTGGAGAGTGATCACCGc‑caaacaaacg                                                                        >  2:735465/2‑146 (MQ=255)                                                                                tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CT‑CCCGCCCCCCCCCCGGCCGCTGCTTTGCAA‑CC‑TTTCAATCCCCCCCCCGC‑GGGGTTTTCCCTCTCGGGAGAGTGTTCACCCACcaaacaccgaa                                                                       >  1:208124/1‑144 (MQ=255)                                                                                 ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATaa                                                                    <  2:925824/150‑1 (MQ=255)                                                                                 ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCACCCCCCCGGCCGGTGCTTTGCAA‑CG‑TTTCAAACCGCCCCCGGC‑GGGTTTTTCCTACTCCGGAGAGTGTTCACCCC‑CAAACAACAGATaa                                                                    >  2:330169/1‑150 (MQ=14)                                                                                   cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCC‑CT‑CCCCCCCCCCCCCCGGGCGTTGGTTTGCAA‑AA‑TTTAAAACCCCCCCCCGCGGGGTTTGT‑TTTCCCCGGAGAGTGGTTACCCA‑CAAAAAAAAAATaaaa                                                                  >  2:1013703/1‑150 (MQ=255)                                                                                    ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GCCC‑CG‑CCCGCCCCCCCCCCGGCCGCTGCTTTACAA‑CGTTTCAAAACC‑CTCCCCGC‑GGGTTTTTCCCCCTCGGGGAAGGGGTCACCGG‑CAAAAAACaaaaaaaaa                                                                 >  2:308296/1‑149 (MQ=255)                                                                                       tGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TT‑CCCCCCACCCCCCCGGCCGTTGCTTCGCAA‑CC‑TTCAAAACCGCTCCCCGC‑GGGTTTTTCCTACTCCGGGGAGTGTTTACCCA‑CAAACAACAGATaaaaaaaa                                                              >  2:1200558/1‑150 (MQ=14)                                                                                        ggCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAg                                                            <  2:1048846/151‑1 (MQ=255)                                                                                           cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCCCCCCCCCGGCCGTTGCTTCGCAA‑AG‑TTCAAATCCGCTCCCCGC‑GGGTTTTTCCTTCTCCGGAGAGGGTTCACCGA‑CCAACAACAGATAAAAAAAAAAGcc                                                         >  2:866297/1‑151 (MQ=14)                                                                                            tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑CT‑CCCCCCACCCTCCCGGCCCTTGCTTCGCAA‑CG‑ATCAAATCCGCTCCCCGC‑GGATTTGTCCTACTCCGGAGAATGGTCCCCCA‑CAAACAACAGAAAAAAAAAAAAGccc                                                        >  2:687879/1‑151 (MQ=14)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCa                                                       <  2:1104151/149‑1 (MQ=255)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGt                                                     <  1:421312/151‑1 (MQ=40)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCCCCACCCCCCCGGCCGGTTCTTTGCAA‑CC‑TTCCAATCCGCCCCCCGC‑GGGGTTGTTCTACCCCGGAGAGGGGTCCCCCA‑CAAACAACACATAaaaaaaaaaggcc                                                        >  1:1058904/1‑143 (MQ=2)                                                                                                  ttgatttgatGTCTGGCAGTTTATGGCGGGC‑GGCC‑TG‑CCCCCCACCCTCCGGGCCGTTGCTTTTCCA‑AG‑TTCAAATACGCTCCCCGC‑GGGTTTTTCCTCCTCGGGAGAGGGTTTCCCCA‑CAAAAAACAGAAAAAAAAAAAAGccc                                                        >  1:577196/1‑145 (MQ=9)                                                                                                   tgatttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAATAGATAAAACAAAAGGCCCAGTCtt                                                  <  2:83796/150‑1 (MQ=255)                                                                                                   tgatttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                  <  2:1000585/150‑1 (MQ=40)                                                                                                   tgatttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCCCCCTCCCGGCCGGTGCTTTGCAA‑CG‑TTCAAATCCGCCCCCCGC‑GGGCTTGTCCTACTCCGGGGAGGGTTTACCCC‑CAAACAAAAGATAAAACAAAAAGACCAGGCtt                                                  >  1:715853/1‑150 (MQ=2)                                                                                                    gatttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTcc                                                <  1:51089/151‑1 (MQ=39)                                                                                                    gatttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                 <  2:662458/150‑1 (MQ=39)                                                                                                      tttgatGTCTGGCAGTTTATGGCGGGC‑GTAC‑TG‑CCCGCCCCCCCCCCGGCCGCTTCCTTGCAA‑CG‑TTCTAATCCCCGCCCCGC‑GGGTTTTTCCTACGCGGGGGAGGGTTTACCCA‑CCAAAAACAAATAAAACAAAAAGACCAGGCtttc                                                >  2:61977/1‑147 (MQ=2)                                                                                                       ttgatGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGa                                              <  1:156423/150‑1 (MQ=38)                                                                                                               tGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATa                                                                     <  2:350915/119‑1 (MQ=35)                                                                                                                ggCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                    <  1:634681/151‑1 (MQ=35)                                                                                                                ggCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                    <  1:686947/151‑1 (MQ=35)                                                                                                                 gCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                   <  1:723279/151‑1 (MQ=37)                                                                                                                 gCAGTTTATGGCGGGC‑GTCC‑CG‑CCCCCCCCCCCCCCGGCCGGTGCCTCTCAA‑CG‑GTCAAAACCCCCCCCCGC‑GGGTTTTTCCTACTCGGGAGAGGGTTTCCCGC‑CAAACAACAGAAAAAAAAAAAAGCCCAGTCTTCCGACTGAGCCtttt                                   >  1:316178/1‑151 (MQ=2)                                                                                                                        aTGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTa                             <  1:1131210/150‑1 (MQ=38)                                                                                                                                   tCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                  >  2:1032508/1‑151 (MQ=40)                                                                                                                                        g‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCTACTGAGCCTTTTGTTTTATTTGATGTCTGGCAg              <  2:971429/151‑1 (MQ=255)                                                                                                                                           ccGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt            <  1:853211/151‑1 (MQ=40)                                                                                                                                           ccGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt            <  2:1096819/151‑1 (MQ=40)                                                                                                                                             gCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc          >  1:1023820/1‑151 (MQ=40)                                                                                                                                                   cTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActc    >  2:920240/1‑151 (MQ=40)                                                                                                                                                     ccGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc  >  1:1066279/1‑151 (MQ=40)                                                                                                                                                       |                                                                                                                                                         AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCC‑TG‑CCCGCCACCCTCCGGGCCGTTGCTTCGCAA‑CG‑TTCAAATCCGCTCCCGGC‑GGATTTGTCCTACTCGGGAGAGTGTTCACCGA‑CAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  NC_000913/3423150‑3423445

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