Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 798,663 C→G 7.3% intergenic (+18/+9) fklB‑I → / ← PP_0685 FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947798,6630CG 7.3% 158.8 / 2.3 53intergenic (+18/+9)fklB‑I/PP_0685FKBP‑type peptidyl‑prolyl cis‑trans isomerase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base G (2/2);  total (32/21)
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.33e-01

GAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACG  >  NC_002947/798522‑798803
                                                                                                                                             |                                                                                                                                            
gaACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGc                                                                                                                                       >  6:265083/1‑149 (MQ=255)
   cGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCt                                                                                                                                     >  5:62849/1‑148 (MQ=255)
        gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTGATGCCCTGGAGGCCGGCAGGGCCGGCCTccag                                                                                                                                 <  1:14887/147‑4 (MQ=255)
        gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTc                                                                                                                               >  1:16681/1‑149 (MQ=255)
        gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                  <  5:242060/146‑3 (MQ=255)
        gcagcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCAGGGCCGGCCTcca                                                                                                                                  >  6:242060/1‑144 (MQ=255)
          agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTc                                                                                                                               >  6:179547/1‑147 (MQ=255)
          agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                             <  6:62849/149‑1 (MQ=255)
          agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                             <  2:16681/149‑1 (MQ=255)
          agcaAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCa                                                                                                                             <  7:243629/149‑1 (MQ=255)
            caAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGt                                                                                                                           <  5:179547/149‑1 (MQ=255)
                   cgcCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccc                                                                                                                    >  2:30060/1‑149 (MQ=255)
                      cTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTccccc                                                                                                                  >  2:103225/1‑148 (MQ=255)
                        gTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCGCAATTCCAGTCTGTCCCCCCAg                                                                                                               >  6:50717/1‑149 (MQ=255)
                               ccAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCaa                                                                                                        >  6:83501/1‑149 (MQ=255)
                               ccAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTCTCCCACCAGGGCGCaa                                                                                                        >  7:94133/1‑149 (MQ=255)
                                           gCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATa                                                                                            <  1:81743/149‑1 (MQ=255)
                                             ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCGTGCCGGCGTCAATTCCAGTCTGT‑CCCCCACGGCGCAACGGCCGCGGATAg                                                                                           >  4:217032/1‑147 (MQ=255)
                                             ttACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                           >  4:180164/1‑148 (MQ=255)
                                             ttACGCCGCCCAAGGCGTTGCCAGCATCCCGCCGCACAGCCTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCACTGGAGTCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCACGGCGCAACGCCCccgctaag                                                                                           >  8:26709/1‑144 (MQ=255)
                                                 ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTATGTCCCCCCa                                                                                                                >  4:183407/1‑123 (MQ=255)
                                                 ggCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGGAATGCCCTGGAGGCCGGCCCTGCCGGCCGCAATTCCAGGCTGTCCCCCCAGGGCGCACCGCCCGCGCGTAGgagac                                                                                      >  4:196249/1‑148 (MQ=255)
                                                       cAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGcagaacagaa                                                                                 >  6:120825/1‑148 (MQ=255)
                                                         aGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCTGCGCATAGCAGaacagaaaca                                                                              <  1:103225/149‑1 (MQ=255)
                                                         aGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGGCGAGCTGCTCGACGTTCTGTACTGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGaacagaaaca                                                                              <  6:234388/149‑1 (MQ=255)
                                                          ggCGTTGTCACCATCCCGCCGCACAGCGTTCTGGTGTTCCACGCCTAGCTGCTCTACGTTCTGTAATTCCCTGGAGGCCGGCCCTGACGGCCTCAAATCCAGTCTGTCCCCCCAGTGCGCAACGCCCGCGCATAGCAGaacagaaacag                                                                             <  7:131734/149‑1 (MQ=255)
                                                             gTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGAtt                                                                          >  1:294024/1‑149 (MQ=255)
                                                                   agcaTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAc                                                                    <  1:30060/149‑1 (MQ=255)
                                                                     aaTCCCGCCCAACAGCGTTTTGGGTTGCGACGTGGTCCTGCTCTACGTTCTTTAATGCCCTGGCGTCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCa                                                                  <  8:170078/148‑1 (MQ=255)
                                                                      aTCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTAc                                                                    >  7:170078/1‑146 (MQ=255)
                                                                              gcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCAt                                                                                             <  6:86181/113‑1 (MQ=255)
                                                                              gcACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCAt                                                                                             >  5:86181/1‑113 (MQ=255)
                                                                                  aGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                           >  7:138770/1‑111 (MQ=255)
                                                                                  aGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATCCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAg                                                                                           <  8:138770/111‑1 (MQ=255)
                                                                                        cTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCg                                                <  5:83501/148‑1 (MQ=255)
                                                                                           gtgtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGtt                                             <  3:180164/148‑1 (MQ=255)
                                                                                           gtgtTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTc                                            >  3:309561/1‑149 (MQ=255)
                                                                                              ttCCACGTAGAGCTGCTCGACCTTCTGTAATGCCCTCGACGCCGGCCCTGCCGGCCTCAATTCCAGTCTGACACCCCAGGGCGCAACCCCCGCGCAAAGCAGAACAGAAACAGATTCCGTACCAATTCATGGAGCACCCCCGGTTCACt                                         >  4:267443/1‑149 (MQ=255)
                                                                                                gaaCGTCGAGCAGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTc                                        >  5:255204/3‑148 (MQ=255)
                                                                                                    gtctttctgctCGAAGTTCTGTAATGCCCTGGTGGCGTGCGCTTGGGGCCTCAATTCCAGTGTTTCCCCCCAGGGCGCAACGCCCGCGCATTTCAGAACAGTAAGAGATTCCGTACCAGTTCCTTGAGCACCCCCGTTTCACTCGAAtt                                   <  7:26709/143‑1 (MQ=255)
                                                                                                      cGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCa                                 <  4:201661/149‑1 (MQ=255)
                                                                                                          ctgatCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGACTGTCCCCCCAGGGCGCAACACCCGCGCATAGCAGAACAGCAACAGATTCCGAACCAGTTCCTAGAGCACCCCCGGTTAACTCGAATTCAGccc                             >  4:225861/1‑149 (MQ=255)
                                                                                                                gACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCAccc                                                   >  4:201565/1‑121 (MQ=255)
                                                                                                                gACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCAccc                                                   <  3:201565/121‑1 (MQ=255)
                                                                                                                  cGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGt                     <  5:265083/149‑1 (MQ=255)
                                                                                                                     tCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTc                    >  6:220944/1‑147 (MQ=255)
                                                                                                                      cTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGATACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAg                 >  7:236355/1‑149 (MQ=255)
                                                                                                                        gTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTACGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGt               >  5:187085/1‑149 (MQ=255)
                                                                                                                              cccTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCt         >  7:23988/1‑149 (MQ=255)
                                                                                                                                 tGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGt      >  3:28606/1‑149 (MQ=255)
                                                                                                                                    aGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAgg                <  7:282186/136‑1 (MQ=255)
                                                                                                                                    aGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAgg                >  8:282186/1‑136 (MQ=255)
                                                                                                                                     ggccggccCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACg  >  6:79101/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                            
GAACGCTGGCAGCAAATGGCGCCTGTACGTGCCAAGCGAGCTGGCTTACGGCGCCCAAGGCGTTGGCAGCATCCCGCCGCACAGCGTTCTGGTGTTCGACGTCGAGCTGCTCGACGTTCTGTAATGCCCTGGAGGCCGGCCCTGCCGGCCTCAATTCCAGTCTGTCCCCCCAGGGCGCAACGCCCGCGCATAGCAGAACAGAAACAGATTCCGTACCAGTTCCTTGAGCACCCCCGGTTCACTCGAATTCAGCCCGCTCATGTCCAGGTCACCCTGGTCACG  >  NC_002947/798522‑798803

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