Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,667,254 C→T 9.2% intergenic (+197/+65) PP_3230 → / ← PP_3231 phosphoribosyl transferase domain‑containing protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,667,2540CT 9.2% 150.7 / 4.7 54intergenic (+197/+65)PP_3230/PP_3231phosphoribosyl transferase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (2/3);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

CTCGTGGCGGGCGTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGCAAGGGCAGGGCAGTCACCGGCCGGGT  >  NC_002947/3667107‑3667399
                                                                                                                                                   |                                                                                                                                                 
cTCGTGGCGGGCGTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCt                                                                                                                                                  >  1:201783/1‑149 (MQ=255)
   gTGGCGGGCGTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAGCGCCGCGGCCTTCGc                                                                                                                                               <  1:311772/149‑1 (MQ=255)
   gTGGCGGGCGTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGc                                                                                                                                               >  2:220274/1‑149 (MQ=255)
      gcgggcgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCg                                                                                                                                             <  2:102926/148‑1 (MQ=255)
        gggcgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAGCGCCGCGGCCTTCGccaa                                                                                                                                            >  7:207869/1‑145 (MQ=255)
        gggcgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAGCGCCGCGGCCTTCGccaa                                                                                                                                            <  8:207869/147‑3 (MQ=255)
        gggcgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGt                                                                                                                                           >  5:86463/1‑148 (MQ=255)
         ggcgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAGCGCCGCGGCCTTCGccaaggt                                                                                                                                         <  2:100775/149‑6 (MQ=255)
           cgTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACt                                                                                                                                       <  6:222339/149‑1 (MQ=255)
            gTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGTTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTg                                                                                                                                      <  1:75652/149‑1 (MQ=255)
            gTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTg                                                                                                                                      <  6:68502/149‑1 (MQ=255)
              cTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGt                                                                                                                                    >  4:5566/1‑149 (MQ=255)
                gAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTc                                                                                                                                  >  2:178548/1‑149 (MQ=255)
                 aGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTc                                                                                                                                  <  2:201783/148‑1 (MQ=255)
                  gcgcCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGt                                                                                                                                >  6:137054/1‑149 (MQ=255)
                       cGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGt                                                                                                                           <  3:45168/149‑1 (MQ=255)
                         caatcaatACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAc                                                                                                                         <  3:5566/149‑1 (MQ=255)
                          aatcaatACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAcc                                                                                                                        >  3:229195/1‑149 (MQ=255)
                                       cAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATc                                                                                                           >  6:243346/1‑149 (MQ=255)
                                        aGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATcc                                                                                                          <  5:137054/149‑1 (MQ=255)
                                          aCGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGt                                                                                                        >  5:60908/1‑149 (MQ=255)
                                                  gcgGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATggg                                                                                                >  3:276661/1‑149 (MQ=255)
                                                   cgGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATgggg                                                                                               <  6:86463/149‑1 (MQ=255)
                                                       aCGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAc                                                                                                                         >  3:78343/1‑119 (MQ=255)
                                                       aCGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAc                                                                                                                         <  4:78343/119‑1 (MQ=255)
                                                       aCGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAgc                                                                                                                  <  3:73789/126‑1 (MQ=255)
                                                       aCGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAgc                                                                                                                  >  4:73789/1‑126 (MQ=255)
                                                                    tCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATa                                                                              <  3:131180/149‑1 (MQ=255)
                                                                        tactCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATacac                                                                           >  3:114322/1‑148 (MQ=255)
                                                                            ctTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAAt                                                                      >  8:289750/1‑149 (MQ=255)
                                                                             ttCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCaaa                                                                       <  2:81095/147‑1 (MQ=255)
                                                                                     gctgAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGtt                                                             >  3:27201/1‑149 (MQ=255)
                                                                                      ctgAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTTGCACCTTGGCGAAGGCCGCGGCACTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGttt                                                            >  1:219434/1‑149 (MQ=255)
                                                                                      ctgAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGttt                                                            <  4:114322/149‑1 (MQ=255)
                                                                                       tgAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTc                                                           >  1:174803/1‑149 (MQ=255)
                                                                                           cTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTcc                                                       <  4:229195/149‑1 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAcc                                                                                                                        <  5:317089/81‑1 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTAcc                                                                                                                        >  6:317089/1‑81 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACAcc                                                                          <  8:270733/127‑1 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACAcc                                                                          >  7:270733/1‑127 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCctct                                                    >  7:314/1‑149 (MQ=255)
                                                                                              aaTCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCctc                                                     <  6:60908/148‑1 (MQ=255)
                                                                                                tCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTgg                                                  <  5:243346/149‑1 (MQ=255)
                                                                                                  gTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCt                                                >  8:92901/1‑149 (MQ=255)
                                                                                                       gTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGt                                                              <  7:291421/130‑1 (MQ=255)
                                                                                                       gTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGt                                                              >  8:291421/1‑130 (MQ=255)
                                                                                                            ccGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCtg                                      <  1:178548/149‑1 (MQ=255)
                                                                                                             cGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCtgt                                     >  4:156000/1‑149 (MQ=255)
                                                                                                                   aaGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGTATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTc                               <  1:220274/149‑1 (MQ=255)
                                                                                                                    agagCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCg                              <  7:92901/149‑1 (MQ=255)
                                                                                                                      agCCTGGCACCTTGGCGAAGGCCGCGGCACTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGc                            <  2:219434/149‑1 (MQ=255)
                                                                                                                        ccTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGCaa                          <  3:145793/149‑1 (MQ=255)
                                                                                                                               caccggCGAGCGCCGCGGCCTTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTCATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCGCTGGCTCAAGATGCTGTTGGTTCGGAAagggcagg                   >  6:304352/5‑149 (MQ=255)
                                                                                                                                              cggcgCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGCAAGGGCAGGGCAGTCAccggccgg    <  8:314/149‑1 (MQ=255)
                                                                                                                                                gcgcTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGCAAGGGCAGGGCAGTCACCGGCCGGGt  <  7:289750/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCGTGGCGGGCGTCTGAGCGCCCGCAATCAATACGTAGCAGACGCTCTCGCGGCACGTGGCCTAAGGTCACTACTCTTCGACCTGCTGAGCTCAATCGTCAGGTGGACCGCTCGAAGAGCCTGGCACCTTGGCGAAGGCCGCGGCGCTCGCCGGTGACTGGTTCGTGAAGTACCTGCAGCGCATATCCGTCTGATGGGGCTTGGCCGAGGCTCATACACCAAATGGATAGGTTTCGTCCTCTGGCTCAAGATGCTGTTGTTTCGGCAAGGGCAGGGCAGTCACCGGCCGGGT  >  NC_002947/3667107‑3667399

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