Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,938,718 T→A 11.5% L408Q (CTG→CAG)  rbbA → ribosome‑associated ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,938,7180TA11.5% 92.9 / 2.3 34L408Q (CTG→CAG) rbbAribosome‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base T (14/16);  new base A (2/2);  total (16/18)
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

CTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5938579‑5938864
                                                                                                                                           |                                                                                                                                                  
cttgcgcTGCAGGCCCGCCTTTTCGACTTGCCCACAGCCGAGAGCCCCTATCGCATCGATGAGCTGATCGACCACTTCTACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCa                                                                                                                                          <  1:451005/144‑1 (MQ=255)
    cccTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTgg                                                                                                                                      >  2:604931/1‑150 (MQ=255)
     ccTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGc                                                                                                                                     <  1:645003/150‑1 (MQ=255)
        gCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTgg                                                                                                                                      >  2:611042/1‑146 (MQ=255)
        gCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTgg                                                                                                                                      <  1:611042/146‑1 (MQ=255)
           tGCCCGCCTGTTCGACATGCAACACGCCACAACCACCCATCCCATTTATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCt                                                                                                                               <  2:606261/150‑1 (MQ=255)
                             gCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTg                                                                                                                                             <  2:501389/118‑1 (MQ=255)
                             gCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTg                                                                                                                                             >  1:501389/1‑118 (MQ=255)
                                cAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGAt                                                                                                          >  2:608540/1‑150 (MQ=255)
                                 aaaGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATc                                                                                                         >  1:329431/1‑150 (MQ=255)
                                      cGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTcga                                                                                                    <  1:608540/150‑1 (MQ=255)
                                                      aTCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGc                                                                                    >  1:558870/1‑150 (MQ=255)
                                                       tCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCg                                                                                   >  1:570183/1‑150 (MQ=255)
                                                              gccgATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTcc                                                                                       <  1:245526/136‑1 (MQ=255)
                                                              gccgATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTcc                                                                                       >  2:245526/4‑139 (MQ=255)
                                                              gCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGACAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGAccc                                                                            >  2:283788/1‑150 (MQ=255)
                                                                   tCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCccgccg                                                                       >  1:103391/1‑150 (MQ=255)
                                                                       gcgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGc                                                                                                                                           >  1:553490/1‑78 (MQ=255)
                                                                       gcgcTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGACAGGCGCTGGc                                                                                                                                           <  2:553490/78‑1 (MQ=255)
                                                                           ttCGACCTCGCCGCTATCGCCGACTAACCGTCTGGCGCCATGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGACCGGAAGTCTTAATCCTCGACGAAACGGCCTGCGGCGTCGACTGCGCTGCACGcgag                                                               >  1:39682/1‑149 (MQ=255)
                                                                              gACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCgacgac                                                            <  2:329431/150‑1 (MQ=255)
                                                                                     cctcTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGgc                                                     <  1:474420/147‑1 (MQ=255)
                                                                                                   cAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGaa                                       <  1:604931/150‑1 (MQ=255)
                                                                                                    aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCg                                                                                   >  1:570125/1‑105 (MQ=255)
                                                                                                    aaCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCg                                                                                   <  2:570125/105‑1 (MQ=255)
                                                                                                     aCCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCTCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACt                                     <  1:612872/150‑1 (MQ=255)
                                                                                                      ccGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTg                                    >  1:738657/1‑150 (MQ=255)
                                                                                                             gcgcCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCTGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCcgcg                             >  2:359334/1‑150 (MQ=255)
                                                                                                              cgcCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGa                            <  1:456152/150‑1 (MQ=255)
                                                                                                                         ttGCTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCCGTTCTGCATCGCCCCGAAGTCTTAATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTTGCGCCTGCTTGTAGAACTGTCCCGCGAGCAGGTCGTGa                 <  2:39682/150‑1 (MQ=255)
                                                                                                                             ccactgcCAGCGACAGGCGCTGGCGGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGcgacga                                                             >  1:275909/6‑102 (MQ=255)
                                                                                                                                ctgcCAGCGACAGGCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTg                  <  2:443309/140‑1 (MQ=255)
                                                                                                                                     aGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCctct     <  2:566222/150‑1 (MQ=255)
                                                                                                                                        gcCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAATTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTcca  <  1:283788/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CTGGCCCTGCATGCCCGCCTGTTCGACTTGCCCAAAGCCGAGAGCGCCCAGCGCATCGATGAGCTGATCGAGCGCTTCGACCTCGCCGCTATCGCCGACCAACCGTCTGGCGCCCTGCCCCTTGGTTTGCGCCAGCGCCTGTCGCTGGCAGTGGCGGTTCTGCATCGCCCGGAAGTCTTGATCCTCGACGAACCGACCTCCGGCGTCGACCCCGCCGCCCGCGACGACTTCTGGCGCCTGCTGGTAGAACTGTCCCGCGAGCAGGGCGTGACCATTTTCCTCTCCA  >  NC_002947/5938579‑5938864

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