Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,383,639	T→C	24.4%	G97G (GGA→GGG)	menD ←	2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,383,639	0	T	C	24.4%	120.7 / 11.2	67	G97G (GGA→GGG)	menD	2‑succinyl‑5‑enolpyruvyl‑6‑hydroxy‑3‑ cyclohexene‑1‑carboxylic‑acid synthase
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base C (2/15);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGC  >  NC_007779/2383401‑2383873                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               cGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTtct                                                                                                                                                                                                                                                <  1:533355/239‑1 (MQ=255) cGCCAGTCAGCGGTTAAGGAAAGGGGCAGTTGATATGGACTCACCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTtct                                                                                                                                                                                                                                                <  2:918032/239‑1 (MQ=255)      tACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑t                                                                                                                                                                                                                                          <  1:980199/239‑1 (MQ=255)         aGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑tt                                                                                                                                                                                                                                         <  1:893233/237‑1 (MQ=255)           cGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGGCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCGGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGGGGA‑G‑CC‑GTTTTTTTccccgtcag                                                                                                                                                                                                                                       >  1:1102416/1‑231 (MQ=255)            ggTTCAGCAAACGGGCAGTTGATATTGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAAcc                                                                                                                                                                                                                                     <  2:2193182/238‑1 (MQ=255)             gTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTG‑AGAGAA‑T‑GAGGTTTTTTTCGGG‑TTAcccc                                                                                                                                                                                                                                    >  1:1104651/1‑239 (MQ=255)              ttCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTATAAGGAA‑T‑AG‑GTTTTTTCCCGT‑TTAACCCg                                                                                                                                                                                                                                   >  1:110870/1‑239 (MQ=255)                cAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTT‑AGTAGA‑T‑CA‑GTTTTTTTCCCG‑GTTACCCGG‑Ct                                                                                                                                                                                                                               >  1:1330055/1‑239 (MQ=255)                cAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTT‑CGGGAA‑G‑CC‑GTTTTTTTTCGGCTGAACCCgg                                                                                                                                                                                                                                  >  1:612371/1‑238 (MQ=255)                cAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGCGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGGGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTG‑AGGGAA‑TACG‑GTTTTTTTCCGT‑TAAACACg                                                                                                                                                                                                                                   >  1:1386736/1‑237 (MQ=255)                   cAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTT‑AGAGAAGA‑AA‑GTTTTTTCCGGG‑GTAACCCgcttcc                                                                                                                                                                                                                              >  1:329593/1‑234 (MQ=255)                     aaCGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAAt                                                                                                                                                                                                                          <  1:2219004/239‑1 (MQ=255)                        gggCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTA‑AGAGAA‑G‑CA‑GATTTTCTCCGG‑GTTACCCGG‑CCTCCAtaac                                                                                                                                                                                                                       >  1:1221077/1‑236 (MQ=255)                                gttatGGACTCCCCCCGCATGAAGCGTGCCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  2:2161513/237‑1 (MQ=255)                                gATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  1:627382/239‑1 (MQ=255)                                gATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  1:598762/239‑1 (MQ=255)                                gATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  1:739011/239‑1 (MQ=255)                                gATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  2:1715607/239‑1 (MQ=255)                                  tatGGCCTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATag                                                                                                                                                                                                             <  1:2032513/239‑1 (MQ=255)                                  tatGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAAATGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGATA‑GGAGGA‑G‑CA‑GTTTTTTTCCGT‑TTAACCCGG‑CGTCCACAAGTGCCGGATaa                                                                                                                                                                                                             >  1:530744/1‑238 (MQ=255)                                           cccccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTTTCCCG‑GTAACCCCG‑CCTCAATCAGTGCCGGATAGA‑GATTTGc                                                                                                                                                                                                    >  1:2017859/1‑238 (MQ=255)                                           cccccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTGA‑AGAGGA‑G‑CA‑GGTTTTTTTCTG‑CTAAACCGG‑CTGTCATCAGTGCCGGGGAGA‑GAatttct                                                                                                                                                                                                   >  1:61627/1‑236 (MQ=255)                                           cccccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTA‑AGAGAA‑A‑CA‑GGTTTTTTCCGG‑GTAACCCCG‑GCTCGATCAGAGCCGGGTAGA‑GAttttg                                                                                                                                                                                                    >  1:11339/1‑236 (MQ=255)                                           cccccGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTA‑AAAGGA‑T‑GC‑GTTTTTTTCCCG‑GTAATCCGG‑GCGCCATCAGTGGCGGGAAGA‑GATTTGc                                                                                                                                                                                                    >  1:674278/1‑238 (MQ=255)                                                   tGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGGTA‑AGAGAA‑G‑CG‑GTTTTTTTCCGG‑GTAACCCGG‑CTTCACTACGTGCCGGGAAGA‑GAATTGCCACCGCCGt                                                                                                                                                                                           >  1:5278/1‑239 (MQ=255)                                                     aaGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑c                                                                                                                                                                                                                                <  1:1562960/201‑1 (MQ=255)                                                     aaGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGACGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGCGGGTTA‑GGGAAA‑C‑CA‑GTTTTTTTCCGG‑TTACCCCGC‑CTCCAA‑CACGGCCGGGTAGAGGATTTGcacctcccgg                                                                                                                                                                                           >  1:160787/1‑229 (MQ=255)                                                       gCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTT‑AGAGAA‑A‑CA‑GGTTTTCTTCGG‑TTAACCCGG‑CTTCAATCAGGGGCGGGTAGA‑AAATTTCCACCCCCGTGCCg                                                                                                                                                                                       >  1:1852000/1‑239 (MQ=255)                                                             ggagtGCGTGGTGGGGGGTAACAGCCCGCCAGCGTGCGGGGATTTCTGGGGTCGGGCGGGGAAATTAAATACGTTGGGTGGGGGGAGGGGCGAACATTCCCGGCTGGGGATTTGCCGGATTCGCGCCGCAGTAAGTT‑AGCTCCGGCGGGCGAGCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCGGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGt                                                                                                                                                                                  <  2:1813867/237‑1 (MQ=255)                                                              gagaGCGTGGTCGATGGGAGAAACCAGCCAGCGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTca                                                                                                                                                                                <  1:863622/239‑1 (MQ=255)                                                                              gTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGAt                                                                                                                                                                                                               <  1:1844117/193‑1 (MQ=255)                                                                                              gTGGGGGGATATACTGGGTCGGGCGCGGAAATGAAATTCTGGGCGGGGGGTGAGAGGCGAACATCCCCGGCTGGGGAATTGCCTGATTCGCGCCGCATTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGAATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACtt                                                                                                                                                 <  2:623024/237‑1 (MQ=255)                                                                                                         tCCGGGGTGGGGCGAGGCAATGAACTCCGGTGGGTGGGGGGGGGGGCGACCTTCCCCGGCGGGGGAATTGCCTGATTCGCGCCGCGGTCAATT‑GGCTCCGGGGGGCGATCGGG‑GGTTA‑AGAGAT‑T‑CG‑TTTTTTCCCCGT‑TTACCCCGG‑CTTCAATCAGGGCGGGATGGG‑GATTGCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAg                                                                                                                                    <  2:2031460/239‑1 (MQ=255)                                                                                                           cggggTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCGGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTGGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAgcg                                                                                                                                  <  1:395688/237‑1 (MQ=255)                                                                                                              ggTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGGA‑T‑CA‑GTTTTTTTCCGG‑GTAACCCCG‑GTTCAAACCGTGCCGGGTAGA‑GAGTTTCCACCGCCGTGCCGGGGGGCACAAACACCCCCACCCGCCGCTTT‑GTTGCTTTTGCCCGCCCCCGCGcc                                                                                                                                >  1:1228926/1‑238 (MQ=255)                                                                                                                               aaagccTGTGCGTGGGGTGAGAGGCGAACATCCCCGGCTGGCGAATTGCCTGATCCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc                                                                                                              <  1:952743/234‑1 (MQ=255)                                                                                                                               aaaTACTGGGCGGGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATGGCCTGATTCGCGCCGCAGTCAATT‑AGCTCGGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGcccc                                                                                                                       <  1:1422132/230‑1 (MQ=255)                                                                                                                                aaTACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTAATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTCCTCGGG‑TTAACCCGG‑CTTCAACCAGGGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACg                                                                                                               <  1:357153/237‑1 (MQ=255)                                                                                                                                      gtgCGGGGGGTGGGAGGCGAACATCCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTACCCCGG‑CTTCAACCAGTGCCGGATAGA‑GTTTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGa                                                                                                       <  1:1951294/239‑1 (MQ=255)                                                                                                                                       tgCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGCGGGTTA‑AG‑GAA‑A‑CA‑GGTTTTTTCCCG‑GTTACCCCG‑CCTCAATCAGTGTCGGGTAGA‑GAATTTCCCCCCCCGTGCCCGGGGGCACAAACACCCCCCCCCGCCGCCTGTCTGCCTTT‑GCCCGCCCCCGCCCCCCACTCCCCAACCCCCCCTCCTCCaa                                                                                                      >  2:444127/1‑239 (MQ=255)                                                                                                                                             gggTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATggg                                                                                                 <  1:92346/238‑1 (MQ=255)                                                                                                                                               gTGAGGGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGgtg                                                                                               <  1:2058646/238‑1 (MQ=255)                                                                                                                                                         tccATCCCCGGGTGGCGAAGTGCCGTGTTCGGCCCGCATTAGATT‑AGCACGGGCGGGGGCTCGGT‑GGTTG‑AGAAAA‑C‑AA‑TTTTTTCTCCGG‑TTACCCCGC‑CTTCAATC‑GTGCCGGAATGAGGTTTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGTCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGtgaatga                                                                                     <  2:1221077/236‑1 (MQ=255)                                                                                                                                                                ccGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGaa                                                                             <  1:394676/239‑1 (MQ=255)                                                                                                                                                                              tccTGGTCCGCGCCGGTGTCATTC‑ACCTCCGGGGGGCGCTGGGC‑GTTTA‑AGAGAA‑T‑CG‑GTTTTTCCCCTG‑TTACCCCGG‑CTTAAAACAGCCCCGGATAGA‑GATTTCCCCCCGCCGTGCCGGAGTTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTccgccgccg                                                                <  2:886348/237‑1 (MQ=255)                                                                                                                                                                               ctgtgtgggcgcCGGAG‑CAATTGAGCGCCGCGGGGGGAGCGGG‑GGTAA‑AGAAAA‑T‑CA‑‑TTTTTCCCCCG‑TTTACCCGGTTTTCAACCGGGGCGGGAGAGA‑GTTTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCt                                                              <  2:1228926/231‑1 (MQ=255)                                                                                                                                                                                cTGATTCGCGCCGCAGTCAATT‑GGCTCCGGCGGGCGATGGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCGGG‑TTACCCCGG‑CTTCAACCAGTGCCGGATAGA‑GATTTCCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCt                                                              <  1:547531/238‑1 (MQ=255)                                                                                                                                                                                cTGATTCGCGCCGCAGTAAATT‑AGCTCGGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCt                                                              <  1:1635177/238‑1 (MQ=255)                                                                                                                                                                                    ttCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAT‑C‑CA‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGGCCCGGAGAGA‑GTTTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGt                                                         <  1:89761/239‑1 (MQ=255)                                                                                                                                                                                       gcgcCGCAGTCAATT‑AGCCCCGGGGGGCGATCGGC‑GGTTA‑AGAGAA‑C‑CA‑GTTTTCCTCCGG‑TTACCCCGG‑CTTCAACCAGTGCGGGATAGA‑GTTTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTaa                                                      <  1:182428/239‑1 (MQ=255)                                                                                                                                                                                       gccccgcAGTCAATT‑CGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGTGCGGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTaa                                                      <  1:1749694/236‑1 (MQ=255)                                                                                                                                                                                              attCAATT‑AGCGCGGGCGGGCGATCGGC‑GTTAA‑AAAAAA‑T‑CG‑GTTTTTCCCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTa                                               <  1:26843/237‑1 (MQ=255)                                                                                                                                                                                                    tttcgcccggccgggcgATGGGG‑GGTAA‑AGAGAA‑T‑CA‑GTTTTGCCCCGG‑TTAACCCGG‑TTCCAATCAGGGCCGGAAAGA‑GATTTGCCGCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAAcg                                          <  1:1335165/228‑1 (MQ=255)                                                                                                                                                                                                          tCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑TTTTTTCCCCGT‑TTACCCCGC‑CTTCAATCAGTGCCGGATAGA‑GATTTCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGcc                                    <  1:664414/239‑1 (MQ=255)                                                                                                                                                                                                          tCCGGCGGGCGACCGGC‑GTTTA‑AGAGAA‑T‑CA‑GTTTTTCCCCGG‑TTACCCCGC‑CTTCATTCAGTCCCGGAAAGA‑GTTTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGcc                                    <  1:1375916/239‑1 (MQ=255)                                                                                                                                                                                                            gggcgggcgACCGGC‑GTTTA‑AGAGAA‑T‑CA‑TTTTTTTCCCGT‑TTACCCCGC‑CTTCAATGGGTGCCGAATAGA‑GATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTg                                  <  1:1782890/238‑1 (MQ=255)                                                                                                                                                                                                                    gATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGT‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGAt                           <  1:45778/238‑1 (MQ=255)                                                                                                                                                                                                                       gggc‑ggTAA‑AGAGAA‑T‑CA‑TTTTTTCCCCGG‑TTAGCCCGG‑CTTCACTCAGTGCGGGATAGA‑GATTTGCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAg                       <  1:275640/238‑1 (MQ=255)                                                                                                                                                                                                                       cggc‑ggTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTACCCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAg                       <  1:444127/239‑1 (MQ=255)                                                                                                                                                                                                                                      aa‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGACACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGGGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGTCg           >  1:1936989/1‑238 (MQ=255)                                                                                                                                                                                                                                       a‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc                                                                                                              >  2:150379/1‑138 (MQ=255)                                                                                                                                                                                                                                       a‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGc                                                                                                              <  1:150379/138‑1 (MQ=255)                                                                                                                                                                                                                                         t‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGa                                             >  1:1891741/1‑202 (MQ=255)                                                                                                                                                                                                                                              ttttttCCCCGT‑TTACCCCGG‑CTTCAACCAGTGCGGGAAAGA‑GATTGCCCCCCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGa     <  1:2072147/238‑1 (MQ=255)                                                                                                                                                                                                                                              gTTTTTCTCCGG‑TATACCCGG‑CTCAACCCAGTGCCGGATGGA‑GTTTTGCCACCGCCGTGCCGGAGGTCACAATCCCCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCGAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGa     <  1:1728721/239‑1 (MQ=255)                                                                                                                                                                                                                                               tttttCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGAc    >  2:4080/1‑239 (MQ=255)                                                                                                                                                                                                                                                 tttCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGAcgc  >  2:846127/1‑239 (MQ=255)                                                                                                                                                                                                                                                 tttCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGGGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTCTCTGACGCCGTGAcgc  >  1:1636261/1‑239 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CGCCATACAGCGGTTCAGCAAACGGGCAGTTGATATGGACTCCCCCCGCATGAAGCGTACCGAGAGCGTGGTCGATGGTAGAAACCAGCCAACGTGCGGGGATATCCTGGGTCGGGCGCGGCAATGAAATACTGTGCGTGGGGTGAGAGGCGAACATTCCCGGCTGGCGAATTGCCTGATTCGCGCCGCAGTCAATT‑AGCTCCGGCGGGCGATCGGC‑GGTTA‑AGAGAA‑T‑CA‑GTTTTTCTCCGG‑TTAACCCGG‑CTTCAATCAGTGCCGGATAGA‑GATTTGCCACCGCCGTGCCGGAGGTCACAATCACCGCCACCGGCTGCTTG‑CTGACTTTCGCCAGCCCCAGCGCCAGATGCCCCAACCCACGCTCATCGAAATGGGTGTGGTGAATGAATGCGGAATTCTCCGCCGCCGCTAACGTTAACGGTGTAGAACGCGAGCCTGGGGCGATACAGATGTGTCTGACGCCGTGACGC  >  NC_007779/2383401‑2383873

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