Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,297	0	G	C	39.8%	152.3 / 42.4	94	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 (16/40);  new base C (37/0);  total (54/40)
Fisher's exact test for biased strand distribution p-value = 1.73e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  NC_000913/3423150‑3423445                                                                                                                                                      |                                                                                                                                                        aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGc                                                                                                                                                      <  2:1333715/151‑1 (MQ=255) aaaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCcgg                                                                                                                                                       >  1:725151/1‑150 (MQ=255)  aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCCCCCC‑CCGGGcc                                                                                                                                                     >  1:266711/1‑151 (MQ=255)   aTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCACCCC‑CCCGGcc                                                                                                                                                     >  1:849782/1‑150 (MQ=255)      agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGCATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCCGCCCGCCACCCT‑CCggg                                                                                                                                                       >  2:709585/1‑145 (MQ=255)      agagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑Gtt                                                                                                                                                 <  2:1020349/150‑1 (MQ=255)       gagCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GCCCCGGCCGCCCCCCC‑CCCGGcccgtg                                                                                                                                                 >  2:904631/1‑146 (MQ=255)          cGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCCGCCCGCCACCCC‑CCcgg                                                                                                                                                       >  2:742935/1‑141 (MQ=255)           gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGTCGCCCCCCAACCC‑CCTGGCC‑Gt                                                                                                                                                  >  2:1299687/1‑144 (MQ=37)              cGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑aa                                                                                                                                       <  2:78504/151‑1 (MQ=255)                aGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACg                                                                                                                                     <  1:104597/151‑1 (MQ=255)                 gggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGt                                                                                                                                    <  2:250902/151‑1 (MQ=255)                    cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑CACAtt                                                                                                                                   >  1:1035337/1‑149 (MQ=25)                    cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTTCCCCCCCCCCC‑CCCGGCC‑GCTGGtt                                                                                                                                             >  2:797947/1‑140 (MQ=37)                     aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑aa                                                                                                                               <  1:1047441/151‑1 (MQ=255)                     aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTATGGCAGTTTATGGCGGGCGGTCCT‑CCCTCCACCCC‑CCCGTCC‑GGGTCTTCGcccagctta                                                                                                                                  >  2:1288614/1‑142 (MQ=255)                         ttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCTCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTC‑AAATcc                                                                                                                           >  1:1346088/1‑151 (MQ=25)                          tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GTCCTGCCCGCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTC‑AAATccc                                                                                                                          >  1:416553/1‑150 (MQ=255)                              gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑CTTTCTTCGC‑AACATTC‑AAAACCGCTcc                                                                                                                      >  1:1135709/1‑151 (MQ=255)                                aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTcccc                                                                                                                    >  1:1388850/1‑150 (MQ=255)                                 ggTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAAcccccccc                                                                                                                     >  1:766768/1‑149 (MQ=21)                                       ccATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGAt                                                                                                              <  1:1466829/150‑1 (MQ=255)                                        cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCGCTCCccgcggggtt                                                                                                            >  1:1272993/1‑147 (MQ=18)                                        cATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GGTGCTTCGC‑AACGTTT‑AAATCCCCTCCccgcggggtt                                                                                                            >  1:1311088/1‑147 (MQ=255)                                         aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGCCCGCCCCCCC‑CCGCGCC‑GTGGGTTCCC‑AACGCTT‑CCATCCGCCCCCCGCCGGttt                                                                                                            >  1:397191/1‑150 (MQ=255)                                         aTTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTg                                                                                                           <  1:323719/151‑1 (MQ=255)                                           tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGTCCGCCACCCC‑CCCGGCC‑GTTTCTTTGC‑TACGTTCAAAACCCGCCCCCGCCGGttt                                                                                                             >  2:539827/1‑148 (MQ=18)                                              ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCACCCC‑CCCGGCC‑GTTGCTTCGC‑AACGGTC‑AAAACCGCCCCCCGCGGAttttttc                                                                                                        >  1:1455465/1‑145 (MQ=18)                                               tGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTGCCGCCCCCCACCCC‑CCCGGCC‑GGTGCTGCGC‑AACGGTC‑ACAACCCCTCCCCCCGCATTTGTTCttc                                                                                                     >  2:148815/1‑149 (MQ=18)                                                gAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCCTCCCCCCACCCC‑CCCGGCC‑GTTTCTTCCC‑AAAGCTC‑AAAATCGCCCCCCGCGGGTTTGTTCTAc                                                                                                     >  2:1246775/1‑150 (MQ=255)                                                    gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCgg                                                                                                 <  1:175471/150‑1 (MQ=255)                                                     ttAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCCGTCCCCCGCCCT‑CCCGGGC‑GTTGCTTTGA‑AACCTTC‑AAAACACCTCCCCGGGGGATTGTCCTACTc                                                                                                   >  2:863213/1‑147 (MQ=17)                                                      tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GGTAGTTCGC‑AACCATC‑CAATTCCCACCCCGCCGCTTTGTCCCACTCcgg                                                                                                >  2:1512243/1‑149 (MQ=17)                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtgt                                                                                         <  1:477220/151‑1 (MQ=255)                                                           cGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAgtg                                                                                          <  2:1324881/150‑1 (MQ=255)                                                              aaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGCTGCCCGCCCCCCC‑CCCGGCC‑GTTGCTTTGC‑AACCTTT‑AAATCCCCTCCCGGCGGAGTTGTTCTACTCcgggg                                                                                              >  1:732653/1‑141 (MQ=17)                                                                   aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGG‑GGCCTCCCCCCCCCCCC‑CCCGGCC‑GGTTTGTTGC‑AACCTTC‑ACATTCGCCCCCCGGGGGTTTGTCCCAccctgg                                                                                                >  2:845058/1‑131 (MQ=255)                                                                   aaaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGA‑GGCCTGCCCGCCACCCC‑CCCGGCC‑GGTGCTTCCC‑CAAGTTT‑ACAACTGCCCCCCGCGGGTTTTTTCTACTCTCGAGGGGGTTCTCCGac                                                                                 >  2:1194788/1‑151 (MQ=255)                                                                    aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGGCCTCCCCCCCCCCT‑CCCGGCC‑GGGTCTTTGC‑AACCATC‑AAATACGCCCCCCGCGGGGTTGTTCTTCTCGGGAGGGGGttttcc                                                                                    >  2:1106600/1‑143 (MQ=255)                                                                    aaaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGCC‑GGCCTGCCCCCCCCCCTCCCGGGCC‑GCT‑CTTCGC‑AAAATTT‑AAATACCCCCCCCGCGGGGTTGTCCTACTCCGGGGAgggtgctccc                                                                                   >  2:1259320/1‑142 (MQ=255)                                                                      aatgCCATCCTTTCGGCTGTCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTCCCCGCCACCCT‑CCGGGTC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                              <  2:237985/148‑1 (MQ=17)                                                                      aaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaa                                                                              <  1:1173829/151‑1 (MQ=255)                                                                        ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                            <  2:884933/151‑1 (MQ=255)                                                                        ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                            <  2:1424295/151‑1 (MQ=255)                                                                        ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGacaaaca                                                                            <  2:109735/151‑1 (MQ=255)                                                                          ccATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCT‑CCCGGCC‑GTTGCTTCCC‑AACGTTC‑AAATCCCCTCCCCGCGGGTTTGTCCTACTCCGGGGAGTGGTCACCGccaaacaa                                                                           >  1:1573775/1‑150 (MQ=17)                                                                           cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                         <  1:323475/151‑1 (MQ=255)                                                                           cATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCC‑CCGGGCC‑GTTGCTTCGC‑CAAGATC‑AAATCCCCCCCCCGCGGGTTTGTCCTACTCCGGGGAGTGTTTACCCacaaacaa                                                                           >  1:954444/1‑149 (MQ=17)                                                                           cATCCTTTCGGATGGCCTTTCGCTCGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAacaaca                                                                         <  1:1283453/151‑1 (MQ=255)                                                                              ccTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGGCCACCAACCC‑CCCAGCC‑GTTGCCTACC‑AAACGTC‑CAATCCCCTCCCCGCGGGTTTTTTCCACTCCGGGCAGAGGTCACCCacaaagaaaaaa                                                                       >  1:963655/1‑146 (MQ=9)                                                                               cTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGGTTTGC‑CACCTTC‑AAATACCCCCCCCGCGGGTTTTTTCTTCTCCGGGGAGGGTTTACCCacaaaaaa                                                                           >  1:468817/1‑145 (MQ=255)                                                                                tttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCCCCCCCCGCGGGTTTTTTCTACTCGGGGGAGTGTTCACCCAAAAAAAACAGATa                                                                     >  1:6686/1‑150 (MQ=14)                                                                                 ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCACCCC‑CCCGGCC‑GTTTCTTCCC‑AACGGTC‑AAATCCCCTCCCCGCGGGTTTTTCCCCCTCGGGGGAGAGGTCACCCACAAACAAACGATaa                                                                    >  1:917828/1‑150 (MQ=14)                                                                                  tCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCT‑CCCGGCCGGTTGTTTC‑C‑AAAGTTT‑AAAACCCCCCCCCGCGGATTTGTCCTACTCCGGAGAGGGTTTACCCACCAACAAAAGAaaaa                                                                   >  1:1456553/1‑150 (MQ=9)                                                                                   cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCCCTCCCCGCGGGTTTTTCCTACTCGGGGGAGTGGTTACCGACAAACAACAGATaaaa                                                                  >  1:28093/1‑150 (MQ=14)                                                                                    ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACa                                                                <  1:1146701/151‑1 (MQ=255)                                                                                    ggATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCCGCGGGTTTTTCCTACTCGGGAGAGGGTTCACCCACAAACAACAGATaaaa                                                                  >  1:603749/1‑149 (MQ=14)                                                                                      aTGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                              <  2:1153615/151‑1 (MQ=255)                                                                                      aTGGCCTTTCCCTTGATTTGATGTCTGGCAGTTTATGGCGGTC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACaaa                                                              <  2:814750/151‑1 (MQ=255)                                                                                           cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTT‑AAAACCCCTCCCCGCGGGTTTTTCCTACTCCGGGGAGGGTTTACCCACAAACAAAAGATAaaaaaaaa                                                             >  1:1148344/1‑147 (MQ=14)                                                                                            tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCCCCCCGCGGGTTTTTCCTACTCCGGAGAGTGTTCACCGACAAACAACAGATAAAAAAAAAAGccc                                                        >  1:1282376/1‑151 (MQ=14)                                                                                            tttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCC‑CCCCGCC‑GTTGCTTCCC‑AACGTTT‑CAAACCCCCCCCCGCGGGTTTTTTCTCCTCCGGGGGGGGTTTGCCCCCCAACAACAGATAaaaaaaaaagac                                                         >  2:1454432/1‑146 (MQ=9)                                                                                             ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCCCCCT‑CCGGGCC‑GGTGCTTCTC‑AACCTTC‑AAAAACGCCCCCCGCGGATCTGTCCTACTCGGGGGGGGGGTCACCGACAAACAACCGATAaaaaaaaa                                                             >  2:1575074/1‑145 (MQ=14)                                                                                             ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACGGATAAAACAAAAGGcc                                                         <  1:802482/149‑1 (MQ=255)                                                                                              tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCa                                                       <  1:1030182/150‑1 (MQ=255)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAg                                                      <  2:101309/150‑1 (MQ=255)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGt                                                     <  1:1300633/151‑1 (MQ=40)                                                                                               cGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GGCCTCCCCCCCACCCC‑CCCGGGC‑GGTGCTTCGC‑AAAGGTC‑AAATTCGCCCCCCGCGGAGTTTTCCCACCCGGGAGAGTGTTGACCCGCAAACAACAGATAAAAAAAAAAGcc                                                         >  1:118277/1‑147 (MQ=255)                                                                                                  ttgatttgatGTCTGGCAGTTTATGGCGGGC‑GGCCTCCCCGCCCCCCT‑CCCGGGC‑GTTGCTTCGC‑CACGTTC‑TAAACCGCCCCCCGCGGGTTTTTCCTCCTCCGGGGGGGGGTTTCCCTCAAACAACAGAGAAAACaaaaaacc                                                         >  2:747326/1‑140 (MQ=2)                                                                                                   tgatttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCt                                                   <  2:1270203/149‑1 (MQ=40)                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCCCCCC‑CCCGGCC‑GTTGCTTCGC‑AACGTTC‑AAAACCCCTCCCCGCGGGTTTTTCCTTCTCCGGGGAGAGTTTACCGACAAACAAAAGATAAAAAAAAAAGGCCCGTCTTcc                                                >  1:458888/1‑150 (MQ=2)                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                               <  1:516253/151‑1 (MQ=38)                                                                                                     atttgatGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCC‑CCGGGCC‑GTTTCTTCGCAAACGTTC‑‑AATCCCCCCCCGGCGGGATTTTCCTCCTCAGGAGAGTGTTTACCCACAAACAAACAATAAAAAAAAAAGACCAGTCTTcc                                                >  1:1618754/1‑150 (MQ=2)                                                                                                        tgatGTCTGGCAGTTTATGGCGGGC‑GTCCTTCCCCCCACCCC‑CCCGGGC‑GGGGCTTCGC‑AACGTTC‑AAAACCGCCCCCGCCGGGGTTTTCCTATTCGGGGGAGGGTTCACCCACCAACAAAAGCTAAAAAAAAAAGccccagcttccccc                                             >  2:1517153/1‑139 (MQ=2)                                                                                                               tGGCAGTTTATGGCGGGC‑GTGCTGCCCCCCCCCCC‑CCCGGCC‑GGTTCTTTTC‑CACGTTC‑AAATTCGCCCCCCGCCGGATTGTTCTTCTTGGGGGAGGGGTTACCGACCAACAACAACAAAAAAAAAAAGCACCGCCTTCCGCCTGCGCCtt                                     >  2:1090978/1‑151 (MQ=2)                                                                                                                  cAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTg                                  <  2:36922/151‑1 (MQ=37)                                                                                                                       tATGGCGGGC‑GTCCTGCCCGCCCCCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGTATTTGTCCTACTCGGGAGAGTGTGCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTTATCCTTTTGtttt                              <  2:557555/150‑1 (MQ=255)                                                                                                                       tATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtttt                              <  1:1388401/150‑1 (MQ=38)                                                                                                                           gcgggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                         <  2:1400080/151‑1 (MQ=38)                                                                                                                            tgggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTg                         >  2:950476/2‑150 (MQ=38)                                                                                                                             gggc‑gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                       <  1:803324/151‑1 (MQ=38)                                                                                                                                  gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                    <  1:1346141/150‑1 (MQ=40)                                                                                                                                  gTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCt                   >  2:323276/1‑151 (MQ=40)                                                                                                                                   tCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTg                  >  2:846037/1‑151 (MQ=40)                                                                                                                                    ccTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                 <  2:624464/151‑1 (MQ=40)                                                                                                                                       gcacGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGACCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTGGATGTCggg                 <  2:334136/145‑1 (MQ=255)                                                                                                                                       gCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                <  1:1483683/149‑1 (MQ=40)                                                                                                                                       gCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa               <  2:724860/150‑1 (MQ=40)                                                                                                                                           gCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc          <  1:1155050/151‑1 (MQ=40)                                                                                                                                           gCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc          >  1:1117816/1‑151 (MQ=40)                                                                                                                                            ccACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGCGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGCCTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc         >  2:1375753/1‑151 (MQ=25)                                                                                                                                              aCCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt        >  1:1591921/1‑150 (MQ=40)                                                                                                                                               cccT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCt        <  2:1392542/149‑1 (MQ=40)                                                                                                                                               cccT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTa       <  2:159768/150‑1 (MQ=40)                                                                                                                                                  t‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct   >  2:279597/1‑151 (MQ=40)                                                                                                                                                  t‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGCCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTActct   >  2:817153/1‑151 (MQ=255)                                                                                                                                                    ccGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACtctc  >  2:525588/1‑151 (MQ=40)                                                                                                                                                      |                                                                                                                                                        AAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGC‑GTCCTGCCCGCCACCCT‑CCGGGCC‑GTTGCTTCGC‑AACGTTC‑AAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTC  >  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 < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑