Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,651,428 C→A 13.2% intergenic (‑185/+28) PP_4959 ← / ← fba diguanylate cyclase/fructose‑bisphosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,651,4280CA13.2% 174.5 / 10.9 59intergenic (‑185/+28)PP_4959/fbadiguanylate cyclase/fructose‑bisphosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (27/24);  new base A (4/4);  total (31/28)
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

ATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGC  >  NC_002947/5651285‑5651569
                                                                                                                                                |                                                                                                                                             
aTCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                          >  2:36180/1‑149 (MQ=255)
aTCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCgg                                                                                                                                          <  1:36180/149‑1 (MQ=255)
aTCCAATAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAAACAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCAGCATCGCTGCggg                                                                                                                                         >  2:539955/1‑150 (MQ=255)
              tAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTAc                                                                                                                           <  1:1168454/150‑1 (MQ=255)
                  ccTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGATGCGGGTTTTTCAGGGTTACCTGt                                                                                                                       <  2:1086138/150‑1 (MQ=255)
                   cTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTc                                                                                                                      >  1:712361/1‑150 (MQ=255)
                     gAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGa                                                                                                                    <  2:712361/150‑1 (MQ=255)
                        gCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 <  2:28887/134‑1 (MQ=255)
                        gCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 >  1:28887/1‑134 (MQ=255)
                        gCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCa                                                                                                                 >  2:955425/1‑150 (MQ=255)
                        gCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCa                                                                                                                 >  1:483953/1‑150 (MQ=255)
                           gCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTATCGATCAGtt                                                                                                              >  2:47898/1‑150 (MQ=255)
                                acGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCt                                                                                                         >  2:520809/1‑150 (MQ=255)
                                  gcaatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTg                                                                                                       >  2:505524/1‑150 (MQ=255)
                                    aatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 <  2:576632/122‑1 (MQ=255)
                                    aatgcaatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAg                                                                                                                                 >  1:576632/1‑122 (MQ=255)
                                         aatCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGCACTTATGTACAGACTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCa                                                                                                >  2:1041340/1‑150 (MQ=255)
                                               aaCGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCAc                                                                                          <  1:535161/150‑1 (MQ=255)
                                                           gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                              <  1:769782/118‑1 (MQ=255)
                                                           gATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGtt                                                                                                              >  2:769782/1‑118 (MQ=255)
                                                               aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTg                                                                          >  1:376230/1‑150 (MQ=255)
                                                               aaaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTg                                                                          <  1:746679/150‑1 (MQ=255)
                                                                aaCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         <  2:27253/150‑1 (MQ=255)
                                                                 aCTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCTGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACTTTGACGGCCAGTTCACCTTTCAGGTAGCGCTGGt                                                                        >  2:591670/1‑150 (MQ=255)
                                                                  cTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgcg                                                                             <  2:359626/144‑1 (MQ=255)
                                                                  cTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAgcg                                                                             >  1:359626/1‑144 (MQ=255)
                                                                   tGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGt                                                                                                                                        <  1:244058/84‑1 (MQ=255)
                                                                   tGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGt                                                                                                                                        >  2:244058/1‑84 (MQ=255)
                                                                       cGGTATCCAACTAACGGTCGGTCTGCCGGCATTATGTACAGGTTTCCATATCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGc                                                                  <  2:1103504/150‑1 (MQ=255)
                                                                           aGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTc                                                              <  2:303237/150‑1 (MQ=255)
                                                                             cccaCTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCACCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGCCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCa                                                            <  1:131727/147‑1 (MQ=255)
                                                                               gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCt                                                                                                                         <  1:109432/87‑1 (MQ=255)
                                                                               gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCt                                                                                                                         >  2:109432/1‑87 (MQ=255)
                                                                               gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttt                                                                                                                                   <  2:869957/77‑2 (MQ=38)
                                                                               gACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCAGCGATGCGGGtttttt                                                                                                                                   >  1:869957/1‑76 (MQ=38)
                                                                                         cggtcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCt                                                <  2:92254/150‑1 (MQ=255)
                                                                                            tcTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGa                                             >  1:837692/1‑150 (MQ=255)
                                                                                               gCCGGAATTATGTACAGGTTTCAATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         >  1:855710/1‑119 (MQ=255)
                                                                                               gCCGGAATTATGTACAGGTTTCAATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTgg                                                                         <  2:855710/119‑1 (MQ=255)
                                                                                                 cGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTc                                        >  2:225669/1‑150 (MQ=255)
                                                                                                  ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCt                                                                >  1:855408/1‑125 (MQ=255)
                                                                                                  ggAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCt                                                                <  2:855408/125‑1 (MQ=255)
                                                                                                        aTGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCAt                                 <  2:483953/150‑1 (MQ=255)
                                                                                                           tACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGc                              >  1:1162802/1‑150 (MQ=255)
                                                                                                             cAGGTTTCCATCTCCCATAAAAAAAACCCGCAGCGATGCGGGt                                                                                                                                        <  2:15950/43‑1 (MQ=17)
                                                                                                             cAGGTTTCCATCTCCCATAAAAAAAACCCGCAGCGATGCGGGt                                                                                                                                        >  1:15950/1‑43 (MQ=17)
                                                                                                                gTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCcgcg                         >  2:44866/1‑150 (MQ=255)
                                                                                                                        ctcCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGaa                 <  1:44866/150‑1 (MQ=255)
                                                                                                                          cccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAgg               >  1:285698/1‑150 (MQ=255)
                                                                                                                           ccATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGc              <  1:468536/150‑1 (MQ=255)
                                                                                                                             aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                     <  1:1072006/125‑1 (MQ=255)
                                                                                                                             aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAg                                     >  2:1072006/1‑125 (MQ=255)
                                                                                                                             aTAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGACCTTCGAGACCATCAGCGCGGTGCCGAAGGCtt            >  1:62105/1‑150 (MQ=255)
                                                                                                                                tgaaaaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg          >  2:294321/3‑150 (MQ=255)
                                                                                                                                aaaaaaaCCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCg          <  1:505524/150‑1 (MQ=255)
                                                                                                                                    aaaCCCGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAc      >  1:883857/1‑150 (MQ=255)
                                                                                                                                        ccGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc  >  1:30357/1‑150 (MQ=255)
                                                                                                                                        ccGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc  <  1:1041340/150‑1 (MQ=255)
                                                                                                                                        ccGCAGCGATGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAAcgcgc  <  1:294321/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
ATCCATTAGTACTCTAGACCTGAGGCTGCAACACGCAATGCAATCAGAACGCTATCAATGATCAAACTGTACGGTAGCCGACTAACGGTCGGTCTGCCGGAATTATGTACAGGTTTCCATCTCCCATAAAAAAACCCGCATCGCTGCGGGTTTTTCAGGGTTACCTGTCGATCAGTTGACCTTGGCGGCCAGTTCACCTTTCAGGTAGCGCTGGTACATGCCTTCCAGAGAGATCGGCTTGATCTTCGAGGCATTGCCCGCGGTGCCGAAGGCTTCGTAACGCGC  >  NC_002947/5651285‑5651569

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