Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,951,150 T→G 9.8% intergenic (‑114/+81) PP_3486 ← / ← PP_3487 cytochrome c/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,951,1500TG 9.8% 110.9 / 2.1 39intergenic (‑114/+81)PP_3486/PP_3487cytochrome c/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/10);  new base G (2/2);  total (27/12)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCC  >  NC_002947/3951004‑3951297
                                                                                                                                                   |                                                                                                                                                   
gccctgggcccggggGGAGGTCTATGTTCCTATGCCAAGCCCGCAGTTATGCAATTTCGAACGCGGTGTGGGTGGGGCCTTGCGTCGCTAATGGGCTGCGCAGCAGCCCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCt                                                                                                                                                     <  6:44249/143‑1 (MQ=255)
gTCCAGGGCCTGGGGTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAATCGCAGAAAGAGGGGAGCGCTgc                                                                                                                                                   <  2:331488/149‑1 (MQ=255)
   cagggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcg                                                                                                                                                >  4:130450/1‑149 (MQ=255)
    agggccagggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATGGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTgcgcgc                                                                                                                                               <  8:318762/149‑1 (MQ=255)
            ggcTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGAAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTcg                                                                                                                                       <  6:412846/149‑1 (MQ=255)
                 aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATTTCGCAACCCGTGTGGGTGCGGCCTGGCGTCGCTAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                   <  1:174252/148‑1 (MQ=255)
                 aGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGAc                                                                                                                                   <  2:328330/148‑1 (MQ=255)
                         gTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCg                                                                                                                          <  3:130450/149‑1 (MQ=255)
                                                                       gAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCGCCCCCAAAAGCGGATGCCCCAGAGGCCAAATACCCCGGCGGccc                                                                            >  2:109853/1‑149 (MQ=37)
                                                                             ccTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATa                                                                      >  2:258137/1‑149 (MQ=255)
                                                                                 gCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGGCAATTCCCCAGGCGGCCCATGATACTCa                                                                  >  5:258141/1‑149 (MQ=255)
                                                                                  cGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCac                                                                 >  7:11983/1‑149 (MQ=255)
                                                                                             ggCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAg                                                      >  7:325500/1‑149 (MQ=255)
                                                                                              gCTGCCCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGc                                                     >  7:16884/1‑149 (MQ=255)
                                                                                                tGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg                                                   >  5:134694/1‑149 (MQ=255)
                                                                                                tGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCgg                                                   >  2:125763/1‑149 (MQ=255)
                                                                                                     agcag‑cCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCCCCGACCGGCCAGCGGCTggg                                              >  4:58827/1‑149 (MQ=255)
                                                                                                      gcag‑cCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                             >  7:402154/1‑149 (MQ=255)
                                                                                                      ccag‑cCCCGGCCGGGCCAGCTCAAGCGCAGAACGCGGGCCGCGCGGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGt                                             <  1:109853/148‑1 (MQ=255)
                                                                                                       cag‑cCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTc                                            <  8:11983/149‑1 (MQ=255)
                                                                                                         g‑cCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCaa                                          >  5:326966/1‑149 (MQ=255)
                                                                                                            cccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAAGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                               >  8:289210/1‑142 (MQ=255)
                                                                                                            cccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAAGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTgg                                               <  7:289210/142‑1 (MQ=255)
                                                                                                             ccGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAg                                       >  2:135827/1‑149 (MQ=255)
                                                                                                              cGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGc                                      >  4:175311/1‑149 (MQ=255)
                                                                                                                         gCTCAAGCGCAGAACGCGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACCGTGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCGGGGTCAACAGCTGCATGCGCTg                           <  7:274120/149‑1 (MQ=39)
                                                                                                                         gCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTg                           >  2:363095/1‑149 (MQ=255)
                                                                                                                               gcgcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAAt                     >  4:196282/1‑149 (MQ=255)
                                                                                                                                 gcAGAAAGAGGGGAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcg                   >  8:113718/1‑149 (MQ=255)
                                                                                                                                        gaaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTg                           <  8:403804/130‑1 (MQ=39)
                                                                                                                                        gaaaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTg                           >  7:403804/5‑134 (MQ=39)
                                                                                                                                          aaggAGCGCGCAGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCg          >  8:216567/3‑149 (MQ=255)
                                                                                                                                           gggAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGTCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCgg         >  5:283928/1‑149 (MQ=255)
                                                                                                                                             gAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg       >  3:127238/1‑149 (MQ=255)
                                                                                                                                             gAGCGCTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTg       >  5:78043/1‑149 (MQ=255)
                                                                                                                                               gcgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCa     >  7:384453/1‑149 (MQ=255)
                                                                                                                                                cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccgc               >  6:422322/1‑138 (MQ=255)
                                                                                                                                                cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATcgccgc               <  5:422322/138‑1 (MQ=255)
                                                                                                                                                cgcTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCAt    >  2:336094/1‑149 (MQ=255)
                                                                                                                                                  cTGCGCGCTCCTTTCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGACAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATcc  >  1:362840/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTCCAGGGCCAGGGCTGAGGTCGAAGTTGCCATGCCAAGCACGCAGGTAAGCCATATCGCAAGCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCTGCGCAGCAG‑CCCCGGCAGTGGCAGCTCAAGCGCAGAAAGAGGGGAGCGCTGCGCGCTCTTATCGCGACGCAAGGCCGCTCCCACAAAAGCGGATGCCACAGAGGCCAATTCCCCAGGCGGCCCATGATACTCACCGACCGGCCAGCGGCTGGGTCAACAGCTGCATGCGCTGGGCAATCGCCGCTGGCGGTGCATCC  >  NC_002947/3951004‑3951297

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