Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,453,608 T→C 10.0% D504G (GAC→GGC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,453,6080TC10.0% 103.3 / 3.6 40D504G (GAC→GGC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (21/15);  new base C (2/2);  total (23/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

CCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453463‑2453747
                                                                                                                                                 |                                                                                                                                           
ccAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGCCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGa                                                                                                                                          >  8:292666/1‑149 (MQ=255)
 cAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGGAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGAt                                                                                                                                         <  5:76152/149‑1 (MQ=255)
    gCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGCCGCCGAATTc                                                                                                                                      <  5:315059/149‑1 (MQ=255)
             cccAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTcc                                                                                                                                     <  6:298153/141‑1 (MQ=255)
             cccAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTcc                                                                                                                                     >  5:298153/1‑141 (MQ=255)
                              tGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGc                                                                                                            <  8:211299/149‑1 (MQ=255)
                                          gtcgCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGt                                                                                                >  1:140597/1‑149 (MQ=255)
                                           tcgCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTc                                                                                               <  1:254469/149‑1 (MQ=255)
                                                 ccGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAtgtg                                                                                         >  4:83254/1‑149 (MQ=255)
                                                 ccGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAtgtg                                                                                         >  4:76319/1‑149 (MQ=255)
                                                         gCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTCGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCg                                                                                 >  1:17581/1‑149 (MQ=255)
                                                           gggCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGgc                                                                               >  1:129916/1‑149 (MQ=255)
                                                            ggCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATTGTCGATGTGCACCACCGgcg                                                                              <  7:292666/149‑1 (MQ=255)
                                                                     gTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcgc                                                                     <  3:83254/149‑1 (MQ=255)
                                                                              ggCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCgg                                                            >  3:218498/1‑149 (MQ=255)
                                                                                ccACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcg                                                                      >  2:55555/1‑137 (MQ=255)
                                                                                ccACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcg                                                                      <  1:55555/137‑1 (MQ=255)
                                                                                ccACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGt                                                           >  1:157632/1‑148 (MQ=255)
                                                                                 cACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGc                                                                                        <  6:197207/118‑1 (MQ=255)
                                                                                 cACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGc                                                                                        >  5:197207/1‑118 (MQ=255)
                                                                                  aCCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                        <  2:140597/149‑1 (MQ=255)
                                                                                             cAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTtgatga                                             >  8:191270/1‑149 (MQ=255)
                                                                                                          tcGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGc                                >  5:140475/1‑149 (MQ=255)
                                                                                                           cGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCg                               >  3:72056/1‑149 (MQ=255)
                                                                                                                    ccAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTcgccg                      >  1:144874/1‑149 (MQ=255)
                                                                                                                      aGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCCTCGTTGGCGGTTTCGCcgcg                    <  2:17581/149‑1 (MQ=255)
                                                                                                                       gCGTCAGGAATTCGGCGGCCTGGCCGTGGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGGAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCgcgc                   <  2:144245/149‑1 (MQ=255)
                                                                                                                          tCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTg                >  8:283331/1‑149 (MQ=255)
                                                                                                                            tggAAATCGACGGCCAGGCCGCCGAATTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGttt                           >  5:225218/2‑136 (MQ=255)
                                                                                                                            tggAAATCGACGGCCAGGCCGCCGAATTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGttt                           <  6:225218/135‑1 (MQ=255)
                                                                                                                            tggAAATCGACGGCCAGGCCGCCGAATTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGt               >  3:89938/2‑148 (MQ=255)
                                                                                                                              gAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCAGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAc            >  2:62011/1‑149 (MQ=255)
                                                                                                                                    ggcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAcggcgcc      >  8:83082/1‑149 (MQ=255)
                                                                                                                                     gcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCg                                                             >  6:96723/1‑93 (MQ=255)
                                                                                                                                     gcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCg                                                             <  5:96723/93‑1 (MQ=255)
                                                                                                                                        gcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCtt                 >  6:316929/1‑134 (MQ=255)
                                                                                                                                        gcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCtt                 <  5:316929/134‑1 (MQ=255)
                                                                                                                                        gcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAcggcgccggcg  <  4:218498/149‑1 (MQ=255)
                                                                                                                                                gTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                        >  6:185890/1‑87 (MQ=255)
                                                                                                                                                gTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                        <  5:185890/87‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
CCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453463‑2453747

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