Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,010,130 G→A 8.1% A435V (GCG→GTG)  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,010,1300GA 8.1% 139.5 / 2.6 49A435V (GCG→GTG) PP_1792group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base A (2/2);  total (24/25)
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.86e-01

GGGCCAGGCCTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAA  >  NC_002947/2009987‑2010262
                                                                                                                                                |                                                                                                                                    
gggCCAGGCCTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGt                                                                                                                                 <  1:159860/149‑1 (MQ=255)
         cTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCACTAGTg                                                                                                                                <  4:167813/141‑1 (MQ=255)
         cTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCACTAGTg                                                                                                                                <  4:164385/141‑1 (MQ=255)
         cTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCACTAGTg                                                                                                                                >  3:167813/1‑141 (MQ=255)
         cTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCACTAGTg                                                                                                                                >  3:164385/1‑141 (MQ=255)
             gCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGa                                                                                                                    >  7:133512/1‑149 (MQ=255)
                     acgccacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCg                                                                                                            >  8:19679/1‑149 (MQ=255)
                     acgccacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGATAGTGCTATACCGGCGATTAGATCg                                                                                                            >  4:190132/1‑149 (MQ=255)
                         cacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTg                                                                                                                                >  1:47685/1‑125 (MQ=255)
                         cacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTg                                                                                                                                <  2:47685/125‑1 (MQ=255)
                         cacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGc                                                                                                        <  1:88072/149‑1 (MQ=255)
                         cacgGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGc                                                                                                        >  4:263932/1‑149 (MQ=255)
                            ggCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGt                                                                                                     <  3:263932/149‑1 (MQ=255)
                                            tttGAGGGAGCTGCTTGCTCCGGGCTTGTTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCt                                                                                      >  1:153569/1‑149 (MQ=255)
                                              tGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCtttt                                                                                   >  6:125494/1‑149 (MQ=255)
                                                aGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCttttt                                                                                  <  1:230460/148‑1 (MQ=255)
                                                aGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTg                                                                                 >  2:54663/1‑149 (MQ=255)
                                                  ggAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGcc                                                                               >  3:263096/1‑149 (MQ=255)
                                                    aGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTg                                                                             <  4:263096/149‑1 (MQ=255)
                                                       tgcttgctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCtt                                                                          >  5:215557/1‑149 (MQ=255)
                                                        gcttgctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCttt                                                                         >  5:130910/1‑149 (MQ=255)
                                                         cttgctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCttt                                                                         <  3:190132/148‑1 (MQ=255)
                                                         cttgctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCtt                                                                          <  2:89148/147‑1 (MQ=255)
                                                         cttgctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCtt                                                                          <  8:133512/147‑1 (MQ=255)
                                                            gctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCTGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAg                                                                      <  3:103183/148‑1 (MQ=255)
                                                            gctCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAg                                                                      <  6:215557/148‑1 (MQ=255)
                                                              tCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAg                                                                      <  5:48630/146‑1 (MQ=255)
                                                              tCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAg                                                                      >  6:48630/1‑146 (MQ=255)
                                                                cGGGCTTGTACAGGTATCGCTAGAGCGATGGTATAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAAt                                                                  <  4:149880/148‑1 (MQ=255)
                                                                 gggCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTg                                                                                 <  6:251476/132‑1 (MQ=255)
                                                                 gggCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTg                                                                                 >  5:251476/1‑132 (MQ=255)
                                                                   gCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCa                                                              <  6:225006/149‑1 (MQ=255)
                                                                       gTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCaaa                                                            <  7:19679/147‑1 (MQ=255)
                                                                          aCAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCct                                                        >  1:87742/1‑149 (MQ=255)
                                                                             accaTCGCTCTAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTg                                                                             <  2:133557/122‑1 (MQ=255)
                                                                             accaTCGCTCTAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTg                                                                             >  1:133557/4‑125 (MQ=255)
                                                                                     tAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCttt                                                                         <  2:88928/121‑1 (MQ=255)
                                                                                     tAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCttt                                                                         >  1:88928/1‑121 (MQ=255)
                                                                                       gagCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATggc                                           <  2:153569/149‑1 (MQ=255)
                                                                                         gCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATggcgg                                         <  5:125494/149‑1 (MQ=255)
                                                                                          cGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATggcgg                                         >  7:190731/1‑148 (MQ=255)
                                                                                                       cccGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTgcg                           <  6:130910/149‑1 (MQ=255)
                                                                                                        ccGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTtcgc                          >  7:191972/1‑149 (MQ=255)
                                                                                                                 gTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTagag                  >  7:102521/1‑148 (MQ=255)
                                                                                                                                gCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTAGAGAAAATGTCTGCGCCaa  >  7:98811/1‑149 (MQ=255)
                                                                                                                                     ttAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCATACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCttt                                                    <  7:110629/94‑1 (MQ=255)
                                                                                                                                     ttAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTATGAGCAATGCCAAACCCTCtta                                                    >  8:110629/1‑93 (MQ=255)
                                                                                                                                        gCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTAGAGAAAATGTCTgc       <  4:144661/136‑1 (MQ=255)
                                                                                                                                        gCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTAGAGAAAATGTCTgc       >  3:144661/1‑136 (MQ=255)
                                                                                                                                                |                                                                                                                                    
GGGCCAGGCCTTGGCATAACCACGCCACGGCGTCCTTGCGTTGCTTTGAGGGAGCTGCTTGCTCCGGGCTTGTACAGGTATCGCTAGAGCGATGGTAGAGAACCCGAGGCAGGTGCGTTACAGACCTGCCGTTTAGCTCCGTCGCTAGTGCTATACCGGCGATTAGATCGTGGCGGTTAATGGCAGACTGCTTTTTGCCTGCTTTGAGCAATGCCAAAGCCTCTTTGACAATGGCGGCAGAGAAAATTGCGCCAGTAGAGAAAATGTCTGCGCCAA  >  NC_002947/2009987‑2010262

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: