Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,007,598 G→A 8.5% G499G (GGC→GGT)  PP_2628 ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,007,5980GA 8.5% 132.2 / 2.1 47G499G (GGC→GGT) PP_2628ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (2/2);  total (22/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 = 5.12e-01

CGCCAGACAGCTTGATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCC  >  NC_002947/3007456‑3007727
                                                                                                                                               |                                                                                                                                 
cgccAGACAGCTTGATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAg                                                                                                                             >  5:303049/1‑149 (MQ=255)
      aCAGCTTGATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAata                                                                                                                       >  7:107260/1‑149 (MQ=255)
           ttGATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTc                                                                                                                  >  7:3579/1‑149 (MQ=255)
             gATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTc                                                                                                                  <  8:107260/147‑1 (MQ=255)
              aTGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTcgc                                                                                                                <  4:312322/148‑1 (MQ=255)
               tGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATa                                                                                                                              <  8:182232/133‑1 (MQ=255)
               tGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATa                                                                                                                              >  7:182232/1‑133 (MQ=255)
                gCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGaa                                                                                                             >  1:187067/1‑149 (MQ=255)
                   gcgTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATc                                                                                                           >  5:314278/1‑148 (MQ=255)
                   gcgTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATc                                                                                                           >  3:250983/1‑148 (MQ=255)
                        cacccaccAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACg                                                                                                      <  5:229352/148‑1 (MQ=255)
                         acccaccAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGt                                                                                                    >  8:267999/1‑149 (MQ=255)
                         acccaccAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGt                                                                                                    <  8:3579/149‑1 (MQ=255)
                             accaGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCaa                                                                                                >  5:333011/1‑149 (MQ=255)
                                    tgtCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAg                                                                                         <  6:303049/149‑1 (MQ=255)
                                    tgtCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAg                                                                                         <  7:267999/149‑1 (MQ=255)
                                     gtCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAgg                                                                                        >  5:258749/1‑149 (MQ=255)
                                         tAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTc                                                                                    >  6:101524/1‑149 (MQ=255)
                                             cGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGa                                                                                <  6:258749/149‑1 (MQ=255)
                                              gCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCt                                                                                  <  4:195964/146‑1 (MQ=255)
                                              gCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCt                                                                                  >  3:195964/1‑146 (MQ=255)
                                              gCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAg                                                                               <  2:187067/149‑1 (MQ=255)
                                                  ttGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGAc                                                                           >  3:8176/1‑149 (MQ=255)
                                                        tGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GACCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGc                                                                      >  3:77223/1‑148 (MQ=255)
                                                        tGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GACCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGc                                                                      >  1:63810/1‑148 (MQ=255)
                                                         ggCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGTCACTGTCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGt                                                                              <  3:315327/139‑1 (MQ=255)
                                                         ggCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGt                                                                              >  4:315327/1‑139 (MQ=255)
                                                          gCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GACCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGc                                                                      <  4:77223/146‑1 (MQ=255)
                                                          gCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GACCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGc                                                                      <  2:63810/146‑1 (MQ=255)
                                                           cGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGAt                                                                  <  6:314278/149‑1 (MQ=255)
                                                                      ctggctgATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCa                                                       <  5:101524/149‑1 (MQ=255)
                                                                         gctgATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAgg                                                     <  6:333011/148‑1 (MQ=255)
                                                                           tgatgaAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTCTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTa                                                  >  1:249255/1‑149 (MQ=255)
                                                                                aaaCCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCt                                              <  4:8176/148‑1 (MQ=255)
                                                                                   ccttCGGCCTGGGCACTGGCCTCAGCGGCACTGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGt                                          <  2:249255/146‑1 (MQ=255)
                                                                                    cATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGACCGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGt                                          >  5:344318/1‑148 (MQ=255)
                                                                                      tCGGCCTGGGCACTGGCCGCAGCGGGACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCTTAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGAGCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGAc                                       <  8:46133/149‑1 (MQ=255)
                                                                                      tCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAg                                                                                         <  4:47615/99‑1 (MQ=255)
                                                                                      tCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAg                                                                                         >  3:47615/1‑99 (MQ=255)
                                                                                      tCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATc                                                <  5:341727/140‑1 (MQ=255)
                                                                                      tCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATc                                                >  6:341727/1‑140 (MQ=255)
                                                                                                 aCTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGCGGGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGc                                            >  1:73585/1‑134 (MQ=255)
                                                                                                 aCTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGc                                            <  2:73585/133‑1 (MQ=255)
                                                                                                    ggCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCt                         >  3:52749/1‑149 (MQ=255)
                                                                                                               cACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATg              <  2:314467/149‑1 (MQ=255)
                                                                                                                aCGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGc             <  4:250983/149‑1 (MQ=255)
                                                                                                                           cTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTccccc  <  4:52749/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                 
CGCCAGACAGCTTGATGCCGCGTTCACCCACCAGGGTGTCGTAACCGCTGTTGCCCTGGCGGTCGCTGAGCTGGCTGATGAAACCATCGGCCTGGGCACTGGCCGCAGCGGCACGGATCTGTGCTTCGGTAGCATCGGGGC‑GGCCATAGGCAATATTGTCGCGAATCGAACGGTGCAACAGCGAGGTGTCCTGAGTGACCATGCCGATGGCGCTGCGCAGGCTATCCTGCGTGACCTTGGCGACGTTCTGCCCGTCGATGCGTATTTCCCCC  >  NC_002947/3007456‑3007727

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