Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	236,541	C→T	R159* (CGA→TGA)	dnaQ →	DNA polymerase III epsilon subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	236,541	0	C	T	100.0%	230.9 / NA	66	R159* (CGA→TGA)	dnaQ	DNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/33);  total (33/33)

GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236445‑236637                                                                                                 |                                                                                                 gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  2:808889/1‑100 (MQ=255) gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               <  1:1374600/100‑1 (MQ=255)   cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaacg                                                                                             >  1:528627/1‑100 (MQ=255)   cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaacg                                                                                             >  2:1563995/1‑100 (MQ=255)    aCCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGc                                                                                            <  2:1291813/100‑1 (MQ=255)     ccGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCt                                                                                           >  2:1077309/1‑100 (MQ=255)       gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         >  1:339052/1‑100 (MQ=255)        aTAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCa                                                                                        <  1:1507514/100‑1 (MQ=255)            ncTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACggg                                                                                    >  2:1143604/2‑100 (MQ=255)            ccTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACggg                                                                                    <  1:384885/100‑1 (MQ=255)              ttGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGc                                                                                  <  1:589959/100‑1 (MQ=255)                 cGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAACAGATAACAGTAAATGAACGCTGCACGGGGCAtt                                                                               >  1:260519/1‑100 (MQ=255)                  ggTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTa                                                                              >  2:1095412/1‑100 (MQ=255)                  ggTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTa                                                                              >  2:977247/1‑100 (MQ=255)                      gCGAGAAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTcn                                                                          <  2:803369/100‑2 (MQ=255)                       cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  2:971968/1‑100 (MQ=255)                        gAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGAt                                                                        <  2:838338/100‑1 (MQ=255)                        gAGGAAAATGTTTCCCGGTAAGCGCAACAGCATCGATGCGTTATGTGCTCGCTACGAAATAGATACCACGAAATGAACGCTGCACGCGGCATTACTCgtt                                                                        <  1:1064860/100‑3 (MQ=255)                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  2:1273392/100‑1 (MQ=255)                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  2:87642/1‑100 (MQ=255)                          ggAAAATGTTTCCCGGTAAGCGCAACAGCCACGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  2:1041797/1‑100 (MQ=255)                             aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                   <  1:1143609/100‑1 (MQ=255)                                tGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGAt                                                                >  1:212980/1‑100 (MQ=255)                                  tttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATcc                                                              >  1:603111/1‑100 (MQ=255)                                   ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             >  1:999333/1‑100 (MQ=255)                                      ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  2:1425541/100‑1 (MQ=255)                                      ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  2:1425509/100‑1 (MQ=255)                                       cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         <  1:120083/100‑1 (MQ=255)                                       cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         >  2:181624/1‑100 (MQ=255)                                        ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCgg                                                        >  1:861920/1‑100 (MQ=255)                                         gTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGa                                                       <  2:901420/100‑1 (MQ=255)                                         gTAAACGCAACAGCCTCGATGCGTTATGTGTTCGCTATCAACTAGCTAACAGTCAATCGACGCTGCGCGGGGCATTAATCGATGCCCACATCCGTGCgg                                                        <  2:1461191/99‑1 (MQ=255)                                           aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  1:1540510/100‑1 (MQ=255)                                            aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    <  2:861915/100‑1 (MQ=255)                                             gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   >  2:167702/1‑100 (MQ=255)                                                 aaCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATc                                               <  1:87642/100‑1 (MQ=255)                                                   cAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCTTTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTg                                             <  1:1113461/100‑1 (MQ=255)                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  1:1304299/100‑1 (MQ=255)                                                     gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:1413529/100‑1 (MQ=255)                                                        tCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGAt                                        >  1:1515339/1‑100 (MQ=255)                                                         cGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                       <  1:1376408/100‑1 (MQ=255)                                                             gCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCg                                   <  1:276139/100‑1 (MQ=255)                                                               gTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTCCTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggt                                 >  2:1087504/1‑100 (MQ=255)                                                                 tATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                               <  1:701016/100‑1 (MQ=255)                                                                  aTGTGCTCGCTGCGAAATAGATAACAGTAAATGAAAGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                  >  2:879949/1‑96 (MQ=255)                                                                  aTGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                              >  2:613536/1‑100 (MQ=255)                                                                     ngCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                             >  2:730785/2‑98 (MQ=255)                                                                       cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                         <  1:737797/100‑1 (MQ=255)                                                                       cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                         >  2:901473/1‑100 (MQ=255)                                                                        tCGCTACGAAATAGATAACAGTAAATTAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACg                        <  2:784260/100‑1 (MQ=255)                                                                          gCTACGAAATAGATAACAGTAAATGAACGTTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGCAGTTTATCTGGCGATGACCGGCGGTCAAACGTc                      <  2:950425/100‑1 (MQ=255)                                                                            tACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                    <  1:1300691/100‑1 (MQ=255)                                                                             aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                    <  1:405324/99‑1 (MQ=255)                                                                             aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGAt                   >  2:1396068/1‑100 (MQ=255)                                                                              cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  1:1353813/1‑100 (MQ=255)                                                                              cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  2:1310715/1‑100 (MQ=255)                                                                                    aGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCtttt            >  1:1261967/1‑100 (MQ=255)                                                                                    aGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCtttt            <  2:29281/100‑1 (MQ=255)                                                                                     gATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg           >  1:536250/1‑100 (MQ=255)                                                                                     gATAACAGTAAATGAACGCTGCACGGGGCATTAATCGATGCCCAGATCCTTGCGGAAGTATATCTGGCGATGACCGGTGGTCACACGTCGATCGATTTTg           >  1:129966/1‑100 (MQ=255)                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAtn      <  2:260519/100‑2 (MQ=255)                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:663086/1‑100 (MQ=255)                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:301952/1‑100 (MQ=255)                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:926348/1‑100 (MQ=255)                                                                                          cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:1474847/1‑100 (MQ=255)                                                                                           aGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     <  2:536249/100‑1 (MQ=255)                                                                                            gTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa    <  2:189109/100‑1 (MQ=255)                                                                                            gTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa    <  1:167702/100‑1 (MQ=255)                                                                                              aaaTGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAg  >  1:732333/1‑100 (MQ=255)                                                                                              aaaTGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAg  >  1:1472264/1‑100 (MQ=255)                                                                                                 |                                                                                                 GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236445‑236637

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