Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,219,104 C→T 8.7% T214M (ACG→ATG)  PP_1064 → alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,219,1040CT 8.7% 127.6 / 2.5 46T214M (ACG→ATG) PP_1064alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (21/21);  new base T (2/2);  total (23/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 = 7.92e-01

CTACACGCTGACTGATATGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGG  >  NC_002947/1218958‑1219250
                                                                                                                                                  |                                                                                                                                                  
cTACACGCTGACTGATATGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGc                                                                                                                                                  <  1:234822/149‑1 (MQ=255)
     cGCTGACTGATATGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTc                                                                                                                                             <  4:312202/149‑1 (MQ=255)
        tgactgaTATGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATg                                                                                                                                          >  5:334125/1‑149 (MQ=255)
              atatGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCa                                                                                                                                      <  2:284257/147‑1 (MQ=255)
               tatGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGtt                                                                                                                                   <  4:108954/149‑1 (MQ=255)
                   gCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAgg                                                                                                                               >  6:188942/1‑149 (MQ=255)
                           cGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGGATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGcggcgg                                                                                                                       >  6:257987/1‑149 (MQ=255)
                           cGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGcggcgg                                                                                                                       <  4:171642/149‑1 (MQ=255)
                                       gATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGgcgc                                                                                                           >  5:295667/1‑149 (MQ=255)
                                                 ttGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGgc                                                                                                 <  1:249865/149‑1 (MQ=255)
                                                 ttGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGgc                                                                                                 <  6:295667/149‑1 (MQ=255)
                                                       gTGCGCCAGTTTCATGTCTTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                           <  3:294650/149‑1 (MQ=255)
                                                       gTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                           >  4:164748/1‑149 (MQ=255)
                                                          cgcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTg                                                                                          >  1:135853/1‑147 (MQ=255)
                                                             cAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                           >  1:134073/1‑143 (MQ=255)
                                                               gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGATGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCg                                                                                                               <  5:161231/121‑1 (MQ=255)
                                                               gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGATGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCg                                                                                                               <  5:166102/121‑1 (MQ=255)
                                                               gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGATGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCg                                                                                                               >  6:166102/1‑121 (MQ=255)
                                                               gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGATGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCg                                                                                                               >  6:161231/1‑121 (MQ=255)
                                                                tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCt                                                                                  >  3:247907/1‑149 (MQ=255)
                                                                     tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                             >  6:292178/1‑149 (MQ=255)
                                                                     tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                             <  6:334125/149‑1 (MQ=255)
                                                                       tATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCgg                                                                           >  2:233704/1‑149 (MQ=255)
                                                                          tGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGgcgc                                                                                                           >  2:203847/1‑114 (MQ=255)
                                                                          tGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGgcgc                                                                                                           <  1:203847/114‑1 (MQ=255)
                                                                            ggTGTGAGCATGGGCAGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcag                                                                      >  2:5561/1‑149 (MQ=255)
                                                                                    cATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCg                                                               >  8:249978/1‑148 (MQ=255)
                                                                                           ggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGa                                                       <  5:188942/149‑1 (MQ=255)
                                                                                           ggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGa                                                       <  5:292178/149‑1 (MQ=255)
                                                                                              aTGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGt                                                    >  5:163359/1‑149 (MQ=255)
                                                                                               tGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTg                                                   <  7:55460/149‑1 (MQ=255)
                                                                                               tGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTg                                                   <  7:272313/149‑1 (MQ=255)
                                                                                                      ccAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                             >  8:112713/1‑148 (MQ=255)
                                                                                                        aGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGag                                          <  5:257987/149‑1 (MQ=255)
                                                                                                             cTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAgcagca                                     <  7:249978/149‑1 (MQ=255)
                                                                                                                  ggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGcc                                 >  8:105017/1‑148 (MQ=255)
                                                                                                                          cgcCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACctgctg                        <  3:164748/149‑1 (MQ=255)
                                                                                                                                       gcgcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGt                                                         >  3:282221/1‑103 (MQ=255)
                                                                                                                                       gcgcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGt                                                         <  4:282221/103‑1 (MQ=255)
                                                                                                                                       gcgcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg             >  8:175244/1‑147 (MQ=255)
                                                                                                                                        cgcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAg          >  8:54691/1‑149 (MQ=255)
                                                                                                                                             ccTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGCACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCt      >  3:27201/1‑148 (MQ=255)
                                                                                                                                              cTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa    <  7:175244/149‑1 (MQ=255)
                                                                                                                                              cTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa    <  1:5561/149‑1 (MQ=255)
                                                                                                                                                gACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAgg  >  7:11447/1‑149 (MQ=255)
                                                                                                                                                 aCGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAg   <  7:54691/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTACACGCTGACTGATATGGCCGATGACGGCTTGCACCTGATGGACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGG  >  NC_002947/1218958‑1219250

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