Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,777,524 C→G 8.6% intergenic (+258/+34) PP_5476 → / ← PP_1585 hypothetical protein/antidote protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,777,5240CG 8.6% 176.9 / 5.8 58intergenic (+258/+34)PP_5476/PP_1585hypothetical protein/antidote protein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base G (2/3);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAG  >  NC_002947/1777386‑1777671
                                                                                                                                          |                                                                                                                                                   
gTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGcc                                                                                                                                           >  4:253533/1‑149 (MQ=255)
 tttAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGccc                                                                                                                                          <  6:44082/149‑1 (MQ=255)
     aCTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGGGCGTGGCGccgcga                                                                                                                                       <  5:236076/148‑5 (MQ=255)
       ttGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGtt                                                                                                                                     >  2:211309/1‑148 (MQ=255)
          gCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTAtt                                                                                                                                 <  3:253533/149‑1 (MQ=255)
               ccAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATc                                                                                                                             <  2:82962/148‑1 (MQ=255)
                cAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAg                                                                                                                           >  2:164652/1‑149 (MQ=255)
                         aTGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTc                                                                                                                  <  2:334637/149‑1 (MQ=255)
                          tGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTc                                                                                                                  >  6:221196/1‑148 (MQ=255)
                            ccGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCaa                                                                                                                <  5:221196/148‑1 (MQ=255)
                                 gagCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGAc                                                                                                           <  4:307423/148‑1 (MQ=255)
                                        tAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCt                                                                                                   <  4:138696/149‑1 (MQ=255)
                                          ggTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTgg                                                                                                 <  3:92143/149‑1 (MQ=255)
                                           gTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                                >  4:286171/1‑149 (MQ=255)
                                                 gtgtACATTGTGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgac                                                                                           <  4:239132/148‑1 (MQ=255)
                                                  tgtACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGAccc                                                                                         >  6:186267/1‑149 (MQ=255)
                                                    tACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCga                                                                                       <  1:164652/149‑1 (MQ=255)
                                                      cATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCCCGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCgat                                                                                      >  5:235286/1‑148 (MQ=255)
                                                       aTTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGAGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCgatga                                                                                    >  1:256275/1‑149 (MQ=255)
                                                            ggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgacc                                                                                          >  8:162529/1‑138 (MQ=255)
                                                            ggTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgacc                                                                                          <  7:162529/138‑1 (MQ=255)
                                                              ttGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATc                                                                             <  5:186267/149‑1 (MQ=255)
                                                                        ccccAATCTCTCCGTTTTTCGATGTCTGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGc                                                                   >  1:186287/1‑149 (MQ=255)
                                                                         cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAtgt                                                                                                                      <  8:183092/97‑1 (MQ=255)
                                                                         cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGAtgt                                                                                                                      >  7:183092/1‑97 (MQ=255)
                                                                         cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCAt                                                                              <  3:296498/137‑1 (MQ=255)
                                                                         cccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCAt                                                                              >  4:296498/1‑137 (MQ=255)
                                                                          ccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgacc                                                                                          <  6:134759/124‑1 (MQ=255)
                                                                          ccAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGgaccgacc                                                                                          >  5:134759/1‑124 (MQ=255)
                                                                             aTCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTc                                                              <  3:286171/149‑1 (MQ=255)
                                                                                  tcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGc                                                                                                    <  7:235725/106‑1 (MQ=255)
                                                                                  tcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGc                                                                                                    >  8:235725/1‑106 (MQ=255)
                                                                                  tcCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTAGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCtg                                                         >  8:19959/1‑149 (MQ=255)
                                                                                   ccGTTTTTCGATGTCAGTCGCAGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGc                                                                   <  2:181223/138‑1 (MQ=255)
                                                                                   ccGTTTTTCGATGTCAGTCGCAGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAACCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGc                                                                   >  1:181223/1‑138 (MQ=255)
                                                                                     gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGttt                                                                                                                                    >  4:78171/1‑71 (MQ=255)
                                                                                     gTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGttt                                                                                                                                    <  3:78171/71‑1 (MQ=255)
                                                                                                tCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                                            <  1:211309/148‑1 (MQ=255)
                                                                                                tCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTAGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGaa                                            <  7:19959/148‑1 (MQ=255)
                                                                                                    tCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAg                                           >  4:109102/1‑145 (MQ=255)
                                                                                                    tCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAg                                           <  3:109102/145‑1 (MQ=255)
                                                                                                       cgGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAg                                           >  5:358285/1‑142 (MQ=255)
                                                                                                       cgGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAg                                           <  6:358285/142‑1 (MQ=255)
                                                                                                          gATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTg                                   >  8:272763/1‑147 (MQ=255)
                                                                                                               tGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGtt                             >  3:284101/1‑148 (MQ=255)
                                                                                                                 cGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGa                           <  2:186287/148‑1 (MQ=255)
                                                                                                                     aaCGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCa                                                   <  4:202185/120‑1 (MQ=255)
                                                                                                                     aaCGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCa                                                   >  3:202185/1‑120 (MQ=255)
                                                                                                                        gATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGaa                    >  6:142262/1‑148 (MQ=255)
                                                                                                                            cTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAg               <  6:235286/149‑1 (MQ=255)
                                                                                                                              aggggcgCCACGGGCGTGGCGCCGCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCa                                                   >  6:193903/4‑111 (MQ=255)
                                                                                                                              aggggcgCCACGGGCGTGGCGCCGCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCa                                                   <  5:193903/108‑1 (MQ=255)
                                                                                                                               ggggcgCCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                                >  6:337173/3‑65 (MQ=38)
                                                                                                                               ggggcgCCACGGGCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGga                                                                                                <  5:337173/63‑1 (MQ=38)
                                                                                                                               gcggcgCCACGCCCGTGGCGCCCCGGTTTATTGATCAGGATGTGTTCAAGCGCCCCGACGCTGGCCCGACCCGATGATCTCCTCGGCGTATTGGTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCATAACTCAgggg            >  3:39200/1‑149 (MQ=255)
                                                                                                                                 ggcgCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGggtg          <  2:256275/149‑1 (MQ=255)
                                                                                                                                     ccACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGgtgt      >  1:126882/1‑149 (MQ=255)
                                                                                                                                         gCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAg  <  5:142262/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GTTTAACTTGGCGTCCCAGGCTTTGATGCCGTAGAGCGTGTAGGTCATTGTGTACATTGGGGTTGCAAGTAGCCCCAATCTCTCCGTTTTTCGATGTCAGTCGCGGGATTATGCGGGAACGATCCTCGCGGCGCCACGCCCGTGGCGCCCCTGTTTATTGATCAGGATGTGTTCAAGCGACCAGACGCTGGACCGACCCGATGATCTCATCGCCGTGTTGCTGTTCGGCTGTGCGCAGATCGAAGGCCGTCTGCAGGTTGAGCCAGAACTCAGGGGTGGTGTCCAG  >  NC_002947/1777386‑1777671

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: ‑