Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,409,965 G→A 8.0% A304V (GCC→GTC)  stcD ← N‑methylproline demethylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,409,9650GA 8.0% 143.7 / 2.6 50A304V (GCC→GTC) stcDN‑methylproline demethylase
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (2/2);  total (27/23)
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.17e-01

CCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTGGCGATGTTGCC  >  NC_002947/5409819‑5410107
                                                                                                                                                  |                                                                                                                                              
ccAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGcc                                                                                                                                              >  8:123392/1‑149 (MQ=255)
  aGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCACCCCCGAGGCAATGGCATGACGCGCCGTGGCCAc                                                                                                                                            >  8:234660/1‑149 (MQ=255)
      aaTAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCg                                                                                                                                        <  5:46453/149‑1 (MQ=255)
                gcCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTtg                                                                                                                              >  1:191123/1‑149 (MQ=255)
                  cTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTtgct                                                                                                                            <  5:123606/149‑1 (MQ=255)
                        aaGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATg                                                                                                                       >  1:153219/1‑148 (MQ=255)
                          gggCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGGCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGaa                                                                                                                    <  5:250804/149‑1 (MQ=255)
                             cGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGt                                                                                                                 >  3:35521/1‑149 (MQ=255)
                              ggATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAAtttt                                                                                                                                <  6:170911/133‑1 (MQ=255)
                              ggATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAAtttt                                                                                                                                >  5:170911/1‑133 (MQ=255)
                                     gTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGc                                                                                                         >  3:18978/1‑149 (MQ=255)
                                     gTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGc                                                                                                         >  5:243389/1‑149 (MQ=255)
                                           aCGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTc                                                                                                    <  4:18978/148‑1 (MQ=255)
                                                    gatgatCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCttt                                                                                           <  2:191123/148‑1 (MQ=255)
                                                    gatgatCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTa                                                                                          >  5:245139/1‑149 (MQ=255)
                                                      tgatCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTAtt                                                                                        <  7:123392/149‑1 (MQ=255)
                                                         tCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGc                                                                                                         <  1:147443/129‑1 (MQ=255)
                                                         tCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGc                                                                                                         >  2:147443/1‑129 (MQ=255)
                                                         tCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGAGCCGTGACCACGTCGGGAATTTATGCTCCATGGAAGGTCGGGAAGACGGTCAGCGTCTTTATTTCg                                                                                     >  5:362336/1‑149 (MQ=39)
                                                                 aCGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCaaa                                                                             <  4:35521/149‑1 (MQ=255)
                                                                        gCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTtg                                                                                                                              <  4:105884/93‑1 (MQ=255)
                                                                        gCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTtg                                                                                                                              >  3:105884/1‑93 (MQ=255)
                                                                               ggtcATGTGCGCGCGGGTCATGCCGACCATGTCCACCACCCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCgggg                                                               >  2:336461/1‑149 (MQ=255)
                                                                                gtcATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTc                                                                                      <  1:205893/125‑1 (MQ=255)
                                                                                gtcATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTc                                                                                      >  2:205893/1‑125 (MQ=255)
                                                                                 tcATGTGCGCGCGTGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGGCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGc                                                              <  2:153219/148‑1 (MQ=255)
                                                                                 tcATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCa                                                             >  1:463/1‑149 (MQ=255)
                                                                                      tgCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTgcg                                                        <  1:336461/149‑1 (MQ=255)
                                                                                      tgCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTgc                                                         <  6:245139/148‑1 (MQ=255)
                                                                                                  aTGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGc                                                              <  2:91302/131‑1 (MQ=255)
                                                                                                  aTGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGc                                                              >  1:91302/1‑131 (MQ=255)
                                                                                                  aTGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGa                                             <  6:243389/148‑1 (MQ=255)
                                                                                                   tGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTc                                                                          >  6:182492/1‑118 (MQ=255)
                                                                                                   tGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTc                                                                          <  5:182492/118‑1 (MQ=255)
                                                                                                    gCCGACCATGACCACCAACCCCGATGCAATGGCATGACGCGCCGTGACAACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCTCCGGAAAAGTCCAGCTTCGGGGCATTACGCATACCCAGGATCg                                          >  8:4246/1‑149 (MQ=255)
                                                                                                       gACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGga                                       <  4:307318/149‑1 (MQ=255)
                                                                                                          cATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGgatga                                    >  4:109137/1‑149 (MQ=255)
                                                                                                               ccaccaACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCg                               >  1:162095/1‑149 (MQ=255)
                                                                                                                caccaACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGgatgat                                   <  5:32250/144‑1 (MQ=255)
                                                                                                                caccaACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGgatgat                                   >  6:32250/1‑144 (MQ=255)
                                                                                                                caccaACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCg                               <  2:162095/148‑1 (MQ=255)
                                                                                                                        ccGAGGCAATGGCATGACGCTCCGTGGCCACTTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGACGGTCAGCGTCTTTATTTCGCCGGCAAAGTCAAGGTGCGGGGCATTGCTCATACCATGAATCGGGATGATGTCGGTCAGccct                      >  4:52537/1‑148 (MQ=255)
                                                                                                                            ggCAATGGCATGACGCGCCGTGGCCACGTCGGGAATGTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGt                  <  6:174734/149‑1 (MQ=255)
                                                                                                                              cAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCg                >  4:73739/1‑149 (MQ=255)
                                                                                                                               aaTGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCgg               >  3:29498/1‑149 (MQ=255)
                                                                                                                                   gCATGACGCGCCGTGACCACGTTGTGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGTTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTgg            <  6:362336/148‑1 (MQ=255)
                                                                                                                                    cATGACGCGCCGTGACCACGTCGTGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCTGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTGGCg          <  7:4246/149‑1 (MQ=255)
                                                                                                                                     aTGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTcgcg          >  1:127473/1‑148 (MQ=255)
                                                                                                                                       ggcgcgCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTGGCGATg       >  4:179701/3‑149 (MQ=255)
                                                                                                                                            gcCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTGGCGATGTTGcc  >  6:267527/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGGTCATGTGCGCGCGGGTCATGCCGACCATGTCCACCAACCCCGAGGCAATGGCATGACGCGCCGTGGCCACGTCGGGAATTTTTGCTGCATGGAAGGTCGGGAAGCCGGTCAGCGTCTTTATTTCGCCGGCAAAGTCCAGGTGCGGGGCATTGCGCATACCCTGGATCGGGATGATGTCGGTCAGCCCGGCGTCGGTGGCGATGTTGCC  >  NC_002947/5409819‑5410107

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