Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,195,773	T→C	E327E (GAA→GAG)	rfaE ←	fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,195,773	0	T	C	100.0%	98.3 / NA	30	E327E (GAA→GAG)	rfaE	fused heptose 7‑phosphate kinase/heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/18);  total (12/18)

GGTCACCCAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCT  >  NC_000913/3195649‑3195908                                                                                                                             |                                                                                                                                        ggTCACCCAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGtcatc                                                                                                                          <  2:530995/140‑1 (MQ=255)    cACCCAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACg                                                                                                                       >  2:547828/1‑140 (MQ=255)       ccAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCa                                                                                                                    <  2:181118/140‑1 (MQ=255)       ccAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCa                                                                                                                    <  2:690081/140‑1 (MQ=255)               cGGGCATTTGCCAGATAGGAGACGTGCCCGGGGTGAAGGTTGTCAAAGTCACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGt                                                                                                            <  2:289948/140‑1 (MQ=255)                gggCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTa                                                                                                           <  1:547828/140‑1 (MQ=255)                         ccAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGt                                                                                                  <  2:76545/140‑1 (MQ=255)                                agagACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTa                                                                                           <  1:680242/140‑1 (MQ=255)                                  agaCGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACa                                                                                         <  2:618057/140‑1 (MQ=255)                                    aCGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGc                                                                                       >  1:84169/1‑140 (MQ=255)                                       tGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCAtt                                                                                    >  1:501336/1‑140 (MQ=255)                                       tGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCAtt                                                                                    >  2:115072/1‑140 (MQ=255)                                       tGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTCTACCTGCACGTCCACGTACAGCAtt                                                                                    >  1:271278/1‑140 (MQ=255)                                        gCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCAttt                                                                                   <  1:115072/140‑1 (MQ=255)                                         cccGGCGTGCAGGATGTCAAATACACCGTTGGTCATGACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAACCCTGTATCTGCACGTCCACGTACAGCAtttt                                                                                  >  2:522501/1‑140 (MQ=255)                                         cccGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCAtttt                                                                                  >  2:793257/1‑140 (MQ=255)                                            ggCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCa                                                                               >  2:601545/1‑140 (MQ=255)                                                 gCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCg                                                                          <  2:501336/140‑1 (MQ=255)                                                   aGGATGTCAAAGACACTTTTGGTCATCACCTCTTTGTCCGCACGTTTCTGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGAt                                                                        <  2:380268/140‑1 (MQ=255)                                                           aaaGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGaaa                                                                <  1:9137/140‑1 (MQ=255)                                                                  ccGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGa                                                         <  1:601545/140‑1 (MQ=255)                                                                                  cTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCAAGCACGATCGGCGAAAGCGTGTAGGTTCCCAGTTTGctg                                         >  2:251072/1‑138 (MQ=255)                                                                                    ttttCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGac                                       <  1:793257/140‑1 (MQ=255)                                                                                     tttCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGacc                                      >  2:535970/1‑140 (MQ=255)                                                                                           cacGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACg                                <  2:84169/140‑1 (MQ=255)                                                                                             cGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAGACCATGGAGGTTACCAGATTGCCGACCACCACGcc                              >  2:234902/1‑140 (MQ=255)                                                                                             cGTTTACGCGCTGCCGCTACGGCCAGCTTCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGcc                              >  2:152536/1‑140 (MQ=255)                                                                                                       cTGCCGCTACGGCCAGCTTCAGCTCCTCTCCGGTCATTACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCa                    <  1:234902/140‑1 (MQ=255)                                                                                                        tGCCGCTACGGCCAGTTTCAGCTCCTCTTCGGTCATCTCGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCTGACCACCACGCCAGCCGCCGCAt                   <  1:251072/140‑1 (MQ=255)                                                                                                                        ttCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAaagaag        >  1:816385/1‑135 (MQ=255)                                                                                                                        ttCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAaagaag        <  2:816385/135‑1 (MQ=255)                                                                                                                         tCAGCTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGGTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGct  <  1:48070/140‑1 (MQ=255)                                                                                                                             |                                                                                                                                        GGTCACCCAGCTTGCGGGCATTTGCCAGATAAGAGACGTGCCCGGCGTGCAGGATGTCAAAGACACCGTTGGTCATCACCACTTTTTCACCACGTTTACGCGCTGCCGCTACGGCCAGCTTCAGTTCCTCTTCGGTCATCACGCCAAAGCCTGTATCTGCACGTCCACGTACAGCATTTTCCAGCTCGATCGGCGAAACCGTGGAGGTTCCCAGTTTGCCGACCACCACGCCAGCCGCCGCATTGGCAAAGAAGCAGGCT  >  NC_000913/3195649‑3195908

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