Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 977,713 A→G Q463R (CAG→CGG)  mukB → chromosome partitioning protein MukB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913977,7130AG96.8% 103.1 / NA 31Q463R (CAG→CGG) mukBchromosome partitioning protein MukB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (14/16);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTA  >  NC_000913/977576‑977828
                                                                                                                                         |                                                                                                                   
ccAGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAATAGCAGGAAGCGACTGAAAAAAGGCTCTCTCTTGAGCGGa                                                                                                                   >  2:235952/1‑140 (MQ=255)
ccAGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAGGCGACTGAAAAAAAGCTCTCTCTTGAGCGGa                                                                                                                   >  2:206508/1‑140 (MQ=255)
ccAGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCAGGAAGCGACTGAAAAAAAGCTCTCTCTTGAGCGGa                                                                                                                   >  2:282337/1‑140 (MQ=255)
  aGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGaaa                                                                                                                 >  1:638638/1‑140 (MQ=255)
    gCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAAt                                                                                                               >  1:137864/1‑140 (MQ=255)
     cGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATg                                                                                                              >  1:447198/1‑140 (MQ=255)
     cGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATg                                                                                                              >  1:491234/1‑140 (MQ=255)
     cGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCCCTTGAGCGGAAAATg                                                                                                              >  1:254033/1‑140 (MQ=255)
       aTTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAg                                                                                                            >  1:26217/1‑140 (MQ=255)
        ttGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGc                                                                                                           <  1:646773/140‑1 (MQ=255)
        ttGCTGCGCTTAATCGTCCAAATGAACTGTGGCATATGCCGGACTTAATCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGc                                                                                                           <  2:185091/140‑1 (MQ=255)
              cgcTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGGAAAAATGGTCTGTCATGAGCGGAAAATGAGCATGgcg                                                                                                     >  1:182802/1‑140 (MQ=255)
                       cagCCAAAGAACTGTCCCATCTTCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACcgcg                                                                                            <  1:222787/138‑1 (MQ=255)
                        tGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCgcgc                                                                                           >  1:517478/1‑140 (MQ=255)
                          ccAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGATAAAATGCTGTCTGTTGGGCGGAAAAAGAGCATGGCGCAAACCGCGcac                                                                                         >  2:211378/1‑140 (MQ=255)
                                    gtgCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTg                                                                               >  2:450794/1‑140 (MQ=255)
                                     tgCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTTTCTTGAGCAGAAAATGTGCATGGCGCAAACCGCGCACAGCCAGttctg                                                                              >  2:222787/1‑137 (MQ=255)
                                      gCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAg                                                                             <  2:137864/140‑1 (MQ=255)
                                        cATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAGCCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCa                                                                           <  2:447198/140‑1 (MQ=255)
                                         aTCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAg                                                                          <  2:491234/140‑1 (MQ=255)
                                                 gACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCa                                                                  <  1:538930/140‑1 (MQ=255)
                                                                   gccgccGAATGGCTGGACACCTACCAGGCGATAGATCTGCAAGCGACTAAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCaa                                                <  2:182802/140‑1 (MQ=255)
                                                                   gccgccGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCaa                                                <  1:597261/140‑1 (MQ=255)
                                                                   gccgccGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCaa                                                <  2:26217/140‑1 (MQ=255)
                                                                          aaTGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCa                                         <  2:638638/140‑1 (MQ=255)
                                                                          aaTGGCTGGAAACCTGCCAGGCTAATGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCa                                         <  1:282337/140‑1 (MQ=255)
                                                                                                       cTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgc            >  1:554715/1‑140 (MQ=255)
                                                                                                         ggAAGCGACTGAAAACATGCTCTCTCTTGAGCGGAAACTGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGtcgctc          <  1:206508/140‑1 (MQ=255)
                                                                                                              cGACTGAAAAAATGCGCTCTCTTGATCGGCAAATGAGTATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGaa     <  1:211378/140‑1 (MQ=255)
                                                                                                                aCTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACt   >  2:340373/1‑140 (MQ=255)
                                                                                                                 cTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTa  <  1:450794/140‑1 (MQ=255)
                                                                                                                 cTGAAAAAATGCTCTCTCTTGAGCGGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTa  <  2:554715/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                   
CCAGGCGATTGCTGCGCTTAATCGTGCCAAAGAACTGTGCCATCTGCCGGACTTAACCGCCGACTGCGCCGCCGAATGGCTGGAAACCTTCCAGGCGAAAGAGCTGGAAGCGACTGAAAAAATGCTCTCTCTTGAGCAGAAAATGAGCATGGCGCAAACCGCGCACAGCCAGTTTGAGCAGGCTTATCAGCTGGTGGTGGCAATCAACGGCCCACTGGCGCGTAACGAGGCGTGGGATGTCGCTCGCGAACTA  >  NC_000913/977576‑977828

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