Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,191	0	A	C	29.1%	58.2 / 27.5	56	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (15/24);  new base C (7/9);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295899‑4296482                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc                                                                                                                                                                                                                                                                                                   <  2:711404/299‑1 (MQ=25)    tatttggttgtacggcttcgctgTTTATGCCGTCTCGGTCTTCATCTCTTGATGCGACGTTTATCGCCTCTTATCAGGCTGACGTCAATCAGCTCACTCGACTGCGTTGGCG‑GATTGTTGTAGGTCGGATAACGCGTTTACTCAGCAT‑CCGACATCACCGCCTCCTGCGACGTTTCCCGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATACTGCGTTTACGCCGCATCCGACATCATCGCCTGATGCTACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                              <  1:567197/279‑1 (MQ=2)                  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACATCGCACTCGCCTGATTTATCTTTATTTTGGTGCTCGGAACAGTCCTTAACGCCGCACCCGCCCTCCAGCCCTGATCTGCCCCTTGACTCGTCTTATCCGGCCGACCCCAg                                                                                                                                                                                                                                                                                 >  1:45971/1‑300 (MQ=11)                   aTAAGGCGTTTAGGCCGCATCCGACATCAACTCCTGATGCGAAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTATGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                               <  2:148255/301‑1 (MQ=17)                    acggaggggtatgacgaatgagcaaTCAACACATGAGGCGACGGTTACCGTGGCTAATCAGGCCTAGGGCAGACGGCGCAATAGCATGAGTTAGCGTGATT‑TTTTATGTCGGAGAAGGCGTTTACGCCGCAG‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                               <  1:163702/276‑1 (MQ=11)                         cttttAGGCCGCACCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                                          <  2:212070/298‑1 (MQ=25)                         cGTTTATGCCGCAGCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                                          <  2:574131/300‑1 (MQ=25)                           tttATGCCTCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                       <  2:329577/301‑1 (MQ=25)                           tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                                                                                                                                                                                       >  2:286086/1‑301 (MQ=35)                                      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGTCGTTTACGCCGCAT‑CCGCCATCAAAGCCTGTTGCGACGCTTGCCGCGTCTTCTCATGCCTACTCCACACGTCCCAATAGTCTGATTTCGCGCGATCTTGTGCGTCGGTTCATGACTTTCGTACGCCTCCGACTCCCACGCCTGATGCCACGCGTACCGCGCCTTATCTGGcttctgcccgacctcgacctcgcactattca                                                                                                                                                                                                                                                            >  1:170037/1‑271 (MQ=2)                                      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGACTGATTTAGCGTGATGTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGACTGATGCAACGCTTAGCGCGTCTTATCAGGCCTACGCCAGCACGCGCCTTACCCTGATTTa                                                                                                                                                                                                                                                            >  1:723730/1‑301 (MQ=11)                                      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                                                                                                                                                                                            >  2:773162/1‑301 (MQ=34)                                         caCAGCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGCCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                                                                                                                                                                                         <  2:371726/300‑1 (MQ=21)                                           cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCCATAGCCTGACTTATCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCACGGCCTGATTCGACCATTACCGCCTCTTCACTGTCCTCCCCTAGACAGCGTACTATCCtctcttaccca                                                                                                                                                                                                                                                        >  1:83790/1‑290 (MQ=11)                                                  gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                                                                                                >  2:600797/1‑301 (MQ=34)                                                   gctgtggctggggttttctcgTTTTGTCAGGCCTACGTGAGACGGCGAATTAGTCTAATTTCGGGTGAGT‑TTGTGGGGCGGGTAAGGCGTTTACTGCGCAT‑CCGAAATCCACGCTGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                                                                                                                                                                                                <  1:500064/281‑1 (MQ=11)                                                    ctggtgacactcttcacgtgggtgctCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATG‑TTGTAGGTCGGATGAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCTGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                              <  2:9269/276‑1 (MQ=11)                                                    cggATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                              <  2:616610/299‑1 (MQ=21)                                                    cTGATGCGAAGCTTAACGCGTCTTGTCAGGCCTACGCCAGAAAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCGGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                                                                                                                                                                                              <  2:271354/301‑1 (MQ=21)                                                                        tCTTATCAGGCCTACGCCAGACAGCTCAATAACCTGATTTAGCGTGATT‑TTGTAGGTCGGAGAAGGCGGTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGc                                                                                                                                                                                                                          <  2:768538/301‑1 (MQ=17)                                                                           tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTccgccg                                                                                                                                                                                                                        >  2:383733/1‑300 (MQ=17)                                                                                     aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACCCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTACGCCGCATCCGACAtt                                                                                                                                                                                                             >  2:561920/1‑300 (MQ=25)                                                                                     aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCGACGCTTCACGCGTCTTATCAGGCCTACGCAAACCCGCGCAATAGCCTGATTTCGCGTCATTTTGTAGGTCGGATAAGTCGTTTACACCGCAACCgcacta                                                                                                                                                                                                             >  1:563175/1‑296 (MQ=11)                                                                                     aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAt                                                                                                                                                                                                              >  2:760801/1‑300 (MQ=25)                                                                                     aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAt                                                                                                                                                                                                              >  2:178609/1‑300 (MQ=25)                                                                                          aGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                                               <  2:780537/230‑1 (MQ=33)                                                                                           tcctgctcatgactcttaTTTAGCGTGATT‑TTGGAGGTCGGATAAGGCTTTTACGCCGCAT‑CCTACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTCTCAGGCCTACGCCATACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                       <  2:624752/284‑1 (MQ=11)                                                                                            agaGCGCAATAGACTGATTTAGCTTGAGG‑TTGTAGGTCGGGTAAGGCGTTTACGCCGCAG‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                       <  2:69241/298‑1 (MQ=255)                                                                                             cAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                              <  2:37237/212‑1 (MQ=34)                                                                                             cAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                                                                                                                                                                                                              >  1:37235/1‑212 (MQ=34)                                                                                             cAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACCCTTAACGCTCATTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTAACCGTGAATATGCAGGTAGAATAAgggtctgaagcctcacccaccaatagtgctaa                                                                                                                                                                                                     >  1:119437/1‑271 (MQ=11)                                                                                             cAGCGCAATAGCCTGATTTAGCGTGATT‑GTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                                                                       <  2:119442/299‑1 (MQ=21)                                                                                                cgcAATAGCCTGATTTGTCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                                                                        <  2:728900/295‑1 (MQ=21)                                                                                                cgcAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                                                                    <  2:754159/299‑1 (MQ=21)                                                                                                cgcAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCATGATGCGACGCTTACCGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGACTTACCGTGATTTTGTCGGTCGGATATGTCGTTTACGccccttaccacccgaatccctgat                                                                                                                                                                                                   >  1:866441/1‑278 (MQ=11)                                                                                                       tggttgatggtgtgttatgtggcaGGTCGGGTCGGGTGGTTACGCCGAAT‑GCGACATCAACGCCTGAGGCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGGCAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAGGGCGTTTATGCCGCATCCGAGATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                                                                            <  1:200072/278‑1 (MQ=11)                                                                                                       gCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCATGCCTACGCCAGACCTCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGCCGTTTACGCCGCATCCGACATCCATGCCtgctgccgtgat                                                                                                                                                                                           >  2:322471/1‑295 (MQ=11)                                                                                                       cacgaaTTTAGCGTGAGT‑TTGTAGGACGATTAACGCATTTAAGGGGCATGCGGATAT‑ATGGCCTATTATGCCGCCTGCCCGGTCTTATCAGGGCTACGCCAAAAGGCCCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTGGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                                                                            <  2:640329/295‑1 (MQ=11)                                                                                                            ttgtactgTGATT‑TTGGAGTGCCGCTACGCCGTCTACGACGATG‑CCTACATCAACGCCGGAGGCGAGTGTTGCATAGGTGTATAAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGGGACCCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCttgccg                                                                                                                                                                                      <  1:561901/294‑5 (MQ=11)                                                                                                                      tttgtggtagttcggacaggggttcttcgccgcAT‑CATACATAGAGGCCGGAGCCGACGCGTTCCGCGGCTTATCAGGCCTAAGCCAGAACGCGCAATAGCCTGATTTATCGTGATTTGGTAGGGCGGAAAAGGCGTTTAGTCCGCACCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                             <  1:476408/274‑1 (MQ=11)                                                                                                                      aTT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCCGTGTCTTAACa                                                                                                                                                                            >  1:132032/1‑301 (MQ=12)                                                                                                                      aTT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATAAATGCCTGATGCGACGCTTTGCGCGttctgtca                                                                                                                                                                            >  1:730850/1‑294 (MQ=255)                                                                                                                                 gcGGATAAGTCGTTTAGGCCGCAT‑CCGACATCACCGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                                                                                   <  2:855080/299‑4 (MQ=11)                                                                                                                                  cGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTTTCCGGCCTACGCAAGCCAGCGCACTAGCCTGTATTCGCGTGCTTTTGTTGGTCTGACAAGGCGTTTACGGCGCCTCCGACCTCCATACCTGTTCCGACTCttgcgatgtcatattctgactgtcttt                                                                                                                                                                 >  1:722071/1‑276 (MQ=11)                                                                                                                                                                                                                                                        taatgtctttagcgcccATCGGACATCAACGCCTGCGGCGACTGTTACCGCTTCTTCGCAGACCTACGCCGGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  1:773132/285‑1 (MQ=11)                                                                                                                                                                                                                                                        aaaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  2:723757/300‑1 (MQ=2)                                                                                                                                                                                                                                                        aaaGGCGATGACGACGAATCCGACATCAAGGACGGATGCGACGCGGGCGGCGTCTTATCAGGCCTACGCCAGGCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAGAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  2:83794/300‑1 (MQ=11)                                                                                                                                                                                                                                                         aaGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                           <  2:141058/299‑1 (MQ=2)                                                                                                                                                                                                                                                              gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                    >  2:456544/1‑301 (MQ=12)                                                                                                                                                                                                                                                                        cATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                           >  2:484974/1‑300 (MQ=255)                                                                                                                                                                                                                                                                          tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTATTATCCGGCCTATCTTTACCGTTGGTTAATTTTTCaaaa                        >  1:366256/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                    aCGCCTGATGCGACGCGTGCCGCGTCTTATCAGGGCTACGCACGACAGCGCAATAGTCTGATTTAGCGTTATTGTGTAGGTCGGATAAGGCGTTTAGGCCGCAGCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc              <  2:733324/301‑1 (MQ=14)                                                                                                                                                                                                                                                                                    aCGCCAGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc              <  2:536797/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                     cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTCGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACCTCCATTCCTGATGCTACCCTTGCCCCGGCTTATCAAGCCTACCTTAACCGTTGGGTACGTTTTTAACCtttttcaacc              >  1:156568/1‑294 (MQ=11)                                                                                                                                                                                                                                                                                               caCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc   <  2:797731/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                aCGCTTGCCGCTTCTTAGCAGGCCGACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  <  2:704273/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                        CGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATT‑TTGTAGGTCGGATAAGGCGTTTACGCCGCAT‑CCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4295899‑4296482

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