Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,038,005	T→C	G34G (GGA→GGG)	recJ ←	ssDNA‑specific exonuclease RecJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,038,005	0	T	C	100.0%	174.2 / NA	50	G34G (GGA→GGG)	recJ	ssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/24);  total (26/24)

CACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCATTTTTTGGTGTT  >  NC_000913/3037872‑3038140                                                                                                                                      |                                                                                                                                        cACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCgg                                                                                                                                   <  2:702014/140‑1 (MQ=255)     gACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTgg                                                                                                                               <  1:749180/140‑1 (MQ=255)        cacaATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGt                                                                                                                            <  2:649016/140‑1 (MQ=255)                tCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACgg                                                                                                                    >  1:661719/1‑140 (MQ=255)                tCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTCCGTTCCAGTTCTTGCGAACTGCGTAGCCCCAGGCTGGCGTATAAACgg                                                                                                                    >  2:326582/1‑140 (MQ=255)                  cgcgTTCCTTCGCGAAAAGCGTTGTAAAGCATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGgcg                                                                                                                  <  1:326582/140‑1 (MQ=255)                     gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACgg                                                                                                                    <  1:161648/135‑1 (MQ=255)                     gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACgg                                                                                                                    >  2:161648/1‑135 (MQ=255)                     gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCgcag                                                                                                               >  1:134106/1‑140 (MQ=255)                     gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCgcag                                                                                                               >  2:571621/1‑140 (MQ=255)                     gTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCgcag                                                                                                               >  1:49281/1‑140 (MQ=255)                       tCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCTGTTCTTGCGCACAGCGGACCCCCCGGCTGGCTTATAAACAGCgcagca                                                                                                             >  2:388097/1‑140 (MQ=255)                           tCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAggg                                                                                                         >  1:695885/1‑140 (MQ=255)                            cgcgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGa                                                                                                        >  2:820623/1‑140 (MQ=255)                             gcgAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAg                                                                                                       <  2:49281/140‑1 (MQ=255)                                        tGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                                                             <  2:503513/139‑1 (MQ=255)                                        tGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGc                                                                                             >  1:503513/1‑139 (MQ=255)                                          tAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCggg                                                                                          <  1:571621/140‑1 (MQ=255)                                               gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGt                                                                                     <  1:593356/140‑1 (MQ=255)                                               gATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGt                                                                                     >  2:714161/1‑140 (MQ=255)                                                    cAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCTAGGCAGACAATTCAGCGGGCAAGTCTGcc                                                                                >  1:398789/1‑140 (MQ=255)                                                          ccTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGCAGGCAATTGAGCGCGCAAGTCTGCCGTTTCa                                                                          >  2:369674/1‑140 (MQ=255)                                                            ttttCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATc                                                                        >  2:1089/1‑140 (MQ=255)                                                                         tCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcga                                                           <  2:661719/140‑1 (MQ=255)                                                                          cAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgac                                                          <  1:820623/140‑1 (MQ=255)                                                                          cAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgac                                                          <  1:115704/140‑1 (MQ=255)                                                                           ggTTGCTGCCAGGGCAGCATACGTTGAACACTGCGTTGCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgacg                                                         <  1:284265/139‑1 (MQ=255)                                                                           aGTTGCTGCCAGGGCCGCATACATTCAACTCTGCGTTCCAGTTCTTGCGCACTGCGTTCCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGcgacg                                                         <  1:305239/140‑1 (MQ=255)                                                                             ttGCTGCCAGGGCAGCATACCTGTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGaa                                                       <  1:388097/140‑1 (MQ=255)                                                                                 tgcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTg                                                   <  2:458028/140‑1 (MQ=255)                                                                                  gcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGt                                                  <  2:453985/140‑1 (MQ=255)                                                                                  gcCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAAATCAGCGGGTAAGTCAGCCGTTTCATAGACTAAACGGCGAAGAAGTTGt                                                  >  2:385253/1‑140 (MQ=255)                                                                                       ggCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTTTGCAGTTACATCGACTCCACGACGACGAAGTTGTATCtc                                             >  2:284265/1‑139 (MQ=255)                                                                                             aTACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGttt                                       >  1:116024/1‑140 (MQ=255)                                                                                               aCCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCa                                     <  1:1089/140‑1 (MQ=255)                                                                                                  tttAACACTGCGTTCCAGTTCTAGCCCACAGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCAcgc                                  <  1:385253/140‑1 (MQ=255)                                                                                                   ttAACACTGCGTTCCTGTTCTTGCGTACTTCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCAcgcg                                 <  1:369674/140‑1 (MQ=255)                                                                                                      acacTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTACACGGCGCCGAAGTTGTATCTTTTGTTTCACGCGAat                              >  1:272818/1‑140 (MQ=255)                                                                                                        acTGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGCGGAAGTTGTATCTGATGTTTCAAGCGAatta                            >  2:398954/1‑140 (MQ=255)                                                                                                          tGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCAcgcg                                 >  1:772718/1‑133 (MQ=255)                                                                                                          tGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCAcgcg                                 <  2:772718/133‑1 (MQ=255)                                                                                                          tGCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCAGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAattatt                          >  1:705122/1‑140 (MQ=255)                                                                                                           gCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTAttt                         >  1:211222/1‑140 (MQ=255)                                                                                                           gCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTAttt                         >  1:746391/1‑140 (MQ=255)                                                                                                           gCGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTAttt                         >  2:848280/1‑140 (MQ=255)                                                                                                            cGTTCCAGTTCTTGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTa                        >  1:77390/1‑140 (MQ=255)                                                                                                               tCCAGTTCTTGCGCAGTGCGTAGCCGCCGTCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCg                     <  1:398954/140‑1 (MQ=255)                                                                                                                        tGCGCACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCAtt            >  1:149051/1‑140 (MQ=255)                                                                                                                           gcACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCAttttt         <  2:24114/140‑1 (MQ=255)                                                                                                                           gcACTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCAttttt         <  1:714161/140‑1 (MQ=255)                                                                                                                              cTGCGTACCCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCATTTTTTgg      <  2:134106/140‑1 (MQ=255)                                                                                                                                  gTACCCCCCGGCTGGCGTATAACCGCCGCGGCAAGGGTGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCATTTTTTGGTGtt  <  1:221643/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        CACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCCTTTTCGACGCCGCTCAGTTGCTGCCAGGGCAGCATACCTTTAACACTGCGTTCCAGTTCTTGCGCACTGCGTACTCCCCGGCTGGCGTATAAACGGCGCAGCAAGGGAGGCAATTCAGCGGGCAAGTCTGCCGTTTCATCGACTTCACGGCGACGAAGTTGTATCTGTTGTTTCACGCGAATTATTTACCGCTGGTCATTTTTTGGTGTT  >  NC_000913/3037872‑3038140

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