Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,745,019	0	C	G	29.0%	85.1 / 17.7	69	G72G (GGC→GGG)	yiaN	2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base C (12/37);  new base G (12/8);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 1.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744783‑3745247                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     gCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑Tgcgggcgg                                                                                                                                                                                                                                    <  1:968087/239‑1 (MQ=255)         atgggtcTTAGTCATGGCTGGTCTGATTTGTCTGGGCTGTGTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCa                                                                                                                                                                                                                             <  2:2460866/232‑1 (MQ=255)           tgAGCTTAATCATGGCGTGGCTGATTTGTCTGGGCTTTCTGTTGGGTGGTGGTGCTATCGGTTTGCCTATTGCCTTTGCACTG‑TTGT‑TGTCCGTGGCGGCGGTAATGTTCTGGCTGGGCA‑TGTTTGATGTCCAGGTCATGGCCCAAACGCTGGTGAATGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑TTGCTGCCGGGTGAAATGATGAA‑TGCGGGCGGCTTGTCaaa                                                                                                                                                                                                                           <  2:73987/237‑1 (MQ=255)             aGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAACCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAATCATTAA‑ATAGGGGGGGGTGTCaaaaa                                                                                                                                                                                                                         >  1:2153877/1‑237 (MQ=255)              gCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                        <  1:817247/238‑1 (MQ=255)              cttgaagcATGGGTGTGTTGATTTTTCTGGGCTGTTGGTTGGGTGGTATTGCTTGCTGTTTGCCTGTTGCCTGGGCAGTG‑TTGT‑TGTGCGTGGCGGCGTTAATGTTCTGGCTGGACT‑TTTTTGATGTCCAGATCGTGGCCCAAACGCTGGTGAACGGAGCGGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                        <  2:952959/231‑1 (MQ=255)                  aaTCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGGGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGGAATCATGAA‑TGCGTGCGGCGTGTGGAAACGc                                                                                                                                                                                                                       >  2:1758910/1‑235 (MQ=255)                       tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGGAATCATGAA‑TTCGGGGGGCTTGTCAAAACGCAttgttcg                                                                                                                                                                                                               >  1:1017140/1‑236 (MQ=255)                       tGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAACTATGAA‑AGCGGGCGGCCTGTCAAAAAGACTTGTTGa                                                                                                                                                                                                               >  1:2669768/1‑238 (MQ=255)                       ggtcTGTTCAGATTTTTCTGGGCTGTCTTGTGGTTGGTTTTGTTATGGGTTGGCCTATTGCCGGGGCACTG‑TTGT‑TGGGCGGGGGGGCGTTAATTTTCTGGCTGGACA‑TGTTTGATTTCCAGATCATGGCACAAACGCTGGGGAAGGGAGCCGATAGCTTCTCCCGGCTGGCGTTTCCGTTC‑TTT‑GGGCGGGGGGGGGAAATCATGAA‑TGGGGGCGGCTTGTCAAAACGCAttgttgcc                                                                                                                                                                                                              <  2:2615660/236‑3 (MQ=255)                         gCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAAACATGAA‑TGGGGGCGGGTTGTCAAAACCCAttgttgt                                                                                                                                                                                                               >  1:1171265/1‑235 (MQ=255)                                tGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGGGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGAt                                                                                                                                                                                                      <  1:2015474/238‑1 (MQ=255)                                         tGGGCTTTCTGTTGGGTGGTGTTGCGTTCGGTTTGCCTAGTGCCGGGGCTCGG‑TTGT‑TGTGCGGGGGGGCGTTAATGTTCGGGGTGGACA‑TGTTTGATGTCCAGCTCAGGGCACAACCGGTGGTGAGCGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTT‑TTT‑GTGCTGGCGGGTGAAATAATGAA‑TGGGGGCGGCTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACtggtg                                                                                                                                                                                            <  2:1202037/239‑1 (MQ=255)                                         tGGGCTGTTTGTTGGGTGGTATTGCTATGGGTTTGCCTATGGCCGGGCAATTG‑TTGT‑TGGGGGGGGCGGCGTTAATGTTCTGGCTGGAAA‑TGTTTGATGTCCAGATCAGGCCACAAACGCTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTAAAATCATGAA‑TGGGGGCGGTTTGTAAAAACGCTTTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                            <  2:547568/239‑1 (MQ=255)                                         ggggCTGTTTGTGGGGTGCTACGGCTATCGGGTGGCCTTTGACCTGGGCAGTG‑TTGT‑GGGGGGGGGCGGCGTTATTGTTCTGGCTGGATA‑TGTTTGATGTCCAGTTCATGGCAAAAGCGCGGGTGAAGGGAGCAGATAGCTTCTCCTGGCGGGCGATTCGGTTT‑TTT‑GTGCGGGCGGGTGAAAACATAAA‑TGCGGGCGGCTTGAAAAAACGCATTGTTGACCTGCCGATGAAACtggtg                                                                                                                                                                                            <  2:533711/238‑1 (MQ=255)                                              tgtgtgtTGGGGGTTGTTGCTTCCGGTTTCCTTATGGCGTGGGGTCTT‑TTGT‑GGTGCGGGGGGGCTTAAATGTTCGGGCTGGACT‑TTTTTGATGTCCAAATGATGGCAAAAACGCGGGTGAACGGGGCCGATAGCTTCTCCCTGCGGGCGTTTCCTTTC‑TTT‑GTGCGGGCGGGTGAAATCATGAA‑GGGGGGGGGTTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCa                                                                                                                                                                                       <  2:2437664/239‑1 (MQ=255)                                                tctgtTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCaaa                                                                                                                                                                                                                           <  1:2157657/201‑1 (MQ=255)                                                     tGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAAcc                                                                                                                                                                                 <  1:2187816/238‑1 (MQ=255)                                                      gggTGGTATGGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAg                                                                                                                                                                               <  1:2681623/239‑1 (MQ=255)                                                       ggtggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACg                                                                                                                                                                                                                        <  1:1106394/197‑1 (MQ=255)                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATTAA‑TGCGGGGGGCTTGTCAAAAAGCATTGTTGACCTGCCGGTGAAACTGGGGGGGCGTAAACCAGGCgg                                                                                                                                                                           >  1:873633/1‑239 (MQ=255)                                                          ggtATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAATCATGAAGGGGGGGGGGCTGGTATAAACGCTTGTTTGACCTTCGGTTGAAATTGTGGGGTCAT‑AACCAGGCgg                                                                                                                                                                           >  2:272585/1‑239 (MQ=255)                                                            tgttgttttCGGTTTGCTTTTGGCCTGGCCCTTG‑TTGT‑TGGGCGGGGCGGCGTTAATGTTCGGGTGGGACT‑TTTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCCCCCTGCTGGCGTTTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGGGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                          <  2:2629885/230‑1 (MQ=255)                                                            tATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGGAAACCTGAA‑TGTGGGGGGCCTCTCTAAACGCATCGGTGACCTCCCCGTGTAAATGTTGGGGCATGCACCAGGCgtgg                                                                                                                                                                         >  2:1399180/1‑236 (MQ=255)                                                             tgtgtcttCGGTTTGCCTTTTGCCTGGGCGTTGTTTGT‑GGGGGGGGGCGGGGTTAATTTTCGGGTTGAACT‑TTTTGGATGCCCGGATAATGGCACAATCGTTGGTGAAGGGAGCGGATACCTTCTCCCTGCT‑GCGGGTTCGTTCGTTG‑GTGGGGGCGGTTGAAAATAAGAG‑TGCGGGGGAAGTGCCAAACCGCCTTGTTGGCCTGGCGATGAAACTGGTGGGGCATAAACCAGGCggg                                                                                                                                                                          <  2:1227064/231‑1 (MQ=255)                                                             aTTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑GGCGGGCGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                         <  1:2244965/238‑1 (MQ=255)                                                                 cTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTTGTGACCTGCCGATGAAACTGGGGGGGCATAAACCAGGCgggctggg                                                                                                                                                                     >  1:499411/1‑238 (MQ=255)                                                                   aTCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGGAATCATTAA‑TGCGGGCGGCTTGTCAAAACACATTGTTGTCCTGCCGGTGAAAATGGTGGGGCATAAAACAGGGgggcgggggt                                                                                                                                                                   >  1:1610233/1‑236 (MQ=255)                                                                    tCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGTA‑TTCGGGGGGCTTTTCAAAACGCATTTTTGATCTGCCGCTGAAACTGGGGGGGCGTAAAACAGGCggggtgggggac                                                                                                                                                                 >  1:2609436/1‑235 (MQ=255)                                                                      ggTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCCGCCGATGCAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGt                                                                                                                                                               >  1:2629885/1‑239 (MQ=255)                                                                      ggTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCCTTGTTGACCTGCCGGTGAAAAAGGGGGGGCGTAAACCAGGCGGGCCGGGCGACg                                                                                                                                                                >  1:2612488/1‑238 (MQ=255)                                                                          tctcTATTGCCGGGGCAGTT‑GTGT‑TGTGCGGGGCGGGGTAATTTTCCTGGCGGACCT‑TGTTTGATGTCCAGCTCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATCCCTTTC‑TTT‑GTGCGGGGGGGTGAAATCATGAA‑GGCGGGCGGTTTGTAAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGc                                                                                                                                                           <  2:2259315/236‑1 (MQ=255)                                                                           gCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAAACATGAA‑TGCGGGGGGCTTGTCAAAACGCACTGTTTACCTGCCGATGAAACTGGGGGGGGATAAACCAGGCGGGCGGGGGTACGgcggcg                                                                                                                                                          >  1:2259315/1‑239 (MQ=255)                                                                            ctttttGCCTGGGCGTTG‑TTGT‑TGGGCGGGGGGCGGTTAATGTTCTGGCTGGACA‑GTTTTGATGTCCAGATCGTGGCACAACCGGTGGGGAACGGAGCCGATAGCTTCTCCCTGCGGGCGTTTCCGTTC‑TTT‑GGGCGGGGGGGTGAAATAATGAA‑TGCGGGGGGTTTGTAAAAACGTATTGTTGAc                                                                                                                                                                                                              <  2:2202616/182‑1 (MQ=255)                                                                            ccTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCaaaa                                                                                                                                                                                                                          <  1:867307/174‑1 (MQ=255)                                                                                      gggCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAATCTTGAA‑TGCGGGCGGCTTGTCAAAAAGCATTGTTGGCCCGCCGCTgg                                                                                                                                                                                                    >  1:2202616/1‑185 (MQ=255)                                                                                        gCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTa                                                                                                                                              <  1:43861/238‑1 (MQ=255)                                                                                            tg‑ttgtgtgGGCGGGGCGGCGTT‑ATGGTTTGGGTGGGCATTGTTGGATGCCCAATTGATGGCAAAACCGCTGGGAAACGGAGCCGATATCTTCCCCCTGCTGGCGATTCCTTTT‑TTT‑GTGCTGGCGGGTGAAATCAGGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCCGGGCTACGTCGGCGTGCTGGCAGCTATGAt                                                                                                                                          <  2:1841596/239‑1 (MQ=255)                                                                                            tg‑ttgt‑tgtGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGAAATCATGAA‑TGCGGGGGGCTTGTCCAAAAGCATTGTTGACCTGCCGATGAAAATGGTGGGGGATaaa                                                                                                                                                                                   >  1:306651/1‑197 (MQ=255)                                                                                                            gcggcgTTAATGTTCTGGCTGAACA‑TGTTTGATTTCCAGATCATGGCACAAACGGTGGTGAACGGGGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTT‑TTT‑GTGCTGGGGGGTAAAATCATGAA‑TGCGGGCGGTTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCttt                                                                                                                            <  2:1909723/238‑1 (MQ=255)                                                                                                                                       gTTTGATGTCCAGATCATGGCACAAACGGTGGTGAACGGAGCGGATAGCTTCCCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑GGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                               <  1:2667753/193‑1 (MQ=255)                                                                                                                                               tCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCt                                                                                                                                                                                                                                  <  1:784124/102‑1 (MQ=255)                                                                                                                                               tCCAGATCATGGCACAAACGCGGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTT‑TTT‑GGGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCg                                                                                                                                   <  1:211976/197‑1 (MQ=255)                                                                                                                                               tCCAGATCATGGCAAAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTAAAATCATGAA‑TGCGGGCGGCTTGTCAAACCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCt                                                                                                                                               <  1:639636/185‑1 (MQ=255)                                                                                                                                                     tCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCa                                                                                                                                                                                                                      <  1:1422750/108‑1 (MQ=255)                                                                                                                                                         ggCACAACCGCGGGTGACGGGGGCCGATAGCTTCTCCCGGCTGGCGATTCCGTTC‑TTT‑GTGTTGGGGGGAGAAATCAGAAA‑GGGGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                             <  2:435816/129‑1 (MQ=255)                                                                                                                                                         ggCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAATCATGAA‑TGCGGGGGGCTTGTCaaaaag                                                                                                                                                                                                                        >  2:784124/1‑100 (MQ=255)                                                                                                                                                               aaCGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAATCATGAA‑AGCGGGCGGCTTGTCAAAAAGCCTTGTTGGCCt                                                                                                                                                                                                            >  2:1422750/1‑108 (MQ=255)                                                                                                                                                                  gtgggTGAACGGGGCCGATAGCTTCCCCCTGCTGGCGATTCCGTTC‑TTG‑GGGCTGGCGGGTGAAATCATGAA‑GGCGGGCGGCTTGTAAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGt                                                                                             <  1:1704108/213‑1 (MQ=255)                                                                                                                                                                   cTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGGGGCTTGTCAAAAAGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATaaa                                                                                                                                                                                   >  1:435816/1‑129 (MQ=255)                                                                                                                                                                      gTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACtggt                                                                                                                                                                                             <  1:1735501/116‑1 (MQ=255)                                                                                                                                                                          aCGGAGCCGATAGTTTCTCCCTGCTGGCGTTTCCGTTT‑TTT‑GTGCTGGGGGGTGAAATCATGAA‑GGCGGGCGGTTTGTCAAACCGTTTTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAg                                                               <  1:2187627/238‑1 (MQ=255)                                                                                                                                                                                ccGATGGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGGGGAATCATGAA‑TGGGGGCGGCGTGTCCAAAAGCCTTGTTTAACTGCCGGTGAAACCGGGGGGGGATaaa                                                                                                                                                                                   >  2:1735501/1‑116 (MQ=255)                                                                                                                                                                                 cgtaggtttCTCCCGGCGGGCAATTCCTTTT‑TTG‑GTGCTGGGGGGTGAAATCATGAA‑GGCGGGGGGCTTGTCAAAACGCATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACt                                                       <  1:64494/232‑1 (MQ=255)                                                                                                                                                                                   aTAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCGGGCGGGTAAAATCAGGAA‑TGCGGGGGGTTTGTAAAAACGCTTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCgggc                                                                                                                                                                         <  2:2590476/123‑1 (MQ=255)                                                                                                                                                                                        ttCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGGGGCTTGTCAAAACGCATTGTTGAcc                                                                                                                                                                                                             >  1:1605728/1‑82 (MQ=255)                                                                                                                                                                                        ttCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGAcc                                                                                                                                                                                                             <  2:1605728/82‑1 (MQ=255)                                                                                                                                                                                             ccTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCCTGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTTACCTTCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGt                                                                                                                                                               >  1:2590476/1‑123 (MQ=255)                                                                                                                                                                                                         ttcGTTT‑TTTGGTGCTGGGGGGTGAAATCAGGAA‑GGCGGGCGGCTTGTCAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCgggg                                <  1:1760634/237‑1 (MQ=255)                                                                                                                                                                                                              tC‑TTT‑GTGCTGGCGGGGGAAATCAGGAA‑GGCGGGCGGCTTGTCAAACCGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGAtt                          <  1:1762101/239‑1 (MQ=255)                                                                                                                                                                                                               c‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGAt                                                                                                                                                                                                      <  1:839506/66‑1 (MQ=255)                                                                                                                                                                                                                      acgtGGGGGGGGAAATCAGTAT‑GGCGGTGGTCTTGTCAAAACCATTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct                    <  1:1566589/236‑1 (MQ=255)                                                                                                                                                                                                                       gCTGGGGGGTAAAATCAGGAA‑GGGGGGCGGTTTGTAAAACCGCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCctct                    <  1:1615875/238‑1 (MQ=255)                                                                                                                                                                                                                           gCGGGTGAAATCATGAA‑TGCGGGGGGCTTGTCCAAACGCATTGTTGACCTGCCGGTGAAACtgggg                                                                                                                                                                                            >  2:839506/1‑64 (MQ=255)                                                                                                                                                                                                                               gTGAAAGCATGAA‑TGCGGGCGGCTTAAAAAAACGCTTTTTGGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTg                                                                                                                                                     <  2:1529952/101‑1 (MQ=255)                                                                                                                                                                                                                                 gAAATCATGAA‑GGGGGGCGGCTTGTAAAACCGAATTTTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTAt         <  1:1399180/239‑1 (MQ=255)                                                                                                                                                                                                                                  aaatattGAG‑TGGGGGCGGTTTGTAAAAACGAATTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGGGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCg                                                                                        <  1:186375/153‑1 (MQ=255)                                                                                                                                                                                                                                        tGAA‑TGCGGGCGGCTTGTCAAACCGCTTTTTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCa  <  1:142714/239‑1 (MQ=255)                                                                                                                                                                                                                                         gAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGa                                                                                                                                           >  1:1529952/1‑101 (MQ=255)                                                                                                                                                                                                                                          aa‑TGGGGGCGGCTGTCCAAACCCCTTTGTTGCCCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCgcgc    <  1:694466/235‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     GCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG‑TTGT‑TGTGCGGGGCGGCGTTAATGTTCTGGCTGGACA‑TGTTTGATGTCCAGATCATGGCACAAACGCTGGTGAACGGAGCCGATAGCTTCTCCCTGCTGGCGATTCCGTTC‑TTT‑GTGCTGGCGGGTGAAATCATGAA‑TGCGGGCGGCTTGTCAAAACGCATTGTTGACCTGCCGATGAAACTGGTGGGGCATAAACCAGGCGGGCTGGGCTACGTCGGCGTGCTGGCAGCTATGATTATGGCGAGCCTTTCTGGTTCTGCCGTTGCGGATACCGCCGCCGTTGCCGCTCTACTGGTGCCGATGATGCGCAGCGCCAACTATCCGGTCAACCGGGCGGCGGGGCTGATTGCCTCTGGCGGCATTATCGCGCCA  >  NC_000913/3744783‑3745247

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