Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,453,596 G→A 14.7% 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,5960GA14.7% 78.0 / 6.0 34A508V (GCC→GTC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (18/11);  new base A (3/2);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

CGCTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCC  >  NC_002947/2453449‑2453743
                                                                                                                                                   |                                                                                                                                                   
cgcTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTcggc                                                                                                                                                    >  1:144245/1‑149 (MQ=255)
        tttCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGACGGCCAGGc                                                                                                                                            >  6:315059/1‑149 (MQ=255)
             gCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCACGTTGGCCACCCGCACATACAGCTTGGCGTTCTAGGCGTAAACGAGCGTCAGGAATTCGGCGGCCTGGcgggcg                                                                                                                                       >  8:66703/1‑144 (MQ=255)
              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                                                                                    >  5:197207/1‑118 (MQ=255)
                                                                                               cACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGc                                                                                    <  6:197207/118‑1 (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                       <  6:225218/135‑1 (MQ=255)
                                                                                                                                          tggAAATCGACGGCCAGGCCGCCGAATTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGttt                       >  5:225218/2‑136 (MQ=255)
                                                                                                                                          tggAAATCGACGGCCAGGCCGCCGAATTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGt           >  3:89938/2‑148 (MQ=255)
                                                                                                                                            gAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCAGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAc        >  2:62011/1‑149 (MQ=255)
                                                                                                                                                  ggcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCAcggcgcc  >  8:83082/1‑149 (MQ=255)
                                                                                                                                                   gcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCg                                                         <  5:96723/93‑1 (MQ=255)
                                                                                                                                                   gcggcCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCg                                                         >  6:96723/1‑93 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGCTTGGCTTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCC  >  NC_002947/2453449‑2453743

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