Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,947,359 G→A 8.7% G6G (GGC→GGT)  ruvC ← component of RuvABC resolvasome, endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,947,3590GA 8.7% 136.1 / 3.5 46G6G (GGC→GGT) ruvCcomponent of RuvABC resolvasome, endonuclease
Reads supporting (aligned to +/- strand):  ref base G (18/24);  new base A (2/2);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

TCACGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCAT  >  NC_000913/1947220‑1947503
                                                                                                                                           |                                                                                                                                                
tCACGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTAGCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGaat                                                                                                                                          <  4:183514/148‑1 (MQ=255)
  aCGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGaataa                                                                                                                                        <  3:53312/148‑1 (MQ=255)
  aCGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGaata                                                                                                                                         <  4:88112/147‑1 (MQ=255)
   cGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGaataat                                                                                                                                       <  7:190990/148‑1 (MQ=255)
         cATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCAt                                                                                                                                 <  1:48103/148‑1 (MQ=255)
          aTAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCAt                                                                                                                                 >  2:48103/1‑147 (MQ=255)
            aGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATACCGAGAATAATAGCCATCAc                                                                                                                              >  8:76942/1‑148 (MQ=255)
            aGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGAAGACGCCGTAGCCGCTCACGCGCTAACCCCGCTCAATACCGAGAATAAGTGCAATCAc                                                                                                                              >  6:239467/1‑148 (MQ=255)
             gATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCAcg                                                                                                                             >  4:259335/1‑148 (MQ=255)
                       agacgagacgGTAATTCATCCACTGTGTTGCGGATGCATCCGCTACCGAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATACCGAGAATAATAGCCATCACGCGTCTCCGtt                                                                                                                   <  5:239467/148‑1 (MQ=255)
                       agacgagacgGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATACCGAGAATAATAGCCATCACGCGTCTCCGtt                                                                                                                   <  7:76942/148‑1 (MQ=255)
                                 gTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGtt                                                                                                           >  1:373081/1‑146 (MQ=255)
                                  tAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAg                                                                                                        <  3:259335/148‑1 (MQ=255)
                                    aaTCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGGGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGGCTCCGTTTTGCTGTTTAGc                                                                                                       <  4:219040/147‑1 (MQ=255)
                                    aaTCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGc                                                                                                       <  2:398921/147‑1 (MQ=255)
                                             cTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTTCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCtcatc                                                                                             <  5:344882/148‑1 (MQ=255)
                                                  gTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCagag                                                                                         <  8:242440/147‑1 (MQ=255)
                                                   tGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCagagag                                                                                       <  2:373081/148‑1 (MQ=255)
                                                        aTGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCt                                                                                  >  4:123867/1‑148 (MQ=255)
                                                                          tAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAg                                                                                                        <  5:108124/108‑1 (MQ=255)
                                                                          tAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCAGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAg                                                                                                        >  6:108124/1‑108 (MQ=255)
                                                                             gACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCa                                                             >  6:205392/1‑148 (MQ=255)
                                                                              aCAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAgaga                                                                  >  6:312862/1‑142 (MQ=255)
                                                                              aCAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAgaga                                                                  <  5:312862/142‑1 (MQ=255)
                                                                              aCAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCAc                                                            >  6:354652/1‑148 (MQ=255)
                                                                                  ttGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGt                                                         >  4:191257/1‑147 (MQ=255)
                                                                                        acctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGt                                                   >  7:88064/1‑147 (MQ=255)
                                                                                         cctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTaa                                                 <  8:88064/148‑1 (MQ=255)
                                                                                         cctacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTaa                                                 <  3:123867/148‑1 (MQ=255)
                                                                                           tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTaa                                                 >  3:52552/1‑146 (MQ=255)
                                                                                           tacctGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTaa                                                 <  4:52552/146‑1 (MQ=255)
                                                                                               tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCtc                                                                          >  2:78527/1‑117 (MQ=255)
                                                                                               tGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCtc                                                                          <  1:78527/117‑1 (MQ=255)
                                                                                                ggCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACtt                                             <  7:134003/145‑1 (MQ=255)
                                                                                                ggCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACtt                                             >  8:134003/1‑145 (MQ=255)
                                                                                                ggCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCt                                          >  2:40888/1‑148 (MQ=255)
                                                                                                                gCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAAt                           <  1:40888/147‑1 (MQ=255)
                                                                                                                   ggTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                       <  5:205392/148‑1 (MQ=255)
                                                                                                                     tCACGCGCGAACCCGGATCAATGCCGTGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                       <  6:297926/146‑1 (MQ=255)
                                                                                                                     tCACGCGCGAACCCGGATCAATGCCGGGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                       >  5:297926/1‑146 (MQ=255)
                                                                                                                     tCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                       <  6:226631/146‑1 (MQ=255)
                                                                                                                     tCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTc                       >  5:226631/1‑146 (MQ=255)
                                                                                                                              aaCCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGa            <  3:191257/148‑1 (MQ=255)
                                                                                                                                     aTCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACg     >  8:29879/1‑148 (MQ=255)
                                                                                                                                        aaTGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCa   <  5:184919/147‑1 (MQ=255)
                                                                                                                                         aTGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCAt  >  7:14421/1‑147 (MQ=255)
                                                                                                                                           |                                                                                                                                                
TCACGCCCGCATAGATGAGTTTCAGACGAGACGGTAAATCATCCACTTTGGTGCGGATGCATCCGCTACCCAGGTAGGACAGTTGCCTACCTACCTGGCGGATGACGCCGTAGCCGGTCACGCGCGAACCCGGATCAATGCCGAGAATAATAGCCATCACGCGTCTCCGTTTTGCTGTTTAGCAGGCCTCATCAGAGAGTCGCTGCGACCTCATCAGAGATTTCACCGTTATGGTAAACTTCCTGCACGTCGTCGCAATCTTCCAGCATATCGATCAGACGCAT  >  NC_000913/1947220‑1947503

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑