Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,372,049	T→C	I342T (ATC→ACC)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,372,049	0	T	C	100.0%	159.1 / NA	46	I342T (ATC→ACC)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/22);  total (24/22)

ACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAC  >  NC_000913/4371971‑4372126                                                                               |                                                                              aCCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAAccc                                                                             <  1:831300/81‑1 (MQ=255) aCCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAAccc                                                                             <  2:831298/81‑1 (MQ=255)   cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAg                                                                           >  1:831374/1‑81 (MQ=255)   cTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAg                                                                           >  1:831365/1‑81 (MQ=255)    tGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAgg                                                                          >  1:831355/1‑81 (MQ=255)     gggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAggg                                                                         >  1:831363/1‑81 (MQ=255)      ggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGc                                                                        <  1:831302/81‑1 (MQ=255)      ggTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGc                                                                        <  2:831300/81‑1 (MQ=255)       gTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGcc                                                                       >  2:831350/1‑81 (MQ=255)         cAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCgt                                                                     >  2:831330/1‑81 (MQ=255)          aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCgt                                                                     <  2:831303/80‑1 (MQ=255)          aGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCgt                                                                     <  2:831302/80‑1 (MQ=255)            gCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGtt                                                                  >  2:831369/1‑81 (MQ=255)              tAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGc                                                                >  1:831369/1‑81 (MQ=255)               aaaCGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCt                                                               >  2:831360/1‑81 (MQ=255)                  cGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAg                                                            <  1:831305/81‑1 (MQ=255)                    tgttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGAt                                                          <  2:831305/81‑1 (MQ=255)                     gttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGACGCTGCAACCCAGGGCCGTGTTTCTCATATc                                                         >  2:831359/1‑81 (MQ=255)                      ttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATcc                                                        >  2:831379/1‑81 (MQ=255)                      ttgtGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATcc                                                        <  1:831307/81‑1 (MQ=255)                       tgtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCg                                                       >  2:831355/1‑81 (MQ=255)                        gtgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGt                                                      >  1:831367/1‑81 (MQ=255)                         tgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTc                                                     <  2:831308/81‑1 (MQ=255)                         tgATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTc                                                     <  1:831308/81‑1 (MQ=255)                            tCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTcagc                                                  <  1:831310/81‑1 (MQ=255)                              aacaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTcag                                                   >  1:831373/1‑78 (MQ=255)                               acaaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGa                                               <  1:831311/81‑1 (MQ=255)                                caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGAt                                              <  2:831311/81‑1 (MQ=255)                                caaAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGAt                                              >  2:831371/1‑81 (MQ=255)                                 aaaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGAtt                                             >  2:831329/1‑81 (MQ=255)                                 aaaGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGAtt                                             <  1:831313/81‑1 (MQ=255)                                        ccaccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGaagaag                                      >  2:831385/1‑81 (MQ=255)                                         caccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGaagaag                                      >  1:831384/1‑80 (MQ=255)                                          accacTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGc                                     >  1:831390/1‑80 (MQ=255)                                           ccacTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCaa                                   <  1:831314/81‑1 (MQ=255)                                              cTATCATCGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACtt                                <  1:831315/81‑1 (MQ=255)                                                 tcatcGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTc                               <  1:831317/79‑1 (MQ=255)                                                 tcatcGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCt                              <  1:831316/80‑1 (MQ=255)                                                   atcGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                           <  1:831318/81‑1 (MQ=255)                                                   atcGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGAc                           >  2:831377/1‑81 (MQ=255)                                                    tcGATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACt                          <  2:831318/81‑1 (MQ=255)                                                      gATGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTAc                        >  2:831376/1‑81 (MQ=255)                                                       aTGGCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACg                       <  2:831319/81‑1 (MQ=255)                                                         ggCGTGGGTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGAc                     >  2:831386/1‑81 (MQ=255)                                                               ggTGAAGAAGCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGa                <  1:831321/80‑1 (MQ=255)                                                                       agCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTg         >  2:831370/1‑79 (MQ=255)                                                                        gCTGCAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAg      >  1:831394/1‑81 (MQ=255)                                                                            cAACCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAc  <  2:831321/81‑1 (MQ=255)                                                                               |                                                                              ACCTGGGTCAGGCTAAACGTGTTGTGATCAACAAAGACACCACCACTATCATCGATGGCGTGGGTGAAGAAGCTGCAATCCAGGGCCGTGTTGCTCAGATCCGTCAGCAGATTGAAGAAGCAACTTCTGACTACGACCGTGAAAAACTGCAGGAAC  >  NC_000913/4371971‑4372126

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