Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,726,562 (T)5→4 coding (971/2574 nt) clpB ← ATP‑dependent chaperone ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,726,5580T.88.9% 50.3 / 4.8 18coding (975/2574 nt)clpBATP‑dependent chaperone ClpB
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (9/7);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA  >  NZ_CP009273/2726470‑2726640
                                                                                        |                                                                                  
cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                                                   <  2:198266/90‑1 (MQ=255)
   tCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttc                                                                                <  1:308475/90‑2 (MQ=255)
       gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATg                                                                           <  2:225916/90‑1 (MQ=255)
        ccACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGt                                                                          <  1:81688/90‑1 (MQ=255)
          aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                        >  2:526581/1‑90 (MQ=255)
            gCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTg                                                                      <  2:192247/90‑1 (MQ=255)
                         ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc                                                          >  2:541406/1‑90 (MQ=255)
                         ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                         >  2:50238/1‑90 (MQ=255)
                         ggTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCa                                                         >  2:256582/1‑90 (MQ=255)
                                  aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                <  1:1227/90‑1 (MQ=255)
                                      gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                >  2:551556/1‑86 (MQ=255)
                                      gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                <  1:551556/86‑1 (MQ=255)
                                        aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                           <  2:29325/90‑1 (MQ=255)
                                        aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                          >  1:272977/1‑90 (MQ=255)
                                         ggCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTa                                         <  1:489812/90‑1 (MQ=255)
                                           cTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACg                                       >  2:464277/1‑90 (MQ=255)
                                                 aacaaacaCTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGc                                 <  2:211609/90‑1 (MQ=255)
                                                                   acgacgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgc               >  2:210963/1‑90 (MQ=255)
                                                                         ttCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGt         >  2:193552/1‑90 (MQ=255)
                                                                                gcagcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCa  >  2:336422/1‑90 (MQ=255)
                                                                                        |                                                                                  
CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCA  >  NZ_CP009273/2726470‑2726640

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