Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,453,596 G→A 10.3% A508V (GCC→GTC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,453,5960GA10.3% 131.9 / 3.8 49A508V (GCC→GTC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (23/20);  new base A (2/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GCTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGG  >  NC_002947/2453450‑2453739
                                                                                                                                                    |                                                                                                                                               
gcTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTcggcg                                                                                                                                               >  8:88315/1‑149 (MQ=255)
gcTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTcggcg                                                                                                                                               <  6:337188/149‑1 (MQ=255)
       tttCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGc                                                                                                                                        >  8:222420/1‑149 (MQ=255)
       tttCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCCTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGc                                                                                                                                        >  1:238303/1‑149 (MQ=255)
           tGCCAGGCCTCCGAACGCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACTGCTTGGCGTTCTCGGCGTTTTCCAGCGTCAGGAAATCGACGGCCAGGCCgcc                                                                                                                                    <  1:24149/149‑3 (MQ=255)
           tGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGcc                                                                                                                                            >  1:312659/1‑141 (MQ=255)
           tGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGcc                                                                                                                                            <  2:312659/141‑1 (MQ=255)
             ccAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCg                                                                                                                                   <  7:82531/148‑1 (MQ=255)
               aGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGAtt                                                                                                                                <  2:66046/149‑1 (MQ=255)
                 gCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGCCGCCGAATTc                                                                                                                              <  7:222420/149‑1 (MQ=255)
                  ccTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGGGTAGCGGGCGATCAGGGGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTcc                                                                                                                             <  7:70666/149‑1 (MQ=255)
                   cTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTc                                                                                                                              >  8:62526/1‑147 (MQ=255)
                           ccAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGCCGCCGaa                                                                                                                                 <  5:214138/136‑2 (MQ=255)
                           ccAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGCCGCCGaa                                                                                                                                 >  6:214138/1‑135 (MQ=255)
                                      aGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCa                                                                                                                   <  7:13243/139‑1 (MQ=255)
                                      aGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCa                                                                                                                   >  8:13243/1‑139 (MQ=255)
                                          gTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCgg                                                                                                     >  5:246040/1‑149 (MQ=255)
                                               aGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGAc                                                                                                >  4:13310/1‑149 (MQ=255)
                                                       gtcgCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGt                                                                                        >  2:66599/1‑149 (MQ=255)
                                                        tcgCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTc                                                                                       <  7:233896/149‑1 (MQ=255)
                                                         cgCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTc                                                                                                                                    >  6:14507/1‑103 (MQ=255)
                                                           cTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAt                                                                                    <  3:17568/149‑1 (MQ=255)
                                                           cTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAt                                                                                    <  8:174575/149‑1 (MQ=255)
                                                                 gtgtAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGcac                                                                              >  4:227234/1‑149 (MQ=255)
                                                                  tgtAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGcac                                                                              >  6:337617/1‑148 (MQ=255)
                                                                   gtAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGcacca                                                                            >  1:68032/1‑149 (MQ=255)
                                                                          gCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCg                                                                                      <  4:266698/132‑1 (MQ=255)
                                                                          gCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCg                                                                                      >  3:266698/1‑132 (MQ=255)
                                                                                aGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTc                                                               <  2:68032/149‑1 (MQ=255)
                                                                                 gctGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGCCGCCCTGCCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAAGTGCACAACCGGCGCACATTcg                                                              >  2:322108/3‑149 (MQ=255)
                                                                                  gTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcgc                                                             <  1:318745/149‑1 (MQ=255)
                                                                                    tgccggtGGCCAACCTGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCTACGCCATGGTCGATGTGCACCACCGGCTCACCTTcgcg                                                            >  8:212656/7‑149 (MQ=255)
                                                                                    gCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcgcg                                                            >  1:234415/1‑148 (MQ=255)
                                                                                        gTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCt                                                       >  7:59397/1‑149 (MQ=255)
                                                                                         ttGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTc                                                      >  3:261702/1‑149 (MQ=255)
                                                                                              cACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCGCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGt                                                   >  1:227568/1‑148 (MQ=255)
                                                                                                aCCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                <  3:13310/149‑1 (MQ=255)
                                                                                                aCCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                <  4:19696/149‑1 (MQ=255)
                                                                                                 ccGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAgg                                               <  7:88315/149‑1 (MQ=255)
                                                                                                      acaTACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCtt                                          <  8:59397/149‑1 (MQ=255)
                                                                                                         tACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTtgat                                       <  5:337617/149‑1 (MQ=255)
                                                                                                          aCAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTtgat                                       >  7:179446/1‑148 (MQ=255)
                                                                                                          aCAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTtgat                                       >  1:332333/1‑148 (MQ=255)
                                                                                                                     ttCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGtt                           <  2:238303/149‑1 (MQ=255)
                                                                                                                          ggCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTgg                         <  6:246040/146‑1 (MQ=255)
                                                                                                                                  ccAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTcgca               <  7:212656/148‑2 (MQ=255)
                                                                                                                                    cAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTcgccgc             >  5:228179/1‑149 (MQ=255)
                                                                                                                                     aGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCTCCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCcgcg            >  3:301786/1‑149 (MQ=255)
                                                                                                                                            ggAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCa     >  3:295930/1‑149 (MQ=255)
                                                                                                                                              aaTTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAcg   <  3:227234/149‑1 (MQ=255)
                                                                                                                                               aTTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAcgg  >  3:114774/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGG  >  NC_002947/2453450‑2453739

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: