Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,142,186 T→C 8.8% S45S (AGT→AGC PP_2758 → ribose ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,142,1860TC 8.8% 121.1 / 2.2 45S45S (AGT→AGCPP_2758ribose ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/17);  new base C (2/2);  total (26/19)
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.21e-01

ATGAGTGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAA  >  NC_002947/3142052‑3142332
                                                                                                                                      |                                                                                                                                                  
aTGAGTGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTc                                                                                                                                      <  5:290710/149‑1 (MQ=255)
   aGTGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCt                                                                                                                                   >  4:290624/1‑149 (MQ=255)
   aGTGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCt                                                                                                                                   >  4:285884/1‑149 (MQ=255)
     tGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGc                                                                                                                                                  >  8:99700/1‑132 (MQ=255)
        cAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAgg                                                                                                                              <  7:22939/149‑1 (MQ=255)
         aaCGACTACTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAggg                                                                                                                             >  1:249281/1‑149 (MQ=255)
          aCGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAggg                                                                                                                             >  3:102034/1‑148 (MQ=255)
          aCGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAggg                                                                                                                             >  3:182810/1‑148 (MQ=255)
            gACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCaa                                                                                                                          <  1:192041/149‑1 (MQ=255)
               tGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAgcg                                                                                                                       >  8:190845/1‑149 (MQ=255)
               tGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAgcg                                                                                                                       >  8:185946/1‑149 (MQ=255)
                    cGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCg                                                                                                                  >  1:323207/1‑149 (MQ=255)
                        aTGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCAt                                                                                                              <  4:182810/149‑1 (MQ=255)
                        aTGTCCCGGCGCTGCCTGCTCGGCGCGGCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCAt                                                                                                              <  7:124771/149‑1 (MQ=255)
                                       ctgcTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGCGCCACAGTACCTTCg                                                                                                                                     >  4:128912/1‑111 (MQ=255)
                                       ctgcTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGCGCCACAGTACCTTCg                                                                                                                                     <  3:128912/111‑1 (MQ=255)
                                       ctgcTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGCGCCACAGTACCTTCg                                                                                                                                     <  3:124172/111‑1 (MQ=255)
                                       ctgcTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCCACAGCGGTGCCGAGGCCAGCGCCACAGTACCTTCg                                                                                                                                     >  4:124172/1‑111 (MQ=255)
                                        tgcTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCAc                                                                                               >  2:109277/1‑148 (MQ=255)
                                                gcgTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTa                                                                                      <  3:285884/149‑1 (MQ=255)
                                                gcgTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTa                                                                                      <  3:290624/149‑1 (MQ=255)
                                                  gTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACt                                                                                    <  4:102034/149‑1 (MQ=255)
                                                   tCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACtt                                                                                   <  1:109277/149‑1 (MQ=255)
                                                     cTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCAt                                                                                                              >  1:137347/1‑120 (MQ=255)
                                                     cTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCAt                                                                                                              <  2:137347/120‑1 (MQ=255)
                                                         gCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAg                                                                             >  3:66610/1‑149 (MQ=255)
                                                                tcgctcTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGcc                                                                      >  1:129455/1‑149 (MQ=255)
                                                                      tcGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGACGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAATGGCAAGCGCATCGCCAACAGCACCGTGGGCCCCCGCATCTACTTCGACAGACGGGCCTTCAAg                                                                >  6:191934/1‑149 (MQ=255)
                                                                           cTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCa                                                           <  2:323207/149‑1 (MQ=255)
                                                                              gCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGt                                                        >  8:87988/1‑149 (MQ=255)
                                                                                cGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTc                                                       >  2:211135/1‑148 (MQ=255)
                                                                                cGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTc                                                       >  4:188198/1‑148 (MQ=255)
                                                                                    ccATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAgg                                                   >  2:263909/1‑148 (MQ=255)
                                                                                              gcgcCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAgcgc                                        >  4:273715/1‑149 (MQ=255)
                                                                                                    tGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTGCGACAGCCGGGCCTGCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGggc                                  <  7:190845/149‑1 (MQ=255)
                                                                                                    gggcggcTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGggc                                  <  7:185946/148‑1 (MQ=255)
                                                                                                       cggcTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGggcggc                               <  7:87988/149‑1 (MQ=255)
                                                                                                             gAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGATCGAgg                                                   <  4:214583/123‑1 (MQ=255)
                                                                                                             gAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGATCGAgg                                                   >  3:214583/1‑123 (MQ=255)
                                                                                                                ccGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGCCGGCACGCCGATc                      >  3:30370/1‑149 (MQ=255)
                                                                                                                            ccGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCATGGATGcc          >  5:175293/1‑149 (MQ=255)
                                                                                                                                 gCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCCCCAGCATCTACTTCGCCAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGgcg     >  6:85291/1‑149 (MQ=255)
                                                                                                                                 gCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGTCACGCCGATCACCTTGGATGCCggg       <  5:71644/147‑1 (MQ=255)
                                                                                                                                 gCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGTCACGCCGATCACCTTGGATGCCggg       >  6:71644/1‑147 (MQ=255)
                                                                                                                                    aGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACATCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTTGGCGGCACGCCGATCACCTTGGATGCCGGGCGCaa  <  5:191934/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
ATGAGTGACAACGACTGCTCCGGCATGTCCCGGCGCTGCCTGCTCGGCGCGTCCTCGGCGATGCTCGCTCTCGGGCTTGCCGCACCATTGACCCGCGCCATGGCGGCTGGAGCCGACAGCGGTGCCGAGGCCAGTGCCACAGTACCTTCGCTGAAGGGCAAGCGCATCGCCATCAGCACCGTGGGCACCAGCATCTACTTCGACAGCCGGGCCTTCAAGGCGCAGGTCGAGGAAGTACAGCGCCTGGGCGGCACGCCGATCACCTTGGATGCCGGGCGCAA  >  NC_002947/3142052‑3142332

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: