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.95.2% 75.5 / ‑1.6 21coding (975/2574 nt)clpBATP‑dependent chaperone ClpB
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (13/7);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

TTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCC  >  NZ_CP009273/2726472‑2726630
                                                                                      |                                                                        
ttCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCtttt                                                                       >  2:573705/1‑90 (MQ=255)
  cAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttttca                                                                     <  2:201034/90‑3 (MQ=255)
     gCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATg                                                                 >  2:107795/1‑90 (MQ=255)
      ccACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGt                                                                <  2:327173/90‑1 (MQ=255)
        aCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTAc                                                              >  2:297121/1‑90 (MQ=255)
                tCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGata                                                      >  2:489459/1‑90 (MQ=255)
                  gCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATAtt                                                    >  1:488621/1‑90 (MQ=255)
                          aTCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGc                                            <  1:196004/90‑1 (MQ=255)
                                aaCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                      <  2:104207/90‑1 (MQ=255)
                                aaCAGAAGGCTCGGCAACAAACACTTCCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTg                                      <  1:182521/90‑1 (MQ=255)
                                 aCAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTgg                                     >  1:255563/1‑90 (MQ=255)
                                      aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCAcc                                 >  1:506133/1‑90 (MQ=255)
                                      aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                >  1:441093/1‑90 (MQ=255)
                                      aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                >  2:299134/1‑90 (MQ=255)
                                      aGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCt                                >  2:284949/1‑90 (MQ=255)
                                          tCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGCTATTCGTCAAGCGTCGTGGCACCTACGc                            <  2:284548/90‑1 (MQ=255)
                                           cGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCa                           >  1:223384/1‑90 (MQ=255)
                                                         ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac             <  1:365045/90‑1 (MQ=255)
                                                         ttCTGGAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTcaccac             >  1:267543/1‑90 (MQ=255)
                                                              gAAACGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGcc        >  2:571399/1‑90 (MQ=255)
                                                                aaCGACGTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAg      <  2:255563/90‑1 (MQ=255)
                                                                 acgacgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAgc     >  1:76200/1‑90 (MQ=255)
                                                                    acgTTCCAGCGCAGCATCTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGcc  >  2:164692/1‑90 (MQ=255)
                                                                                      |                                                                        
TTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGGCACCTACGCAGTGCAATTCACCACGCGCCAGCGCC  >  NZ_CP009273/2726472‑2726630

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