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.79.2% 45.5 / 20.8 24coding (975/2574 nt)clpBATP‑dependent chaperone ClpB
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base . (10/9);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCAT  >  NZ_CP009273/2726470‑2726641
                                                                                        |                                                                                   
cTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtt                                                                                    >  1:720220/1‑90 (MQ=255)
          aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                         <  1:237880/90‑1 (MQ=255)
          aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                                         >  2:204283/1‑90 (MQ=255)
                            aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                       >  2:571199/1‑90 (MQ=255)
                            aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                                       >  2:517732/1‑90 (MQ=255)
                                 cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtcg                                                   >  2:552923/1‑90 (MQ=255)
                                 cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                                                  >  2:68021/1‑90 (MQ=255)
                                 cAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGcgtcgt                                                  >  2:543306/1‑90 (MQ=255)
                                  aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                 <  1:67023/90‑1 (MQ=255)
                                  aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                                 <  2:377587/90‑1 (MQ=255)
                                    cAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGc                                               <  1:72946/90‑1 (MQ=255)
                                      gAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAc                                             >  2:197398/1‑90 (MQ=255)
                                       aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                            <  1:277904/90‑1 (MQ=255)
                                       aaGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                            <  1:68021/90‑1 (MQ=255)
                                        aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                           >  2:299894/1‑90 (MQ=255)
                                           cTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACg                                        >  1:377748/1‑90 (MQ=255)
                                             cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGc                                       <  1:519494/90‑1 (MQ=255)
                                             cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCa                                      >  2:766360/1‑90 (MQ=255)
                                               gCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGt                                    <  1:77732/90‑1 (MQ=255)
                                                            tCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac                        >  2:461847/1‑90 (MQ=255)
                                                                   acgacgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgc                >  1:415114/1‑90 (MQ=255)
                                                                      acgTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgcgc              <  2:715115/90‑1 (MQ=255)
                                                                      acgTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCACTGCAATTCACCACGCGCCAgcgc              <  1:552923/90‑1 (MQ=255)
                                                                      acgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGcc             <  2:508406/90‑1 (MQ=255)
                                                                                 cagcaTCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCAt  <  2:377748/90‑1 (MQ=255)
                                                                                        |                                                                                   
CTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCCGGTTTCAGCAT  >  NZ_CP009273/2726470‑2726641

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