Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,826,047	A→T	I280F (ATC→TTC)	recG →	ATP‑dependent DNA helicase RecG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,826,047	0	A	T	100.0%	189.9 / NA	52	I280F (ATC→TTC)	recG	ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (17/35);  total (17/35)

AATGACACGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGATCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCTGGTC  >  NC_000913/3825969‑3826124                                                                               |                                                                              aaTGACACGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTc                                                                             <  1:753396/81‑1 (MQ=255) aaTGACACGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTc                                                                             >  2:753467/1‑81 (MQ=255)  aTGACACGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCg                                                                            <  2:753397/81‑1 (MQ=255)     acacGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAgc                                                                         >  1:753469/1‑81 (MQ=255)     acacGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAg                                                                          <  2:753398/80‑1 (MQ=255)        cGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAgcgc                                                                       >  1:753458/1‑80 (MQ=255)          cTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGat                                                                    <  1:753399/81‑1 (MQ=255)             aaaaaTAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATg                                                                 >  2:753480/1‑81 (MQ=255)                  tAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCt                                                            <  1:753400/81‑1 (MQ=255)                       tcctcGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGTAtg                                                       >  2:753478/1‑81 (MQ=255)                        cctcGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGAtgt                                                      >  2:753501/1‑81 (MQ=255)                         ctcGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGAtgt                                                      <  2:753401/80‑1 (MQ=255)                            gccgccTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGcc                                                   >  1:753498/1‑80 (MQ=255)                            gccgccTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCg                                                  <  1:753402/81‑1 (MQ=255)                              cgccTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCgat                                                <  1:753403/81‑1 (MQ=255)                               gccTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCgatg                                               <  1:753404/81‑1 (MQ=255)                                 cTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCgatgat                                             >  2:753497/1‑81 (MQ=255)                                 cTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCgatgat                                             >  2:753474/1‑81 (MQ=255)                                 cTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCgatgat                                             >  1:753488/1‑81 (MQ=255)                                    aCCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGatgatg                                            <  1:753407/79‑1 (MQ=255)                                    aCCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGatgatg                                            <  2:753406/79‑1 (MQ=255)                                    aCCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGc                                           <  1:753405/80‑1 (MQ=255)                                     ccGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGt                                         >  2:753492/1‑81 (MQ=255)                                     ccGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGt                                         <  1:753408/81‑1 (MQ=255)                                        ttCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTg                                      <  2:753410/81‑1 (MQ=255)                                        ttCAAGCCAACAGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTg                                      <  2:753409/81‑1 (MQ=255)                                         tCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTgg                                     >  1:753483/1‑81 (MQ=255)                                            aGCCAACGGGCGTACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGc                                  <  2:753411/81‑1 (MQ=255)                                              ccAACGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAg                                <  1:753412/81‑1 (MQ=255)                                                  cGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTTCAGGGCg                            <  2:753413/81‑1 (MQ=255)                                                  cGGGCGCACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCg                            <  1:753414/81‑1 (MQ=255)                                                      cgcACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGt                        <  2:753416/81‑1 (MQ=255)                                                      cgcACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGt                        <  1:753415/81‑1 (MQ=255)                                                       gcACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTa                       <  1:753417/81‑1 (MQ=255)                                                       gcACAGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTa                       <  1:753418/81‑1 (MQ=255)                                                        cacaGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAg                      >  2:753496/1‑81 (MQ=255)                                                        cacaGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAg                      <  1:753419/81‑1 (MQ=255)                                                          caGGCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGt                    >  1:753491/1‑81 (MQ=255)                                                            ggCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTc                  <  2:753420/81‑1 (MQ=255)                                                             gCACGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTcc                 <  1:753421/81‑1 (MQ=255)                                                               aCGCGTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCgg               <  2:753422/81‑1 (MQ=255)                                                                cgcgTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGt              <  1:753424/81‑1 (MQ=255)                                                                cgcgTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGt              <  1:753423/81‑1 (MQ=255)                                                                  cgTAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTaa            <  2:753425/81‑1 (MQ=255)                                                                    tAGTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTaaaa          <  1:753426/81‑1 (MQ=255)                                                                      gTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAAc         <  1:753427/80‑1 (MQ=255)                                                                      gTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACg        >  1:753494/1‑81 (MQ=255)                                                                      gTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACg        >  2:753489/1‑81 (MQ=255)                                                                      gTGGCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACg        <  1:753428/81‑1 (MQ=255)                                                                        ggCGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCt      >  1:753502/1‑81 (MQ=255)                                                                         gCGGAGTTCGAGCGCGGTATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCt      <  1:753429/80‑1 (MQ=255)                                                                          cGGAGTTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCTgg    <  1:753430/81‑1 (MQ=255)                                                                            gagtTCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCTGGTc  <  2:753431/81‑1 (MQ=255)                                                                               |                                                                              AATGACACGCTGAAAAATAAACTCCTCGCCGCCTTACCGTTCAAGCCAACGGGCGCACAGGCACGCGTAGTGGCGGAGATCGAGCGCGATATGGCGCTGGATGTGCCGATGATGCGTCTGGTGCAGGGCGATGTAGGTTCCGGTAAAACGCTGGTC  >  NC_000913/3825969‑3826124

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: ‑