Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,040,850 T→A 12.5% intergenic (+19/+15) PP_5292 → / ← PP_5293 catabolite repression control protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,040,8500TA12.5% 147.8 / 10.9 55intergenic (+19/+15)PP_5292/PP_5293catabolite repression control protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/24);  new base A (3/4);  total (26/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 = 1.00e+00

GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGA  >  NC_002947/6040705‑6040995
                                                                                                                                                 |                                                                                                                                                 
gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                               >  2:1065114/1‑150 (MQ=255)
gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                               >  2:7498/1‑150 (MQ=255)
gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                               >  1:1123207/1‑150 (MQ=255)
gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                               <  1:7498/150‑1 (MQ=255)
gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCa                                                                                                                                               >  2:486795/1‑150 (MQ=255)
 ttCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCGTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAg                                                                                                                                              >  1:1010823/1‑150 (MQ=255)
       tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCaa                                                                                                                                        <  2:174853/150‑1 (MQ=255)
       tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCaa                                                                                                                                        <  2:724155/150‑1 (MQ=255)
         ccAGATCCTCACTCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                      <  2:970331/150‑1 (MQ=255)
         ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg                                                                                                                                      >  1:326208/1‑150 (MQ=255)
         ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg                                                                                                                                      >  1:778506/1‑150 (MQ=255)
         ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                      >  2:61738/1‑150 (MQ=255)
         ccAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                      <  1:61738/150‑1 (MQ=255)
          cAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATgg                                                                                                                                     >  1:164993/1‑150 (MQ=255)
          cAGATCCTCACCACGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTGCGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCCATgg                                                                                                                                     >  1:137086/1‑150 (MQ=255)
                cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                >  1:354205/1‑149 (MQ=255)
                cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                <  2:354205/149‑1 (MQ=255)
                cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGCGCCTGCATTGCAGCCCCAATGGTTAttt                                                                                                                               <  1:759326/150‑1 (MQ=255)
                 tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTg                                                                                                                              >  1:772055/1‑150 (MQ=255)
                     cccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                          >  2:643648/1‑150 (MQ=255)
                     cccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                          <  1:812005/150‑1 (MQ=255)
                       cGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgc                                                                                                                        <  2:208178/150‑1 (MQ=255)
                            cTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCa                                                                                                                   <  2:952270/150‑1 (MQ=255)
                                gccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                           <  1:740261/122‑1 (MQ=255)
                                gccgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                           >  2:740261/1‑122 (MQ=255)
                                  cgcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcat                                                                                                             >  1:160639/1‑150 (MQ=255)
                                   gcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                            <  2:1060076/150‑1 (MQ=255)
                                   gcTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGatcatc                                                                                                            >  1:1113574/1‑150 (MQ=255)
                                            gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTANCGa                                                                                                   <  2:699755/150‑1 (MQ=255)
                                            gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGa                                                                                                   >  2:719857/1‑150 (MQ=255)
                                                cGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGACCATCGTGGAGCGATAAg                                                                                               <  2:1123207/150‑1 (MQ=255)
                                                  cTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGc                                                                                             <  1:725671/150‑1 (MQ=255)
                                                  cTCGCCTGCCGCGACAGCCACGCTTCTCCCAGCATGCACCGCTGACCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCAGTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGc                                                                                             <  1:500951/150‑1 (MQ=255)
                                                   tCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCt                                                                                            >  2:1048847/1‑150 (MQ=255)
                                                       cTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGCCCGCCAGATCATCGGGTAGCGATAAGGCTTGcc                                                                                        <  1:486795/150‑1 (MQ=255)
                                                           cgcgTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCa                                                                                    <  2:518228/150‑1 (MQ=255)
                                                                 aGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGcc                                                                              <  2:326208/150‑1 (MQ=255)
                                                                    cACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGc                                                                           >  2:720240/1‑150 (MQ=255)
                                                                               cAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCgat                                                                >  2:337805/1‑150 (MQ=255)
                                                                                         cGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat                                                      <  2:1113574/150‑1 (MQ=255)
                                                                                          gCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAata                                                     <  1:1065114/150‑1 (MQ=255)
                                                                                          gCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGACCTCATTCGCCCGCACCGCCGCGATGATGATCAata                                                     <  2:457088/150‑1 (MQ=255)
                                                                                                  tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACt                                             >  2:766976/1‑150 (MQ=255)
                                                                                                  tcgACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCACAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACt                                             >  2:599208/1‑150 (MQ=255)
                                                                                                    gACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGc                                           <  1:1048847/150‑1 (MQ=255)
                                                                                                                aCGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCAt                               <  1:766976/150‑1 (MQ=255)
                                                                                                                      aCCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCaaa                         <  1:337805/150‑1 (MQ=255)
                                                                                                                          tCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGc                     <  2:1010823/150‑1 (MQ=255)
                                                                                                                                    cATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTc                                                                                 <  2:202708/80‑1 (MQ=255)
                                                                                                                                    cATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATAATCGGGTAGCGATAAGGCTTGCCCTCATTc                                                                                 >  1:202708/1‑80 (MQ=255)
                                                                                                                                     aTTGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGc                                                                               <  2:549086/81‑1 (MQ=255)
                                                                                                                                     aTTGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGc                                                                               >  1:549086/1‑81 (MQ=255)
                                                                                                                                            cTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCAt                                                                                   <  1:1089259/70‑1 (MQ=255)
                                                                                                                                            cTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCAt                                                                                   >  2:1089259/1‑70 (MQ=255)
                                                                                                                                             tGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGa  >  2:602893/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGCCGCCAGATCATCGGGTAGCGATAAGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGA  >  NC_002947/6040705‑6040995

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