Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,807,326	T→C	D258G (GAT→GGT)	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,326	0	T	C	100.0%	118.7 / NA	33	D258G (GAT→GGT)	waaQ	lipopolysaccharide core heptosyltransferase 3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/11);  total (22/11)

GAATACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAATCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTA  >  NC_000913/3807250‑3807403                                                                             |                                                                              gAATACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCcac                                                                           <  2:740079/81‑1 (MQ=255) gAATACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCca                                                                            >  2:740127/1‑80 (MQ=255)  aaTACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCcaca                                                                          >  2:740121/1‑81 (MQ=255)  aaTACGTGATCCATTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCcaca                                                                          >  1:740112/1‑81 (MQ=255)    tACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACAcc                                                                        <  1:740080/81‑1 (MQ=255)    tACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACAcc                                                                        <  2:740081/81‑1 (MQ=255)         gATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATa                                                                    >  1:740107/1‑80 (MQ=255)                 gCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACt                                                           >  1:740103/1‑81 (MQ=255)                          aGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATc                                                   >  2:740116/1‑80 (MQ=255)                           gACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCaa                                                 >  1:740115/1‑81 (MQ=255)                            aCTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAAt                                                >  1:740123/1‑81 (MQ=255)                              tAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAAtt                                               >  1:740109/1‑80 (MQ=255)                                  gACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATg                                           >  1:740131/1‑80 (MQ=255)                                   ccTGGCGTTTTCACTGCCGCTGCCCTCTTTCCCGGCGCCGCCCCCACACCAATAAAAAGCACTGCATGATCAATTCATGCa                                         >  1:740117/2‑81 (MQ=255)                                    cTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCAc                                        >  2:740147/1‑81 (MQ=255)                                    cTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCAc                                        >  1:740125/1‑81 (MQ=255)                                          ttttCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGtt                                  <  1:740082/81‑1 (MQ=255)                                           tttCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGtt                                  >  1:740110/1‑80 (MQ=255)                                            ttCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGtt                                  >  2:740111/1‑79 (MQ=255)                                            ttCACTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTc                                 <  2:740083/80‑1 (MQ=255)                                               aCTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAgg                              >  2:740114/1‑80 (MQ=255)                                                cTGCCGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGa                             >  2:740128/1‑80 (MQ=255)                                                   ccGCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAAc                         >  1:740113/1‑81 (MQ=255)                                                     gCTGCAATATGTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACgt                       >  1:740118/1‑81 (MQ=255)                                                               gTCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTg             >  1:740119/1‑81 (MQ=255)                                                                tCCCGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGc            <  1:740084/81‑1 (MQ=255)                                                                 cccGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGc            >  1:740124/1‑80 (MQ=255)                                                                   cGGCGCAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAg         <  1:740085/81‑1 (MQ=255)                                                                        cAGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAg    <  1:740086/81‑1 (MQ=255)                                                                         aGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAg    >  2:740138/1‑80 (MQ=255)                                                                         aGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGt   <  1:740088/81‑1 (MQ=255)                                                                         aGAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGt   <  1:740087/81‑1 (MQ=255)                                                                          gAACCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTa  <  2:740089/81‑1 (MQ=255)                                                                             |                                                                              GAATACGTGATCCGTTGCACCAAATAGACTAATGACTGGCGTTTTCACTGCCGCTGCAATATGTCCCGGCGCAGAATCCACACCAATAAAAAGCACTGCATGATCAATTAATGCACCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTA  >  NC_000913/3807250‑3807403

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