Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	236,552	C→A	H162Q (CAC→CAA)	dnaQ →	DNA polymerase III subunit epsilon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	236,552	0	C	A	100.0%	114.0 / NA	33	H162Q (CAC→CAA)	dnaQ	DNA polymerase III subunit epsilon
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/12);  total (21/12)

ATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTT  >  NC_000913/236475‑236626                                                                              |                                                                           aTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAggg                                                                         <  1:47097/81‑1 (MQ=255)  tGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAgggg                                                                        >  1:47128/1‑81 (MQ=255)   gTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAgggg                                                                        >  1:47143/1‑80 (MQ=255)    tttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCa                                                                      >  1:47129/1‑81 (MQ=255)       cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCAtt                                                                    <  2:47098/80‑1 (MQ=255)       cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCAtt                                                                    <  1:47099/80‑1 (MQ=255)       cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTa                                                                   >  2:47138/1‑81 (MQ=255)        ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTAc                                                                  >  1:47135/1‑81 (MQ=255)           gTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTc                                                                >  1:47134/1‑80 (MQ=255)              aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATg                                                            >  1:47140/1‑81 (MQ=255)                  cAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGcc                                                          >  2:47141/1‑79 (MQ=255)                       gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATcc                                                   <  1:47100/81‑1 (MQ=255)                        ccTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCt                                                  >  1:47148/1‑81 (MQ=255)                            gATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTg                                                <  1:47101/79‑1 (MQ=255)                            gATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCg                                              >  2:47156/1‑81 (MQ=255)                            gATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCg                                              >  2:47142/1‑81 (MQ=255)                                      gtgCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATc                                    >  2:47159/1‑81 (MQ=255)                                         cTCGCTACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTg                                  <  1:47102/80‑1 (MQ=255)                                              tACGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGAt                             >  1:47158/1‑80 (MQ=255)                                                cGAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                            <  2:47103/79‑1 (MQ=255)                                                 gAAATAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGAcc                         >  2:47150/1‑81 (MQ=255)                                                     tAGATAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtg                     <  1:47104/81‑1 (MQ=255)                                                        aTAACAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                  <  2:47105/81‑1 (MQ=255)                                                          aaCAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCa                 >  1:47157/1‑80 (MQ=255)                                                          aaCAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCa                 >  1:47147/1‑80 (MQ=255)                                                            cAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc              <  2:47107/81‑1 (MQ=255)                                                            cAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc              <  2:47108/81‑1 (MQ=255)                                                            cAGTAAACGAACGCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc              <  1:47106/81‑1 (MQ=255)                                                                 aacgaacgCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa         >  2:47152/1‑81 (MQ=255)                                                                    gaacgCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg      >  1:47164/1‑81 (MQ=255)                                                                    gaacgCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg      >  1:47144/1‑81 (MQ=255)                                                                      acgCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCt    >  1:47155/1‑81 (MQ=255)                                                                        gCTGCAAGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCttt  >  1:47149/1‑81 (MQ=255)                                                                              |                                                                           ATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTT  >  NC_000913/236475‑236626

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