Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,859,998 T→C I305V (ATT→GTT)  rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,859,9980TC92.6% 71.6 / ‑0.6 27I305V (ATT→GTT) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (11/14);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC  >  NZ_CP009273/2859915‑2860084
                                                                                   |                                                                                      
tacttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCg                                                                                  >  1:212053/1‑90 (MQ=255)
  cttactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGaa                                                                                <  2:56808/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                              <  1:409387/90‑1 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                              >  1:474565/1‑90 (MQ=255)
    tactCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                              <  2:575539/90‑1 (MQ=255)
     actCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAacac                                                                             >  1:33154/1‑90 (MQ=255)
                      tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                              >  2:524975/1‑72 (MQ=255)
                      tCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAaca                                                                              <  1:524975/72‑1 (MQ=255)
                           atTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTc                                                       <  1:199751/90‑1 (MQ=255)
                             tCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCAc                                                     <  2:214323/90‑1 (MQ=255)
                                 ccccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                                                 <  2:8104/90‑1 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                                                 <  1:566126/88‑1 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGAcc                                                 >  2:566126/1‑88 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa                                               >  2:505805/1‑90 (MQ=255)
                                   ccTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTa                                               >  1:515520/1‑90 (MQ=255)
                                     tGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACa                                             <  2:224400/90‑1 (MQ=255)
                                      gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACAt                                            <  2:234272/90‑1 (MQ=255)
                                      gCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACAt                                            <  2:360558/90‑1 (MQ=255)
                                          ttGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTc                                        >  2:90861/1‑90 (MQ=255)
                                                   ttCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc                               <  1:195925/90‑1 (MQ=255)
                                                   ttCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGc                               <  1:426332/90‑1 (MQ=255)
                                                    tCGCGCAAACGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGcc                              >  2:471979/1‑90 (MQ=255)
                                                          aaaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc                            <  2:147312/86‑1 (MQ=255)
                                                          aaaCGGCGCAGGCCTTCAACCTGAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGc                            >  1:147312/1‑86 (MQ=255)
                                                                                gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc  <  2:289705/90‑1 (MQ=255)
                                                                                gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc  <  2:212053/90‑1 (MQ=255)
                                                                                gAACCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAc  >  1:510471/1‑90 (MQ=255)
                                                                                   |                                                                                      
TACTTACTCGCGGAACAGCGCTTCGATATTCAGCCCCTGCGTTTGCAGGATTTCGCGCAAACGGCGCAGGCCTTCAACCTGAATCTGGCGAACACGTTCACGGGTGAGGCCAATTTCACGACCTACATCTTCCAGTGTTGCCGCTTCGTACCCCAGCAAACCGAATCGAC  >  NZ_CP009273/2859915‑2860084

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