Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,370,599 T→G T61P (ACC→CCC)  sspA ← stringent starvation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,370,5990TG100.0% 120.0 / NA 32T61P (ACC→CCC) sspAstringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (17/15);  total (17/15)

CAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACC  >  NZ_CP009273/3370516‑3370683
                                                                                     |                                                                                    
cAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAAc                                                                                <  2:282867/90‑1 (MQ=255)
    ggCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGctc                                                                            <  1:4635/90‑1 (MQ=255)
    ggCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGctc                                                                            <  1:711059/90‑1 (MQ=255)
     gCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGctct                                                                           >  2:518679/1‑90 (MQ=255)
     gCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGGCAGCTCACGATC‑‑CACCAGGGGCGGAACGctct                                                                           >  1:363901/1‑90 (MQ=255)
                   cGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAgg                                                             <  1:852163/90‑1 (MQ=255)
                       cATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTc                                                           <  1:57992/90‑1 (MQ=255)
                           cAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCAGAACGCTCTGATTCGGGTTGAGGTCAATCAg                                                     <  1:165024/90‑1 (MQ=255)
                                   ccATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAg                                             >  1:293863/1‑90 (MQ=255)
                                   ccATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAg                                             >  2:492565/1‑90 (MQ=255)
                                   ccATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAg                                             >  2:403774/1‑90 (MQ=255)
                                   ccATAATGATGCGAGATTCCCACAGGGTCAGATCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAg                                             >  2:436390/1‑90 (MQ=255)
                                    cATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAgg                                            <  1:518679/90‑1 (MQ=255)
                                    cATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAgg                                            <  1:97104/90‑1 (MQ=255)
                                       aaTGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCgg                                         >  2:648892/1‑90 (MQ=255)
                                        atgatgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGa                                        >  2:703505/1‑90 (MQ=255)
                                          gatgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                      >  1:275721/1‑90 (MQ=255)
                                           atgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGctc                                                                            <  1:459033/51‑1 (MQ=255)
                                           atgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGctc                                                                            >  2:459033/1‑51 (MQ=255)
                                           atgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTg                                     >  1:382113/1‑90 (MQ=255)
                                           atgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTg                                     >  2:525637/1‑90 (MQ=255)
                                            tgCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGt                                    >  2:235942/1‑90 (MQ=255)
                                                            ggTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCa                           <  1:509251/83‑1 (MQ=255)
                                                            ggTCAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCa                           >  2:509251/1‑83 (MQ=255)
                                                               cAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGAt                 >  2:361559/1‑90 (MQ=255)
                                                               cAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGAt                 >  1:297196/1‑90 (MQ=255)
                                                               cAGCTCACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGAt                 >  2:574801/1‑90 (MQ=255)
                                                                    cACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGa            <  1:648892/90‑1 (MQ=255)
                                                                    cACGATC‑‑CACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGa            <  1:266661/90‑1 (MQ=255)
                                                                                caGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACAcc  <  1:703505/90‑1 (MQ=255)
                                                                                caGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACAcc  <  2:206605/90‑1 (MQ=255)
                                                                                caGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACAcc  <  1:228197/90‑1 (MQ=255)
                                                                                     |                                                                                    
CAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATC‑‑CACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACC  >  NZ_CP009273/3370516‑3370683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: