Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,191	0	A	C	25.9%	40.8 / 35.6	59	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (25/18);  new base C (8/7);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTA  >  NC_000913/4295893‑4296475                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAAcg                                                                                                                                                                                                                                                                                                >  1:428659/1‑301 (MQ=25)        gTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGACGCTTAACGCGTCTTa                                                                                                                                                                                                                                                                                         >  1:714946/1‑301 (MQ=25)                ttgggccagtgactgCGTTTAGGCCGCATCCTGCATGAACGCCTGATGAGACGCTTAACGCGTCTTATAAGGCCGACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTAGGCCG‑GAGCAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc                                                                                                                                                                                                                                                                                  <  2:5423/286‑1 (MQ=11)                   tgtcggatgcggcataaacgccttaTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCCACGCCTGATGCTAAGCTTCACGCGTCTTATCAGGCCTACg                                                                                                                                                                                                                                                                              >  2:18118/25‑301 (MQ=11)                             ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATGCCAACCAGCGCAATACTCTGCTTCAGCGTGATGTTGTAGGTCTGAGAAGGCGGTTACTCTGCATCCGACATCCGTGCCTGATGCTACGCTTAACGAGTCTTATCAGCCCTATGCCATTCAgc                                                                                                                                                                                                                                                                     >  2:16035/1‑300 (MQ=11)                             ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGCTTTAGCGTGCTTTTGCAGCCCGGCTAAGACGTTTATGCCGCATTCGACTTCAACGCCTGATGCGCAGCCTCACACGTCGTATGAGGCCTACGCCAGaagtcg                                                                                                                                                                                                                                                                    >  2:255561/1‑296 (MQ=11)                             ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGTCTTCTTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTGTGGCCTCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCACACAgcg                                                                                                                                                                                                                                                                    >  2:255559/1‑301 (MQ=255)                             ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGatcgga                                                                                                                                                                                                                                                                    >  1:710143/1‑296 (MQ=25)                               cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGACGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCCTGATTTAGCGTGACTTTGTTGGCGGGATAAGGCGTTTATGCCGCATCCGATATCAACGCCTGATGCGGCGCTTAAAGCTTCTTATCAGGCCTACGCCAGACAGCGCa                                                                                                                                                                                                                                                                  >  2:247018/1‑301 (MQ=255)                                 aagacgtgtgctcttccgatgtcGCCTAATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGAATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                <  1:685353/279‑1 (MQ=255)                                     cgccttaTCCGCCATCACCGCCTGTTGCGACGCTTAACGCGTCTTATAAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                                                                                                                                                                                            <  1:18119/295‑1 (MQ=255)                                          cATCAGACAGCACCGACTGATGCTACGCGTATGGCGATGTAGAAGGCCGTAGACAGAAATCGAAATAGCCGGATTTAACGGAAGTTGGTGAGTCGGGTAAGGAGTTGATGTCG‑CATCCGACATCAACGGATGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAtt                                                                                                                                                                                                                                                       <  2:130230/301‑1 (MQ=11)                                            tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                                                                     >  1:339872/1‑301 (MQ=34)                                            tCCGACATCAACGACTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                                                                     <  1:247020/301‑1 (MQ=21)                                                 aatcaaggccgaggcagaccctgagtgggTCGTATCAGGCCAAGGACAGAAGGAGCAAGAGAGGTATTGAGAGTAATTTTGTAGGTCGGATAAGGCGTTTAAGCCG‑CATCCGACATCAACGCCGGGTGCGACGCTTGCCGCGGCTTATCAAGACTACGCCAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTg                                                                                                                                                                                                                                                <  2:714941/273‑1 (MQ=11)                                                  aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGa                                                                                                                                                                                                                                               >  1:699397/1‑301 (MQ=34)                                                    cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCCATCGCCTGATTTAGCGTGAtt                                                                                                                                                                                                                                             >  2:203048/1‑301 (MQ=21)                                                    cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtt                                                                                                                                                                                                                                             >  1:406893/1‑301 (MQ=34)                                                    cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGTCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATCTTGTACGCAGGCTAAGGCGTTTTTGCCGCATCCAACATCAACGACTGTTGCGATGCTTAACGCGTCTTTTCAGGCATACGCCAGACAGCCCAATAGCACGCTTTAGACTGAtt                                                                                                                                                                                                                                             >  2:396728/1‑301 (MQ=11)                                                     aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAAAAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAttt                                                                                                                                                                                                                                            <  1:396730/301‑1 (MQ=21)                                                       cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCTGCGCAATAGCCTGTTTTAGCGTGATTTTGTAGGTCGGATAATGCATTTATGCCGCTTCCGATATCAACGCCTGATGCGATGCTTAACGCGTCTTATCAGGCATACGACAGTCAGCGCAATatatcggcacagcttcgtgtcg                                                                                                                                                                                                                                          >  2:685349/1‑280 (MQ=11)                                                        cacgagctctcccgatgtcgcgTCTTATCAGCCCCACGCCATACAGCGCAAGAGCCTGATTTAGAGTGATTGTGTAGGTCGGATAAGGCGTTTACGCCG‑CATTCGACATCAACGCCTGGTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                                                                                         <  2:579473/283‑1 (MQ=11)                                                         ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                                                                                         <  1:760615/300‑1 (MQ=35)                                                             aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                                                                                                                                                                                     >  1:213003/1‑300 (MQ=35)                                                                      tgctcttccgatctCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                                                                                                                                                                                           <  1:297056/288‑1 (MQ=17)                                                                          cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGGaagagcacacg                                                                                                                                                                                                                       >  1:579476/1‑291 (MQ=17)                                                                                    cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCAGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGcagctcggaagaggg                                                                                                                                                                                                             >  2:297054/1‑287 (MQ=12)                                                                                            cGCCAGACAGCGCAATAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGTCAGACAGCGCAATAGCCTGCTTTAGCGTGATCTTGTCAGTCGGATAATGCGTTTCTGCCGCATCCTACCTCATTGCCCGATGCTACGCCTCGCGCTTCTCATCCGTCCTACTTCTGTCCTCACAATAGCCTTCTCTATCGTTTTTTTGTAGGTCCGAtccggtatctttgctccgtcggccctc                                                                                                                                                                                                      >  2:172170/1‑274 (MQ=2)                                                                                                    aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTCATTTTGTAGGTCGGATAAGGCGATTCTGCCGCATCCGAACTCAACGCCTGATGCGACGCTTAACGCGGCGTATCAGGCCTACCACAGCCCGCGGAATAGCCTGATTTAGCGTGATTTTGTAGATCGGAAAAGGCGTTTACGCCGCATCCGAAATCAATGCATGa                                                                                                                                                                                             >  2:368446/1‑301 (MQ=11)                                                                                                       tgtgctcttccgatctgCGTGATTTTGTAGGTCGGATAAGGCGTATACGCCG‑CATCCGACATCAGCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                                                                          <  1:224116/285‑1 (MQ=255)                                                                                                              ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                                                                                   <  1:203049/301‑1 (MQ=21)                                                                                                              ccTGATTTAGCGAGATTTTGTAGGTCGGATACGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                                                                    <  1:306904/300‑1 (MQ=21)                                                                                                                 gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCttgc                                                                                                                                                                                 >  1:560068/1‑298 (MQ=18)                                                                                                                     tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCTATGCCTGAGTCGACGCTTGCCGCGt                                                                                                                                                                            >  2:265130/1‑301 (MQ=12)                                                                                                                       gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCCTCTTATCAGGCCTACGCAAGACAGCGCAACAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTCGCCGCATCCGGCATCAATGCCTGATGcagatcggaagagcgt                                                                                                                                                                           >  2:224115/1‑285 (MQ=11)                                                                                                                       gCGTGATATTGTGAATCAGATAATGGATATACGCCG‑CACATGAAAGCAAAGCCTGATGCGACCATTGCCACCTCTTCGCAGTCCTACGCGAGACAACGCAATAAGCTGATTTCGCGTGATTTTGTAGTTCGGATAAGGCGTTTATGCCGCATAGGGCATCACAGCCAGATACGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                                                                           <  2:92119/300‑1 (MQ=2)                                                                                                                        ctcgtgccacaagcgacgctagtcagagcaCGCCG‑CATCCGACAACAACGACTGATGCGGAGGTTGCCGCGTCTTATCAGGCCTACGCCGGACAGCGCAATAGCCTGATTTAACGTGATTTTGTAGGGCGGATAAGGCGATTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                                                                         <  2:668168/272‑1 (MQ=11)                                                                                                                          tGATTTTGTAGGTGGGATAAGGGGTTTACTGAG‑AAGCGTACATGAACGCCTGATGCGAGGCTTGCCGCGACTGATCGGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                                                                        <  2:192388/300‑1 (MQ=11)                                                                                                                                 gTAGGTCGGATAAGG‑GGTAACGGCGAAAGCCGTAATGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCTTAAGCCAGAGAGCGAATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCAACCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                 <  2:698360/300‑1 (MQ=11)                                                                                                                                    ggTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                              >  1:386878/1‑299 (MQ=14)                                                                                                                                    ggTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                                                              >  1:389659/1‑299 (MQ=9)                                                                                                                                        ggATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttaa                                                                                                                                                         >  1:239357/1‑295 (MQ=2)                                                                                                                                                      cgccg‑cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTagatcggaagagcacacgtctgaactccag                                                                                                                                           >  1:668172/1‑272 (MQ=11)                                                                                                                                                                                                                                                                ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTACGCCAGTAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                 <  2:239356/300‑1 (MQ=255)                                                                                                                                                                                                                                                                  cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                               <  1:265132/300‑1 (MQ=9)                                                                                                                                                                                                                                                                  cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                              >  2:86370/1‑301 (MQ=12)                                                                                                                                                                                                                                                                    tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                            >  1:194253/1‑301 (MQ=12)                                                                                                                                                                                                                                                                        cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGTCGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCACTAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATAAATGCCTGATGCGACGTTTGCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                        >  2:283377/2‑301 (MQ=11)                                                                                                                                                                                                                                                                        cgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATCTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                        >  2:162267/2‑301 (MQ=255)                                                                                                                                                                                                                                                                             cgcgcgggagcAACGCCTGATGCGAGGTGTGCCGCGGATTATCAGGCATACGCAAGACAGCGAAAAAGCCAGAGTTAGCGTGATTTTGTAGGTCGGATGAGGCGTTTATGCCGCACCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                    <  2:213002/290‑1 (MQ=11)                                                                                                                                                                                                                                                                             cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                   <  1:16036/301‑1 (MQ=255)                                                                                                                                                                                                                                                                               tccgcaaTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaat                 <  1:255563/295‑2 (MQ=255)                                                                                                                                                                                                                                                                               cccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                 <  1:255561/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                ccGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGCGATTTTGTAGGTTTGATAAGGCGTTTA‑‑CGCTGCATCTGACATCCATGCCTGATGCGACGCTTGTCGCCTCCTATCAGGCCTATCATAACCGTTGGTTAATTTTTCaaa                  >  2:51161/1‑299 (MQ=11)                                                                                                                                                                                                                                                                                    aaTCAACGCCTTATTGAATGCTTGACGCGTCTTCTCAGGCCTACGCCAGAAAGAGGCATACCCTGATTTAGGGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCGTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAAGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc            <  2:699392/300‑1 (MQ=14)                                                                                                                                                                                                                                                                                           gCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTGTAGGTCGGACAAGGCGTTTACGC‑C‑GCATCTGACATCAATGCCTGATGCGACGCTTGCCGCGTCGTCTCACGCCTGTTCTAACCGTTGGTTAATTTTTCAAAATCTTCACATCCg     >  2:126106/1‑301 (MQ=11)                                                                                                                                                                                                                                                                                            ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg    <  1:86371/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                            ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg    >  1:164634/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                              tGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa  >  1:174523/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                 ATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCG‑CATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTA  >  NC_000913/4295893‑4296475

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑