Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,453,612 T→A 10.5% I503F (ATC→TTC)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,453,6120TA10.5% 103.1 / 2.9 38I503F (ATC→TTC) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base T (20/14);  new base A (2/2);  total (22/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 = 1.00e+00

CAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453464‑2453747
                                                                                                                                                    |                                                                                                                                       
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)
                                                                                                                                               gTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                        >  6:185890/1‑87 (MQ=255)
                                                                                                                                               gTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAg                                                        <  5:185890/87‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
CAGGCCTCCGAACCCAGCCGGTGCAGCGGTGCCAGGGCGTCGTCGCTGCCGGTGTAGCGGGCGATCAGGTGCAGGTTGGCCACCGGCACATACAGCTTGGCGTTCTCGGCGTATTCCAGCGTCAGGAATTCGGCGGCCTGGCCGTCGATTTCCAGCGTGGCCAGGCCCAGGTAGCGGCCGACGCCATGGTCGATGTGCACCACCGGCGCACCTTCGCGCAGCTCGGTGAGGTTCTTGATGACCGCGTCGTTGGCGGTTTCGCCGCGCTTGTCACGGCGCCGGCG  >  NC_002947/2453464‑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: ‑