Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,378,968	0	C	T	26.3%	64.3 / 11.4	38	E98K (GAA→AAA)	menD	2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylate synthase; SEPHCHC synthase
Reads supporting (aligned to +/- strand):  ref base C (3/25);  new base T (10/0);  total (13/25)
Fisher's exact test for biased strand distribution p-value = 6.05e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGT  >  NC_000913/2378730‑2379198                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          tCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTtc                                                                                                                                                                                                                                           <  2:526106/239‑1 (MQ=255)             ggTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAAcc                                                                                                                                                                                                                                <  1:816980/238‑1 (MQ=255)               ttCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCg                                                                                                                                                                                                                              <  1:68240/238‑1 (MQ=255)                                tGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTg                                                                                                                                                                                                                 <  1:378725/234‑1 (MQ=255)                                        aCTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATaga                                                                                                                                                                                                        <  1:153389/235‑1 (MQ=255)                                         cTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTTAGAGAATCAGG‑T‑TTTTTTCCGTTTAA‑CCGGCTTCCACCAATTCCgggcaaaaattt                                                                                                                                                                                                   >  1:452008/1‑229 (MQ=255)                                                      aaGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAAGGAATCA‑G‑G‑TTTTTTCCGGTTAAACCGGCTTCCATCAAGTCCGGGTAGAGAATTGCCCCCCCCGTg                                                                                                                                                                                       >  1:745989/1‑238 (MQ=255)                                                          gTACCGAGAGCGTGGTAGACGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCg                                                                                                                                                                                    <  1:523010/237‑1 (MQ=255)                                                                  agCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGGGAATCC‑G‑TGTTTTTTCCGTTTAACCGGGC‑TCAATCAGTGACGGGTATAGGTTTTGCACCGCCCTGccccgcgtcccac                                                                                                                                                                          >  2:469356/1‑228 (MQ=255)                                                                      tGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTAAGAGGATCA‑GTT‑TTTTTTCCGGTTAA‑CCGGGTTCAATACGTGCCCGATAGAGATTTTTCACCGCCCTGCCCGAGGGCACAATCa                                                                                                                                                                       >  1:464563/1‑238 (MQ=255)                                                                                tAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTtct                                                                                                                                                                                                                                          <  1:748892/161‑1 (MQ=255)                                                                                      ccAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGcc                                                                                                                                                                                     <  1:599645/208‑1 (MQ=255)                                                                                          ccAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGCGGGTTAAAGGATTAG‑G‑T‑TTTTTTCGGGTTAAc                                                                                                                                                                                                                                 >  2:748892/1‑161 (MQ=255)                                                                                                tGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCg                                                                                                                                                                    <  1:205159/215‑1 (MQ=255)                                                                                                       atatCGTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑Gcc                                                                                                                                    <  1:277431/238‑1 (MQ=255)                                                                                                        tttCCTGGGTGGGGGTAGGCAATGAATGTCTGTGCGTGGGATGAG‑G‑GGACAACTTTCCCGGCGGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑Gccc                                                                                                                                   <  2:745989/236‑1 (MQ=255)                                                                                                          gctggggTCGGGCGGGTCATTGAAATACTG‑GGGTTGGGGGAGAA‑GGCAAACATCCCCGGCAGGCAAATTGCCTGATTCGCGCCGCGTTCAATTAGCCCCGGCGGGCGATCGGC‑GTTTAAAAAAATCA‑T‑T‑TTTTCCCCGGTTAACCGGGCTTCAATCAGGGCCGGAAAGAGTTTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAg                                                                                                                                <  2:239815/235‑1 (MQ=255)                                                                                                                  cGGGCGTGGAAAGTAATTCCTGTGC‑TGGGGGTAG‑AGGGCAAACATGCCCGGCTGGCGAATTGCCTGTTCCGCCCCGCAGTCAATTGGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGTTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGa                                                                                                                         <  2:464563/238‑1 (MQ=255)                                                                                                                      cgcgGCAATGAATTGCTGGGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTTATCCGCACCGAAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCCCCGGTTACCCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGc‑c                                                                                                                    <  2:452008/238‑1 (MQ=255)                                                                                                                              tGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑Gccc                                                                                                                                   <  1:511392/216‑1 (MQ=255)                                                                                                                                      tgtgCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTTTCCTGGTAAACCCGCGTCAATCAGTGCCGGATAGAGATTTTTCACCCCCCTGCCGGAGGGGACTAACACCCGCACCCGCCGGCTGCGTACCTTTTGCC‑‑GCCCCCGCGCCCGATACACCCAAcccccagcattg                                                                                                    >  2:447646/1‑230 (MQ=255)                                                                                                                                       gtgCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGGTTAGAGAATCA‑G‑G‑TTTTTTTCGGTTTAACCCGCTTCAATCAGTGGCGGGTAGAGATTTTGCCCCGCCGTGCCGGCGGGCCGCAACACCGCCACCGGCCGCCTGT‑TCGCTTTTCCCACGCCCCCGCCCCAGATCC‑CCACAcccccccgcctcc                                                                                                   >  2:685081/1‑229 (MQ=255)                                                                                                                                         gCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGaaa                                                                                                <  1:416420/239‑1 (MQ=255)                                                                                                                                          cGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTtct                                                                                                                                                                                                                                          >  2:123631/1‑103 (MQ=255)                                                                                                                                          cGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTtct                                                                                                                                                                                                                                          <  1:123631/103‑1 (MQ=255)                                                                                                                                                gtaag‑a‑gGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGAAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGgtg                                                                                          <  1:469356/235‑1 (MQ=255)                                                                                                                                                          cGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGtgtg                                                                                        <  1:561192/232‑1 (MQ=255)                                                                                                                                                           gAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAAcc                                                                                                                                                                                                                                >  1:127437/1‑98 (MQ=255)                                                                                                                                                           gAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAAcc                                                                                                                                                                                                                                <  2:127437/98‑1 (MQ=255)                                                                                                                                                            aaCATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTTTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAgaga                                                                                                                                                                                                      >  1:520003/1‑123 (MQ=255)                                                                                                                                                            aaCATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAgaga                                                                                                                                                                                                      <  2:520003/123‑1 (MQ=255)                                                                                                                                                                ttCCGGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCg                                                                           <  1:447646/239‑1 (MQ=255)                                                                                                                                                                                  ccTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCt                                                         <  1:685902/239‑1 (MQ=255)                                                                                                                                                                                       ttCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGt                                                    <  1:685081/239‑1 (MQ=255)                                                                                                                                                                                                  gTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCt                                                                                                                                                                                                                           <  1:9903/64‑1 (MQ=255)                                                                                                                                                                                                     aaTTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACg                                                                                                          <  1:758322/171‑1 (MQ=255)                                                                                                                                                                                                            tCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑G‑TTTTTTCCGGGTAAACCGGCTTCAATCAggg                                                                                                                                                                                                                 >  2:9903/1‑62 (MQ=255)                                                                                                                                                                                                               ggcgggcgATCGGC‑GGTTAAGAGAATTA‑G‑T‑TTTTTTCCGGTTAAACCGGCTTCAATCAGTGCCCGATAAAGAATTTCCACCCCCGGGCCGGGGGGCACAAAAACCCCCACCGGCCTCTTTG‑TGACTTTTTCCC‑ACCCCAGCGCCAGAAGC‑CCCCACCCACGCCCATCGACATAGGGGTGGTGGATGAAAGCGGAAATTTTCGCCGCCCCCAACGGTAACGGTGTAGAActaggacctcg                            >  2:874088/1‑229 (MQ=255)                                                                                                                                                                                                                         cggc‑ggTTAAGGGAATCA‑G‑T‑TTTCCTCCGGTTACCCCGGCTTCAATCAGTGCGGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAg                  <  1:874088/239‑1 (MQ=255)                                                                                                                                                                                                                                         tCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCt                                                                                                        >  1:684125/1‑138 (MQ=255)                                                                                                                                                                                                                                         tCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCt                                                                                                        <  2:684125/138‑1 (MQ=255)                                                                                                                                                                                                                                          cA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGGGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGGGTTGAACGCGAGCCTGGGGCGATCCAGATGTGTCTGACGCCGt  >  2:937378/1‑239 (MQ=255)                                                                                                                                                                                                                                                    tctcCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCa                                                                                                                                                                 >  2:16779/1‑76 (MQ=255)                                                                                                                                                                                                                                                    tctcCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCa                                                                                                                                                                 <  1:16779/76‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          TCGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAG‑A‑GGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATTAGCTCCGGCGGGCGATCGGC‑GGTTAAGAGAATCA‑G‑T‑TTTTCTCCGGTTAACCCGGCTTCAATCAGTGCCGGATAGAGATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTGC‑TGACTTTCGCCA‑GCCCCAGCGCCAGATGC‑CCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGT  >  NC_000913/2378730‑2379198

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