Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,447,320	T→C	T52A (ACA→GCA)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,447,320	0	T	C	100.0%	163.2 / NA	46	T52A (ACA→GCA)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/17);  total (29/17)

CTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGTCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTT  >  NC_000913/2447187‑2447451                                                                                                                                      |                                                                                                                                    cTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAg                                                                                                                                >  2:1190749/1‑139 (MQ=255)    tAATCCTGAGAGCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCAGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGc                                                                                                                             <  1:1949668/139‑1 (MQ=255)        ccTGAGAGCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGcgacg                                                                                                                         <  1:955441/139‑1 (MQ=255)              agCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGccc                                                                                                                                   >  2:2019862/1‑123 (MQ=255)              agCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGccc                                                                                                                                   <  1:2019862/123‑1 (MQ=255)              agCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTc                                                                                                                   >  2:1756823/1‑139 (MQ=255)                  cTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAgg                                                                                                               >  1:1902154/1‑139 (MQ=255)                   tGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGt                                                                                                              >  2:703636/1‑139 (MQ=255)                   tGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGt                                                                                                              >  1:996965/1‑139 (MQ=255)                        aGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGt                                                                                                         >  1:918789/1‑139 (MQ=255)                               tCGATCAACAAGCCAATGCTGGTGCCGGTTGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGt                                                                                                  <  2:1391035/139‑1 (MQ=255)                                    caacaaGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCt                                                                                             >  2:171511/1‑139 (MQ=255)                                                   ggTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATg                                                                              <  2:996965/139‑1 (MQ=255)                                                    gTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGc                                                                             >  1:22512/1‑139 (MQ=255)                                                    gTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGc                                                                             >  1:1713405/1‑139 (MQ=255)                                                     tGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGcc                                                                            <  1:1756823/139‑1 (MQ=255)                                                         ggTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGa                                                                                 <  2:1595376/130‑1 (MQ=255)                                                         ggTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGa                                                                                 >  1:1595376/1‑130 (MQ=255)                                                         ggTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGc                                                                        <  2:1902154/139‑1 (MQ=255)                                                              taACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGGTCCGGATCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCAGGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAac                                                                   >  2:1375617/1‑139 (MQ=255)                                                                     ttCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTg                                                                                                       >  1:763935/1‑96 (MQ=255)                                                                     ttCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTg                                                                                                       <  2:763935/96‑1 (MQ=255)                                                                      tCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATcga                                                           >  2:577645/1‑139 (MQ=255)                                                                              acCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAg                                                   >  1:1284049/1‑139 (MQ=255)                                                                              acCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAg                                                   >  1:1221269/1‑139 (MQ=255)                                                                                cGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGcc                                                                            <  2:89878/112‑1 (MQ=255)                                                                                cGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGcc                                                                            >  1:89878/1‑112 (MQ=255)                                                                                cGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGcc                                                 >  1:1535787/1‑139 (MQ=255)                                                                                cGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGcc                                                 <  2:1713405/139‑1 (MQ=255)                                                                                             gACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCg                                    >  1:1756587/1‑139 (MQ=255)                                                                                                             cggcgCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCg                                                                                                                  >  1:382425/1‑45 (MQ=255)                                                                                                             cggcgCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCg                                                                                                                  <  2:382425/45‑1 (MQ=255)                                                                                                             cggcgCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAgg                    >  1:39535/1‑139 (MQ=255)                                                                                                             cggcgCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAgg                    >  1:261678/1‑139 (MQ=255)                                                                                                             cggcgCTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAgg                    >  2:949011/1‑139 (MQ=255)                                                                                                                  cTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGtt               >  1:199111/1‑139 (MQ=255)                                                                                                                  cTGGGTGGTATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGtt               <  2:918789/139‑1 (MQ=255)                                                                                                                   tGGGTGGTATAGCGCGATGCCCCAGGGCGACGACgg                                                                                                                     <  2:595524/36‑1 (MQ=255)                                                                                                                   tGGGTGGTATAGCGCGATGCCCCAGGGCGACGACgg                                                                                                                     >  1:595524/1‑36 (MQ=255)                                                                                                                     ggtggtATAGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTAcg            <  1:1190749/139‑1 (MQ=255)                                                                                                                          tataGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGaaa       >  2:1898591/1‑139 (MQ=255)                                                                                                                             aGCGCGATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAagag    <  1:1815691/139‑1 (MQ=255)                                                                                                                               cgcgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAagag    <  1:1326036/137‑1 (MQ=255)                                                                                                                               cgcgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAagag    >  2:1326036/1‑137 (MQ=255)                                                                                                                               cgcgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGtt  >  1:1244986/1‑139 (MQ=255)                                                                                                                               cgcgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGtt  >  2:205820/1‑139 (MQ=255)                                                                                                                                 cgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTa              <  2:1691486/125‑1 (MQ=255)                                                                                                                                 cgATGCCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTa              >  1:1691486/1‑125 (MQ=255)                                                                                                                                      |                                                                                                                                    CTGTAATCCTGAGAGCGCTGGTCAGTGTTTTCGATCAACAAGCCAATGCTGGTGCCGGTAGTAACGCCTTCAAAAACACCGGAGAGAATTTTGACCTGATCCGGCTCGCGGCGCTGGGTGGTATAGCGCGATGTCCCAGGGCGACGACGGTCGAGGTCATGTTGCAGGTCCGCTTCCGTCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTT  >  NC_000913/2447187‑2447451

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