Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,370,587 G→T R65S (CGT→AGT)  sspA ← stringent starvation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,370,5870GT94.7% 58.3 / NA 19R65S (CGT→AGT) sspAstringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (6/12);  minor base A (1/0);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGA  >  NZ_CP009273/3370505‑3370673
                                                                                  |                                                                                      
tAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATccacc                                                                                 <  1:80843/90‑1 (MQ=255)
 aaaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATccacca                                                                                <  2:437419/90‑1 (MQ=255)
 aaaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATccacca                                                                                <  2:296578/90‑1 (MQ=255)
        cATCAGTGGCGGATGCGGGAAAAGCACATCCAGAAATTCCATAATGATGCGAGATTCCCACCGGGTCCGCTCACACTTCACCAGGTTCgg                                                                         >  2:332543/1‑90 (MQ=255)
          tCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGaa                                                                       <  1:543436/90‑1 (MQ=255)
                gCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGt                                                                            >  2:482828/1‑79 (MQ=255)
                gCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGt                                                                            <  1:482828/79‑1 (MQ=255)
                       cGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGctc                                                                  <  2:554998/82‑1 (MQ=255)
                       cGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGctc                                                                  >  1:554998/1‑82 (MQ=255)
                                cTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTc                                                 <  1:584035/90‑1 (MQ=255)
                                            ttCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                     >  1:589048/1‑90 (MQ=255)
                                                aTAATGATGCGAGATTCCCACAGGGTCAGCTCACTATccacca                                                                                >  1:207200/1‑43 (MQ=255)
                                                aTAATGATGCGAGATTCCCACAGGGTCAGCTCACTATccacca                                                                                <  2:207200/43‑1 (MQ=255)
                                                     gatgCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                            <  1:332543/90‑1 (MQ=255)
                                                        gCGAGATTCCCACAGGGTCAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTc                         <  1:164392/90‑1 (MQ=255)
                                                                      gggTCAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGtt           >  1:66728/1‑90 (MQ=255)
                                                                          cAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAt                             >  2:425588/1‑68 (MQ=255)
                                                                          cAGCTCACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAt                             <  1:425588/68‑1 (MQ=255)
                                                                               cACTATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGa  <  2:204702/90‑1 (MQ=255)
                                                                                  |                                                                                      
TAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGA  >  NZ_CP009273/3370505‑3370673

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