Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 36,886 T→C 10.5% Y256C (TAT→TGT)  PP_0033 ← undecaprenyl‑glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294736,8860TC10.5% 96.6 / 2.7 38Y256C (TAT→TGT) PP_0033undecaprenyl‑glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (18/16);  new base C (2/2);  total (20/18)
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

CCACAGCCAGAGGTAGGCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATCCAGG  >  NC_002947/36747‑37033
                                                                                                                                           |                                                                                                                                                   
ccACAGCCAGAGGTAGGCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTccc                                                                                                                                            <  5:340576/149‑1 (MQ=255)
      ccAGAGGTAGGCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGccc                                                                                                                                      <  6:24504/149‑1 (MQ=255)
                gCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTgac                                                                                                                            <  4:188583/149‑1 (MQ=255)
                gCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGAAGCGCAAATCCGAGCTGGATCCAGATGCTGCCGTCCGGCTTGGCATGAGCACTACCATACAGGTCCCCTGCCCAGTTACTgac                                                                                                                            >  7:281347/1‑149 (MQ=255)
                  tttGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATag                                                                                                                                                   <  6:76519/124‑1 (MQ=255)
                  tttGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATag                                                                                                                                                   >  5:76519/1‑124 (MQ=255)
                       aGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTa                                                                                                                      <  6:340947/148‑1 (MQ=255)
                               aTCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCa                                                                                                             <  8:49164/149‑1 (MQ=255)
                                      cACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTc                                                                                                       <  4:46142/148‑1 (MQ=255)
                                                       aGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCACAGAGGTCCCCTGCCCAGTTACTg                                                                                                                              >  3:105761/1‑108 (MQ=255)
                                                       aGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCACAGAGGTCCCCTGCCCAGTTACTg                                                                                                                              <  4:105761/108‑1 (MQ=255)
                                                        gTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTc                                                                                     >  2:124483/1‑148 (MQ=255)
                                                         tGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCa                                                                                   >  2:37620/1‑149 (MQ=255)
                                                            tAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATg                                                                                 >  5:64603/1‑148 (MQ=255)
                                                            tAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATg                                                                                 >  1:97616/1‑148 (MQ=255)
                                                            tAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCACAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATg                                                                                 >  8:32811/1‑148 (MQ=255)
                                                                   aaGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCggg                                                                         <  6:64603/149‑1 (MQ=255)
                                                                          aTGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCCCCAACGAAGAAGTACTCCATGAAGCCGGGCGtttt                                                                   >  8:204813/1‑148 (MQ=255)
                                                                                gTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCg                                                            >  4:212173/1‑149 (MQ=255)
                                                                                           ccGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATa                                                 <  3:212173/149‑1 (MQ=255)
                                                                                             gAGCTGGATCCAGATGCTGCGGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATAc                                                >  5:240929/1‑148 (MQ=255)
                                                                                                      ccAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATAACATGTCaag                                       >  8:229272/1‑146 (MQ=255)
                                                                                                       cAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGa                                     <  7:229272/149‑1 (MQ=255)
                                                                                                        aGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGc                                                                             >  4:196033/1‑108 (MQ=255)
                                                                                                           tgctgcCGTACGGCTTGGCATGAGCACTACCACAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGAtata                                 <  7:32811/149‑1 (MQ=255)
                                                                                                                ccGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGaa                            <  1:124483/149‑1 (MQ=255)
                                                                                                                  gTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACa                          >  8:32671/1‑149 (MQ=255)
                                                                                                                     cGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGat                                    <  2:196050/136‑1 (MQ=255)
                                                                                                                     cGGCTTGGCATGAGCACTACCATAGAGGACCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGat                                    >  1:196050/1‑136 (MQ=255)
                                                                                                                           ggCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACgaag                                                                                             <  6:39083/73‑1 (MQ=255)
                                                                                                                           ggCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACgaag                                                                                             >  5:39083/1‑73 (MQ=255)
                                                                                                                           ggCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTAc                   <  8:204419/147‑1 (MQ=255)
                                                                                                                                gAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCg            >  7:45574/1‑149 (MQ=255)
                                                                                                                                 aGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGc           <  7:32671/149‑1 (MQ=255)
                                                                                                                                      tACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATc      >  3:38253/1‑149 (MQ=255)
                                                                                                                                        ccATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATc      >  2:133296/1‑147 (MQ=255)
                                                                                                                                        ccATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATc      <  8:45574/147‑1 (MQ=255)
                                                                                                                                          aTAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATCCAgg  >  1:212386/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
CCACAGCCAGAGGTAGGCTTTGAAGCGCTGCATCGATCCACGCTTGCGGATCAACAGTGGTAGAAACAAGGCCCATGGCAGTAGCGACAATCCGAGCTGGATCCAGATGCTGCCGTACGGCTTGGCATGAGCACTACCATAGAGGTCCCCTGCCCAGTTACTGACGACGTAGCGGCTCCAGTGCTCACCAACGAAGAAGTACTCCATGAAGCCGGGCGTTTTCATCTCGGCCGCGACATACCATGGCACGGATATACCGAACATGAGTACAAGGCCGCTAATCCAGG  >  NC_002947/36747‑37033

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: