Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,423,297	0	G	C	47.0%	98.9 / 40.0	66	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 (4/31);  new base C (31/0);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 8.34e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423151‑3423446                                                                                                                                                   |                                                                                                                                                          aaTCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGcc                                                                                                                                                       <  1:709101/151‑1 (MQ=255)    cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCGCCCCCCCCCcgg                                                                                                                                                         >  2:517812/1‑146 (MQ=255)    cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGcc                                                                                                                                                       >  1:895403/1‑148 (MQ=255)    cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGtt                                                                                                                                                    <  2:653704/151‑1 (MQ=255)    cAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGt                                                                                                                                                     <  1:554128/150‑1 (MQ=255)        gCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCCCCCCCCgg                                                                                                                                                          >  2:317628/1‑141 (MQ=255)          gTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCCGCCCGCCCCCCCCCcgg                                                                                                                                                         >  2:1001562/1‑140 (MQ=255)           tACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCGCTGCCCCCCACCCCCCCGGGCGTTGC‑tt                                                                                                                                               >  2:138959/1‑148 (MQ=255)            aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCa                                                                                                                                           <  2:422810/151‑1 (MQ=255)            aCGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCGCCCCCCCCCCCGCCGCTg                                                                                                                                                   >  2:704226/1‑144 (MQ=255)                 ggCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑tt                                                                                                                                    <  2:187051/151‑1 (MQ=255)                   cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCACCCCCCCGGCCGTTGC‑TTCGCCA‑AG‑tt                                                                                                                                    >  2:887520/1‑149 (MQ=255)                   cATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGACCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCa                                                                                                                                  <  2:309982/151‑1 (MQ=255)                    aTTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGC‑TTCGCAA‑CG‑tt                                                                                                                                    >  1:855323/1‑148 (MQ=255)                     tttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGTTGC‑TTCGCAA‑AG‑TTTaa                                                                                                                                 >  1:189493/1‑150 (MQ=25)                      ttttATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCTCCCGGCCGTTTC‑TTTTCAA‑CA‑tttcaa                                                                                                                                >  2:369011/1‑146 (MQ=25)                         tATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATcc                                                                                                                             <  2:35649/150‑1 (MQ=255)                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGc                                                                                                                           <  1:156024/151‑1 (MQ=255)                          aTCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCCCCCGGGCGTTGC‑TTCGCAA‑CG‑TTCaaaaaccc                                                                                                                           >  1:470214/1‑146 (MQ=255)                             gCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCGCCCCCCTCCCGGCCGCTGC‑TTCGCCA‑CG‑GTTAAATCCGcccc                                                                                                                        >  2:515833/1‑151 (MQ=21)                              cAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCCCCCCCCCGGCCGGTGT‑TTCGCAA‑AG‑TTCAAAAcccccccc                                                                                                                       >  1:626875/1‑151 (MQ=255)                               aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCcg                                                                                                                      <  1:643351/151‑1 (MQ=255)                               aGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCcg                                                                                                                      <  2:111847/151‑1 (MQ=255)                                 gTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCCCCCCCCCGGCCGCTGC‑TTCGCAA‑AG‑TTCAAAACCCCTccc                                                                                                                       >  1:142490/1‑148 (MQ=21)                                     tCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCGCCACCCCCCCGGCCCTTGC‑TTCGCAA‑CG‑TTCAAATCCCCCCCcggcgg                                                                                                                  >  1:370160/1‑149 (MQ=21)                                          tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCGCCCCCCCCCCGGCCCGTGC‑TTCGCAA‑CC‑GTCAAAACCGCCCCCCGCGGGGTTGt                                                                                                            >  2:432435/1‑150 (MQ=255)                                          tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCCGGCCGTTGC‑TTCGCAA‑CG‑GTCAAAACCGCCCCCCGCGGGttttt                                                                                                            >  1:50038/1‑148 (MQ=18)                                          tAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCTTCCCCCCACCCCCCCGGCCGTTGC‑TTCGCAA‑AG‑TTCAAAACCCCCCCCCGCGGGTTTGt                                                                                                            >  1:716214/1‑150 (MQ=255)                                             ttGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCGCCCCCCCCCCCCCCGGCCGGTGCTTTTCCAA‑CG‑‑TCAAATCCCCTCCccgcggggtttt                                                                                                            >  2:212751/1‑141 (MQ=18)                                                 aTGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCTCCCGGCCGTTGC‑TTCGCCA‑CG‑TTCAAAATCGCTCCCCGCGGGTTTTTCCttc                                                                                                       >  1:779721/1‑146 (MQ=18)                                                   gTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑Cgg                                                                                                  <  2:727752/150‑1 (MQ=255)                                                     tAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGGCCCTCCCCCCCCCCTCCCGGGCGGTTC‑TTTGCAA‑AC‑GTCAAAACCCCCCCCGGCGGGTTTTTCCTtct‑ctg                                                                                                  >  2:200062/1‑146 (MQ=12)                                                      aGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGgaga                                                                                              <  2:604937/151‑1 (MQ=255)                                                       gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGagag                                                                                             <  2:230071/151‑1 (MQ=255)                                                       gTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCCGCCCCCCACCCCCCCGGCCGTTGC‑TTTGCAA‑CG‑TTCAAAACCGCCCCCCGCGGGTTTTTTCTACT‑CCGgag                                                                                               >  1:305163/1‑149 (MQ=17)                                                            aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCGGCCGGTGC‑TTCTCAACCT‑TTCAAACCCCCCCCCGCCGGGTTT‑TCCTTCT‑CCGGggagag                                                                                           >  2:467333/1‑146 (MQ=255)                                                            aaaaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGC‑TTCGCAA‑CG‑TTCAAATCCCCTCCCGGCGGGGTTTTCCTtct‑c                                                                                                    >  2:19745/1‑139 (MQ=17)                                                              aaGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGTTGC‑TTCGCAA‑CG‑TTCAAAACCCCTCCCCGCGGGTTTTTTCTACT‑CCGGAGAGGGtt                                                                                         >  1:846504/1‑148 (MQ=17)                                                               aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGGGTCCTGCCCCCCCCCCCCCCGGCCGCTGG‑TTCGCAA‑CC‑CTCAAAACCGCTCCCCGCGGGGTTTTTCTACC‑CGGGGGAGGGTTCAcc                                                                                     >  2:584361/1‑151 (MQ=255)                                                               aGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCAcc                                                                                     <  1:844509/151‑1 (MQ=255)                                                                 cAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTGCTCCCCCCCACCCTCCCGGCCCCTGC‑TTTGCTA‑CG‑TTCAAATCTCCTCCCCGCGGGTTTTTCCCACT‑CGGGggaggggttacc                                                                                     >  1:1004562/1‑141 (MQ=17)                                                                    aaaGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGacaa                                                                                <  1:1007384/151‑1 (MQ=255)                                                                       ggCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGacaaaca                                                                             <  1:782447/151‑1 (MQ=255)                                                                           aTCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAg                                                                         <  1:614458/151‑1 (MQ=255)                                                                                  cGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCCCCACCCCCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCCCCCCCCGCGGGTTTTTCCTACT‑CGGGAGAGGGTTTACCGACAAACAACAGATaaaa                                                                   >  1:324593/1‑150 (MQ=14)                                                                                         ccttctgCTTAATTTGATGCCTGGCAGTTTATGGCGGGCGGCCTCCCCCCCCCCCCCCCCGCCCCTGC‑TTCCCCA‑CC‑TTCAAATACCCCCCCCCCGGGCTTTTCCCTCC‑CCCGGGAGGGGTCCCCCACCAACAACAGACAaaaaaaaa                                                              >  2:546068/7‑148 (MQ=0)                                                                                         ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTTCCCCCCACCCCCCCGGCCGTTGC‑TTCGCAA‑CC‑TTCCAAACCGCTCCCCGCGGGTTTTTCCTACTCCGGGGGAGGGTT‑ACCCACCAACAACAAATAaaaaaaaaa                                                             >  1:762322/1‑148 (MQ=9)                                                                                         ccTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCCGGCCGTTTC‑TTCGCAA‑CG‑TTCAAATCCCCTCCCCGCGGGTTTTTCCTACT‑CGGGAGAGTGTTTACCGACAAACAACAGATAAAACaaaaa                                                             >  1:29303/1‑148 (MQ=14)                                                                                             tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCCCCCCCCGGCCGTGGC‑TTTGCCA‑CG‑TTCCAAAACGCCCCCCGCGGGTTTTTTCTACC‑CGGGGGAGTGTTTACCCACCAAAAACAAAAAAAAAAAAAAGccc                                                         >  1:790872/1‑149 (MQ=9)                                                                                             tCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCa                                                        <  2:91706/150‑1 (MQ=255)                                                                                                 ttgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCCCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCCCCCCGCGGGTTTGTCCTTCT‑CGCGAGAGTGTTTACCCACAAACAACAGATAAAAAAAAAAGCCCACTCtt                                                   >  1:803422/1‑151 (MQ=14)                                                                                                  tgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                   <  1:998698/150‑1 (MQ=40)                                                                                                   gatttgatGTCTGGCAGTTTATGGCGGGCGGCCTGCCCCCCACCCTCCCGGCCGTTGC‑TTCGCAA‑AG‑TTCAAATCCGCTCCCCGCGGGTTTGTCCTACT‑CCGGGGAGTGGTCACCCACAAACAACAAATAAAAAAAAAAGCCCAGTCTTcc                                                 >  1:628019/1‑151 (MQ=255)                                                                                                         atGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa                                           <  1:682361/151‑1 (MQ=38)                                                                                                           gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTg                                            <  2:567473/148‑1 (MQ=38)                                                                                                               ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                     <  1:251668/151‑1 (MQ=35)                                                                                                                gCAGTTTATGGCGGGCGTCCTGCCCCCCCCCCCCCCGGCCGGTGC‑TTTCCAA‑AGCTTCAAA‑ACGCCCCCCGCGGGTTTTTCCTACT‑CCGGGGGGGGTTTACCGCCAAAAAAAAAAAAAAAAAAAAAGCCCAGGCTTCCCACTGAGCCcttt                                    >  2:639026/1‑151 (MQ=2)                                                                                                                       aTGGCGGGCGTCCTGCCCCCCACCCTCCCGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTTTCCTACT‑CGGGAGAGGGTTCACCCACAAACAACAGATAAAACAAAAAGCCCAGTCTTCCGCCTGAGCCTTTTGTTTTAt                             >  2:190001/1‑151 (MQ=11)                                                                                                                         ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAttt                           <  2:971136/151‑1 (MQ=38)                                                                                                                         ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                            <  1:350461/150‑1 (MQ=38)                                                                                                                            gggcgTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGAt                        <  1:550828/151‑1 (MQ=38)                                                                                                                              gcgTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                      <  1:517827/151‑1 (MQ=39)                                                                                                                              gcgTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGt                      <  1:503883/151‑1 (MQ=39)                                                                                                                                     gCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAg               <  2:626859/151‑1 (MQ=40)                                                                                                                                      cccGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGt              >  1:516878/1‑151 (MQ=40)                                                                                                                                       ccGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGtt             <  1:766255/151‑1 (MQ=40)                                                                                                                                         gCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc           <  1:137618/151‑1 (MQ=40)                                                                                                                                          ccACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc          >  2:817718/1‑151 (MQ=40)                                                                                                                                                  cGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCg  <  2:554650/151‑1 (MQ=40)                                                                                                                                                   |                                                                                                                                                          AATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGC‑TTCGCAA‑CG‑TTCAAATCCGCTCCCGGCGGATTTGTCCTACT‑CGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCG  >  NC_000913/3423151‑3423446

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