Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,807,402	T→C	T233A (ACT→GCT)	waaQ ←	lipopolysaccharide core heptosyltransferase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,807,402	0	T	C	100.0%	132.9 / NA	38	T233A (ACT→GCT)	waaQ	lipopolysaccharide core heptosyltransferase 3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/19);  total (19/19)

ATCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAG  >  NC_000913/3807325‑3807477                                                                              |                                                                            aTCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAAt                                                                          <  1:681334/81‑1 (MQ=255)  tCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATg                                                                         <  1:681335/81‑1 (MQ=255)       acCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGttt                                                                    <  1:681336/81‑1 (MQ=255)            tAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGttt                                                               <  1:681337/81‑1 (MQ=255)              aaaaGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTc                                                              <  2:681337/80‑1 (MQ=255)               aaaGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCg                                                             >  2:681373/1‑80 (MQ=255)               aaaGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGc                                                            <  1:681339/81‑1 (MQ=255)                     cTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTc                                                      >  2:681390/1‑81 (MQ=255)                              cAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAAtct                                             >  2:681388/1‑81 (MQ=255)                                 ttAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCAt                                          >  2:681385/1‑81 (MQ=255)                                   aaTGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATc                                         <  1:681340/80‑1 (MQ=255)                                    aTGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCt                                        <  1:681341/80‑1 (MQ=255)                                    aTGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTa                                       >  2:681384/1‑81 (MQ=255)                                        aCCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTacac                                    >  1:681397/1‑80 (MQ=255)                                           cAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGc                                >  1:681401/1‑81 (MQ=255)                                               tCAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGc                                >  1:681390/1‑77 (MQ=255)                                                cAGGAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGc                                >  1:681398/1‑76 (MQ=255)                                                  ggAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGatca                         <  1:681342/81‑1 (MQ=255)                                                   gAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGatcat                        >  2:681381/1‑81 (MQ=255)                                                   gAAAACGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGatcat                        <  2:681342/81‑1 (MQ=255)                                                    aaaaCGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGatcatc                       >  2:681406/1‑81 (MQ=255)                                                     aaaCGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCt                      <  2:681343/81‑1 (MQ=255)                                                       aCGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGc                    <  1:681345/81‑1 (MQ=255)                                                        cGTGTTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGcc                   >  2:681408/1‑81 (MQ=255)                                                         gtgtTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGcc                   >  1:681404/1‑80 (MQ=255)                                                          tgtTTTACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGcc                   >  1:681402/1‑79 (MQ=255)                                                             tttACCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAggg              >  2:681392/1‑81 (MQ=255)                                                                aCCTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGcc            >  1:681408/1‑80 (MQ=255)                                                                 ccTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCca          <  2:681345/81‑1 (MQ=255)                                                                 ccTGCAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCca          >  2:681409/1‑81 (MQ=255)                                                                     cAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCcaca        <  1:681348/79‑1 (MQ=255)                                                                     cAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCcaca        <  1:681347/79‑1 (MQ=255)                                                                     cAAGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAgg      >  1:681396/1‑81 (MQ=255)                                                                      aaGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGt     <  2:681348/81‑1 (MQ=255)                                                                       aGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGtt    <  2:681350/81‑1 (MQ=255)                                                                       aGGCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGtt    <  1:681350/81‑1 (MQ=255)                                                                        ggCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGt     >  2:681410/1‑79 (MQ=255)                                                                         gCCAGCAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAg  <  1:681352/81‑1 (MQ=255)                                                                              |                                                                            ATCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAG  >  NC_000913/3807325‑3807477

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