Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,287,796 C→T L196L (CTG→CTA ccmF ← cytochrome c‑type biogenesis heme lyase CcmF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,287,7960CT87.0% 67.4 / 1.5 23L196L (CTG→CTAccmFcytochrome c‑type biogenesis heme lyase CcmF
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (12/8);  total (13/10)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCG  >  NZ_CP009273/2287667‑2287920
                                                                                                                                 |                                                                                                                            
ccaGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAaag                                                                                                                       >  2:71948/1‑137 (MQ=255)
  aGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAaagaag                                                                                                                    >  2:190808/1‑138 (MQ=255)
        cccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATgg                                                                                                              >  1:55589/1‑138 (MQ=255)
        cccAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGc                                                                                                             <  1:190808/139‑1 (MQ=255)
                            gCCCATGCGGAACCAAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCgag                                                                                          >  1:235552/1‑138 (MQ=255)
                                   cGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCATCATAGAATCAATGGCAAAAGCGAACGCCACCGAGAAACCCAc                                                                                  <  2:155027/139‑1 (MQ=255)
                                        ccGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGTAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             <  2:340007/139‑1 (MQ=255)
                                        ccGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             <  2:429989/139‑1 (MQ=255)
                                                                         gCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             >  2:362968/1‑106 (MQ=255)
                                                                         gCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             <  1:362968/106‑1 (MQ=255)
                                                                          cTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAACCGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGccc                                           >  2:316455/1‑139 (MQ=255)
                                                                                         gACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAg                            <  2:116329/139‑1 (MQ=255)
                                                                                              gTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGt                       >  2:72533/1‑139 (MQ=255)
                                                                                               tAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGt                       >  1:298608/1‑138 (MQ=255)
                                                                                                 aaaCGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTa                     <  2:55589/138‑1 (MQ=255)
                                                                                                  aaCGCGCATAAGTGCTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAg                   <  2:427221/139‑1 (MQ=255)
                                                                                                                cTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             >  2:30707/1‑67 (MQ=255)
                                                                                                                cTGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTAcc                                                                             <  1:30707/67‑1 (MQ=255)
                                                                                                                 tGTCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAAt    >  2:89549/1‑139 (MQ=255)
                                                                                                                   tCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                 >  2:52981/1‑92 (MQ=255)
                                                                                                                   tCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGAt                                                 <  1:52981/92‑1 (MQ=255)
                                                                                                                   tCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCg  >  2:356668/1‑139 (MQ=255)
                                                                                                                   tCCAGACGCCCGCTTAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCg  >  1:55039/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                            
CCAGCCACCCCAGCCGAGTTCGTAATAGGCCCATGCGGAACCGAGCACGATGCCGAGCGTCAGGAAGATCCACGCTGCCAGCGTCCACGGACGAGTAAAACGCGCATAAGTGCTGTCCAGACGCCCGCTCAGCAAAGAAGCAATGGCAAAAGCAAACGCCACCGAGAAACCCACGTACCCCATATAAAGCAGAGGCGGATGGAAGATCAGCCCCGGATCCTGCAATAGCGGGTTAAGATCGCGACCTTCAATCG  >  NZ_CP009273/2287667‑2287920

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