Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,417,350 T→G K13Q (AAA→CAA)  recE ← Rac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,417,3500TG100.0% 151.1 / NA 44K13Q (AAA→CAA) recERac prophage; exonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (20/24);  total (20/24)

GTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGC  >  NC_000913/1417220‑1417475
                                                                                                                                  |                                                                                                                             
gTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCGAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTGACCGGATGATTTTTGCg                                                                                                                             <  2:484560/133‑1 (MQ=255)
             tAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACAACGTCAGGATCAACCGGTGATTTTTGCGCTTGCCGTAAcgg                                                                                                                >  2:283394/1‑131 (MQ=255)
              aaaaTAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAgg                                                                                                               <  1:137468/133‑1 (MQ=255)
                   aGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAagag                                                                                                          <  1:613588/133‑1 (MQ=255)
                    gctgctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGt                                                                                                         >  1:425982/1‑133 (MQ=255)
                    gctgctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGt                                                                                                         >  2:56775/1‑133 (MQ=255)
                     ctgctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTg                                                                                                        >  1:251635/1‑133 (MQ=255)
                     ctgctCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTGTCAGTAACACGAAGAGTg                                                                                                        >  1:262286/1‑133 (MQ=255)
                          cAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGtttt                                                                                                   <  2:454868/133‑1 (MQ=255)
                             ttttttACATGACTTAACGATCAGGTAGTCTAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTtgtg                                                                                                <  1:283394/133‑1 (MQ=255)
                               ttttACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCt                                                                                              >  1:62877/1‑133 (MQ=255)
                                  tACCTGACTTAACGAGCAGGTAGTCCAGAGAGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCAt                                                                                           <  1:244724/133‑1 (MQ=255)
                                          ttAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGtttttt                                                                                   <  2:251635/133‑1 (MQ=255)
                                          ttAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGtttttt                                                                                   <  1:615904/133‑1 (MQ=255)
                                          ttAACGATCACGTAGGCCAGAGTGACACAGGTGGATTCAAAGTCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGtttttt                                                                                   <  2:262286/133‑1 (MQ=255)
                                                           cAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTTATTGTTTTTTAACCTCAACTCAGATTa                                                                  >  1:272172/1‑133 (MQ=255)
                                                               gTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAAt                                                              >  2:432004/1‑133 (MQ=255)
                                                               gTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAAt                                                              >  1:54446/1‑133 (MQ=255)
                                                               gTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTGTAAGCTAAACTCAAATTAAAAt                                                              >  1:222044/1‑133 (MQ=255)
                                                                   cacaGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTAGTGCTCATCGTTGTTTCACCTCAACTCAGATTAAAATTCGt                                                          >  2:242430/1‑133 (MQ=255)
                                                                        gTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGt                                                     >  2:527547/1‑133 (MQ=255)
                                                                        gTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGt                                                     >  1:344707/1‑133 (MQ=255)
                                                                        gTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGt                                                     >  2:593281/1‑133 (MQ=255)
                                                                         tCGATTCAAAGTCGTCGCTTGCCCACAGTACGACGTCAGGTCCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGtt                                                    <  2:222044/133‑1 (MQ=255)
                                                                         tCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGtt                                                    <  1:527547/133‑1 (MQ=255)
                                                                         tCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGtt                                                    <  2:54446/133‑1 (MQ=255)
                                                                            aTTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTc                                                   >  1:579453/1‑131 (MQ=255)
                                                                            aTTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTc                                                   <  2:579453/131‑1 (MQ=255)
                                                                                  aaTCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATg                                           >  1:487835/1‑133 (MQ=255)
                                                                                       ttgcttgcCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCtt                                      >  1:465785/1‑133 (MQ=255)
                                                                                       ttgcttgcCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCtt                                      >  1:332472/1‑133 (MQ=255)
                                                                                         gcttgcCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGc                                    <  2:447367/133‑1 (MQ=255)
                                                                                         gcttgcCCACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGc                                    <  2:62877/133‑1 (MQ=255)
                                                                                              cccACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTTATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGAt                               <  2:272172/133‑1 (MQ=255)
                                                                                              cccACAGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGAt                               <  2:446726/133‑1 (MQ=255)
                                                                                                cacaGGACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATac                             <  2:425982/133‑1 (MQ=255)
                                                                                                cacaGGACGACGTCAGGTACACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATac                             <  1:242430/133‑1 (MQ=255)
                                                                                                    ggACGACGTCAGGTTTACTGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATacacac                         <  1:208192/133‑1 (MQ=255)
                                                                                                    ggACGACGTCAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATacacac                         >  1:342691/1‑133 (MQ=255)
                                                                                                             cAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATacacac                         <  2:327686/124‑1 (MQ=255)
                                                                                                             cAGGTTCACCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATacacac                         >  1:327686/1‑124 (MQ=255)
                                                                                                                    aCCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTgcgc         <  2:487835/133‑1 (MQ=255)
                                                                                                                    aCCGGATGATTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGACACATACTGTTCATTGCCTgcgc         >  1:223286/1‑133 (MQ=255)
                                                                                                                          tgatTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACg   <  1:593281/133‑1 (MQ=255)
                                                                                                                           gatTTTTGCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGc  <  1:56775/133‑1 (MQ=255)
                                                                                                                                  |                                                                                                                             
GTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTCACCGGATGATTTTTTCGCTTTCCGTAACAGGAAGAGTGGTTTTGTGCTCATTGTTTTTTAACCTCAACTCAGATTAAAATTCGTTTTGTTCAGTGAATGATCTTGCCGGATACACACTGTTCATAGCCTGCGCCATACGC  >  NC_000913/1417220‑1417475

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