Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,011,276 A→T 8.5% V53D (GTT→GAT)  PP_1792 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,011,2760AT 8.5% 133.5 / 3.4 47V53D (GTT→GAT) PP_1792group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (21/22);  new base T (2/2);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GTATTAGCAATTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGCTCGTTGAGGACGTT  >  NC_002947/2011132‑2011424
                                                                                                                                                 |                                                                                                                                                    
gTATTAGCAATTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCATCg                                                                                                                                                    >  7:198798/1‑147 (MQ=255)
   ttAGCAATTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTTGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGaaaa                                                                                                                                                <  6:286109/148‑1 (MQ=255)
      gCAATTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGATACTGGCCCCGTCAACGAAAACTTc                                                                                                                                            >  3:9715/1‑149 (MQ=255)
        actttGAGCGGTTATTCGGTTCGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGa                                                                                                                                          <  8:264134/147‑1 (MQ=255)
         aTTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGAt                                                                                                                                         <  5:322332/149‑1 (MQ=255)
              aGCGGTTAATCGGCTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAACGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACa                                                                                                                                     <  8:328248/148‑1 (MQ=255)
               gCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCATCGAAGTTTTCGTTGACggg                                                                                                                                   <  8:198798/149‑1 (MQ=255)
                                 tcctccAGCAGGCTCTGTCGCAGTTTCACGGTAAAGGCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGt                                                                                                                 <  8:243544/149‑1 (MQ=255)
                                 tCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGt                                                                                                                 >  8:282550/1‑149 (MQ=255)
                                   cGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTg                                                                                                                <  2:30105/148‑1 (MQ=255)
                                    gCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGat                                                                                                              <  1:244234/149‑1 (MQ=255)
                                                   cgcAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCAccc                                                                                               >  6:304260/1‑149 (MQ=255)
                                                    gcAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCAcccc                                                                                              >  1:307036/1‑149 (MQ=255)
                                                     cAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCt                                                                                             <  1:246759/149‑1 (MQ=255)
                                                          gCACGGTAAAGCCATGGAAGTGATCGGCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTg                                                                                            >  1:186747/1‑146 (MQ=255)
                                                           cACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAAc                                                                                                                                               >  7:193243/1‑93 (MQ=255)
                                                           cACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAAc                                                                                                                                               <  8:193243/93‑1 (MQ=255)
                                                            aCGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGaa                                                                                       >  3:268820/1‑148 (MQ=255)
                                                              ggTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCAGTGGGCTCATGTTGGTCTGCCCTGGCCCATGCGCCACTGGACCCGTCAACGACAACTTCGATGACAGGGCGCTGATATGGCTCCGTGATAGCGATAGCCCACCCCTGCAGAAAAt                                                                                    >  2:279514/1‑149 (MQ=255)
                                                                          gAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGATACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGt                                                                         <  4:9715/148‑1 (MQ=255)
                                                                           aaGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTcc                                                                       >  4:183741/1‑149 (MQ=255)
                                                                             gTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGt                                                                     <  7:282550/149‑1 (MQ=255)
                                                                                 tC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGAt                                                                 >  8:251041/1‑149 (MQ=255)
                                                                                    gCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATg                                                                >  4:216575/1‑148 (MQ=255)
                                                                                                aTTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGaaaa                                                    <  7:251041/148‑1 (MQ=255)
                                                                                                 ttGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAaga                                                  <  2:186747/149‑1 (MQ=255)
                                                                                                   gggCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGtt                                                                    <  5:168208/129‑1 (MQ=255)
                                                                                                   gggCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGtt                                                                    >  6:168208/1‑129 (MQ=255)
                                                                                                   gggCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 <  5:304260/148‑1 (MQ=255)
                                                                                                    ggCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAgagga                                               <  3:216575/149‑1 (MQ=255)
                                                                                                     gCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 >  3:252669/1‑146 (MQ=255)
                                                                                                       tCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGagag                                            >  6:357955/1‑149 (MQ=255)
                                                                                                                  tGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGTTATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAgcg                                          <  6:361574/140‑1 (MQ=255)
                                                                                                                  tGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAgcg                                          >  5:361574/1‑140 (MQ=255)
                                                                                                                    cccTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCATGCAGAACAAACCCGTAAATTCCGTTGAAGCAGCCACGAAAAGAGGAGAGCGCTACAGGGTTc                               >  4:330147/1‑149 (MQ=255)
                                                                                                                       tGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 >  7:143376/1‑128 (MQ=255)
                                                                                                                       tGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 <  8:143376/128‑1 (MQ=255)
                                                                                                                         gCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGa                            <  2:315451/147‑1 (MQ=255)
                                                                                                                            aaTGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCTGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACccgc                       >  7:172247/1‑149 (MQ=255)
                                                                                                                                cGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCgccggc                   <  2:307036/149‑1 (MQ=255)
                                                                                                                                 gACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCgccggcc                  >  7:159976/1‑149 (MQ=255)
                                                                                                                                  acacTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCAGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCgccggccg                 >  1:264879/1‑149 (MQ=255)
                                                                                                                                   cacTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGc                >  3:311535/1‑149 (MQ=255)
                                                                                                                                        gCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGCTCGtt           <  4:304041/149‑1 (MQ=255)
                                                                                                                                        gCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGCTCGtt           <  3:183741/149‑1 (MQ=255)
                                                                                                                                         ccctGTCATCGAAGTTTTCGTTGACGGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 <  4:228693/110‑1 (MQ=38)
                                                                                                                                         ccctGTCATCGAAGTTTTCGTTGACGGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAagag                                                 >  3:228693/1‑110 (MQ=38)
                                                                                                                                                 aCGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGCTCGTTGAGGACGtt  <  3:330147/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTATTAGCAATTTGAGCGGTTATTCGGTTGGCATCCGCCAGCAGGCTCTGGCGCAGTTGCACGGTAAAGCCATGGAAGTGATC‑GCCGTGTGGAGCATTGGGCTCATGTTGGTCTGCCCTGGCCAATGCGACACTGGCCCCGTCAACGAAAACTTCGATGACAGGGCGCTGATCTGGCTGCGTGATATCGATAGCCCACCCCTGCAGAACATCCCCGTAAAGTCCGTTGATGCAGCCACGAAAAGAGGAGAGCGCTACAGGTTTCTGACCCGCCGGCCGCTCGTTGAGGACGTT  >  NC_002947/2011132‑2011424

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑