Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,826,614	A→G	T469A (ACC→GCC)	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,614	0	A	G	100.0%	159.5 / NA	45	T469A (ACC→GCC)	recG	ATP‑dependent DNA helicase RecG
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/17);  total (28/17)

TGACCGCCACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGG  >  NC_000913/3826486‑3826749                                                                                                                                 |                                                                                                                                        tGACCGCCACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCAtt                                                                                                                              >  2:480854/1‑140 (MQ=255)    ccGCCACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGAc                                                                                                                           >  1:718454/1‑140 (MQ=255)    ccGCCACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGAc                                                                                                                           >  1:729833/1‑140 (MQ=255)       ccACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACcgc                                                                                                                        >  1:317903/1‑140 (MQ=255)       ccACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACcgc                                                                                                                        >  1:915617/1‑140 (MQ=255)        cACGCCGACCCCCCGCCCGCTGGCAATGACTGCGTATGCCGATCTCGATCCCTGGGTGATAGATGAGCGGCCGCCGGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATCCCCGCCGTGCCGACATCATTGACcgcg                                                                                                                       <  2:729833/140‑1 (MQ=255)        cACCCCGACCCCCCGCCGGCTGGAAATGCCTGCGTATGCCGATCTCGATCCCGGGGTGATAGAGGAGCTCCCGCCGGGCCGCACGCCAGTGACTACGGTCGCTATCCCTGATACCCGCCGTGCCGACATCATTGACcgcg                                                                                                                       <  2:422577/140‑1 (MQ=255)               aTCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGcacca                                                                                                                >  2:994594/1‑140 (MQ=255)                ccccccGCACGCTGCCAATCGCCGCGAAGGCGGATCTCGATACCGCGTTGAAGGATGAGCCCCCCCAAGGCCGCACGCCAGTGACTAGGTTCGTTTTCCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGcaccac                                                                                                               <  2:317903/139‑1 (MQ=255)                  ccccGCCCGCGGGCAATGAGTGCGTATGCCGATCCGGATCCCTCGGTGATAGATGGCCTGCCGCCGGGCCGCACGCGAGTGACTACGGTCGCTATCCCTGATCCCCGCCGTGCCGACATCATTGACCGCGTGCACCACGc                                                                                                             <  2:51428/140‑1 (MQ=255)                      gCACGCTGGCAATGCCTGCGTATGCCGATCTCGATCCCTGGGTGATAGATGAGCGGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACACCATTGACCGCGTGCACCACGCCTGc                                                                                                         <  2:281166/140‑1 (MQ=255)                       cACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTTGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCa                                                                                                        >  1:633021/1‑140 (MQ=255)                            tGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACt                                                                                                   >  1:322928/1‑140 (MQ=255)                            tGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACt                                                                                                   >  2:344334/1‑140 (MQ=255)                                      tGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGgtcgtc                                                                                         >  1:963452/1‑140 (MQ=255)                                       gCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCa                                                                                        >  1:145676/1‑140 (MQ=255)                                        ggTATGCCGTTCTCGATCCCTGGGTGATAGATGGGCTCCCGCCGGGCCGCACGCCAGTGACTAGGTTTGCTATTCCTGATCCCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAg                                                                                       <  2:633021/139‑1 (MQ=255)                                        cGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGgtctgtct                                                                                       >  1:179225/1‑135 (MQ=255)                                          tATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGc                                                                                     >  1:349237/1‑140 (MQ=255)                                             gCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATa                                                                                  <  1:301527/140‑1 (MQ=255)                                               cGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACt                                                                                <  2:718454/140‑1 (MQ=255)                                                  tctcGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGgtc                                                                                            >  2:543699/1‑125 (MQ=255)                                                  tctcGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGgtc                                                                                            <  1:543699/125‑1 (MQ=255)                                                      gATCCCTCGGTGATAGATGAGCTGCCCCCGGGCCGCACGCCAGTGACTCCGGTCGTTTTTCCTGATCCCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTg                                                                         <  2:349237/140‑1 (MQ=255)                                                           cTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGt                                                                    <  2:963452/140‑1 (MQ=255)                                                             cGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGttg                                                                  >  2:377667/1‑140 (MQ=255)                                                                tGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGttgatt                                                               >  1:400978/1‑140 (MQ=255)                                                                  aTAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGttgattga                                                             <  1:377667/140‑1 (MQ=255)                                                                    aGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAag                                                           >  2:467723/1‑140 (MQ=255)                                                                               cgccAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTg                                                >  2:258297/1‑140 (MQ=255)                                                                               cgccAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTg                                                >  2:382852/1‑140 (MQ=255)                                                                                        gCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAgcg                                       <  1:382852/140‑1 (MQ=255)                                                                                        gCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAgcg                                       <  1:258297/140‑1 (MQ=255)                                                                                         cACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAgcgc                                      >  1:944795/1‑140 (MQ=255)                                                                                                   aCTACGGTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGa                            >  2:170195/1‑140 (MQ=255)                                                                                                        ggTCGCTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTa                       >  2:74226/1‑140 (MQ=255)                                                                                                             cTATTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTgg                  >  2:369721/1‑140 (MQ=255)                                                                                                               aTTCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGa                <  1:74226/140‑1 (MQ=255)                                                                                                                ttCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGaa               <  2:145676/140‑1 (MQ=255)                                                                                                                ttCCTGATACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGaa               <  1:83006/140‑1 (MQ=255)                                                                                                                      aTACCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTg         >  1:105257/1‑140 (MQ=255)                                                                                                                        aCCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGaa       >  2:870342/1‑140 (MQ=255)                                                                                                                        aCCCGCCGTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGaa       >  1:725352/1‑140 (MQ=255)                                                                                                                          ccgccgTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAAc     >  2:394417/1‑140 (MQ=255)                                                                                                                             ccgTGCCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTgg  >  1:456235/1‑140 (MQ=255)                                                                                                                                 |                                                                                                                                        TGACCGCCACGCCGATCCCCCGCACGCTGGCAATGACTGCGTATGCCGATCTCGATACCTCGGTGATAGATGAGCTGCCGCCAGGCCGCACGCCAGTGACTACGGTCGCTATTCCTGATACCCGCCGTACCGACATCATTGACCGCGTGCACCACGCCTGCATAACTGAAGGTCGTCAGGCATACTGGGTTTGTACGTTGATTGAAGAGTCGGAATTGCTGGAAGCGCAGGCGGCGGAAGCTACCTGGGAAGAGTTGAAACTGG  >  NC_000913/3826486‑3826749

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