Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,453,602 T→A 12.5% Q506L (CAG→CTG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,453,6020TA12.5% 100.3 / 6.6 40Q506L (CAG→CTG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base A (3/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453457‑2453747
                                                                                                                                                 |                                                                                                                                                 
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                                                                                                                                     >  5:298153/1‑141 (MQ=255)
                   cccAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTcc                                                                                                                                     <  6:298153/141‑1 (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:76319/1‑149 (MQ=255)
                                                       ccGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGAtgtg                                                                                         >  4:83254/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                                                                      <  1:55555/137‑1 (MQ=255)
                                                                                      ccACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTcg                                                                      >  2:55555/1‑137 (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)
                                                                                                                                                 |                                                                                                                                                 
TTTCTGCCAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453457‑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: ‑