Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,944,185 G→A 8.5% A17V (GCC→GTC)  PP_5211 ← ChaC‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,944,1850GA 8.5% 132.8 / 2.6 47A17V (GCC→GTC) PP_5211ChaC‑related protein
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

GCCAAACACCAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGACATACTTTCAAGTGCCGACATAAG  >  NC_002947/5944046‑5944330
                                                                                                                                           |                                                                                                                                                   
gCCAAACACCAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGcc                                                                                                                                            <  7:11139/149‑1 (MQ=255)
 ccAAACACCAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCca                                                                                                                                           >  1:203435/1‑149 (MQ=255)
       acCAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAgg                                                                                                                                     >  6:212576/1‑149 (MQ=255)
         cAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCACGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGgcc                                                                                                                                   >  8:29811/1‑149 (MQ=255)
              cGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGCGGCATAGGCAAACAGCCACACAGgcccgccc                                                                                                                              >  6:57732/1‑149 (MQ=255)
               gggCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCt                                                                                                                             >  8:125769/1‑149 (MQ=255)
                             tACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATg                                                                                                              <  4:149105/149‑1 (MQ=255)
                             tACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATg                                                                                                              >  3:195654/1‑149 (MQ=255)
                                   gATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTc                                                                                                        <  7:235533/149‑1 (MQ=255)
                                   gATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTc                                                                                                        >  3:57875/1‑149 (MQ=255)
                                      gTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCtt                                                                                                      >  2:206791/1‑148 (MQ=255)
                                          aTGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATg                                                                                                  >  4:125467/1‑148 (MQ=255)
                                            gggACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAg                                                                                               >  2:16345/1‑149 (MQ=255)
                                             ggACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGTGTGCCTGGACATGGTTGTCTTCATGGAg                                                                                               <  1:16345/149‑1 (MQ=255)
                                                  acaGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGaac                                                                                        <  6:47471/149‑1 (MQ=255)
                                                  acaGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGaac                                                                                        <  8:48152/149‑1 (MQ=255)
                                                   caGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGaac                                                                                        >  6:38439/1‑148 (MQ=255)
                                                     ggTACAAGCCACGGTGATAGCCATGTACCCACGAGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGaacaac                                                                                     >  6:250684/1‑149 (MQ=255)
                                                             ccACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTcgc                                                                              >  4:172865/1‑148 (MQ=255)
                                                              cACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGa                                                                            <  7:125769/149‑1 (MQ=255)
                                                              cACGGTGATAGCCATGTACCCACGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGa                                                                            <  7:29811/149‑1 (MQ=255)
                                                               aCGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAg                                                                           >  8:54631/1‑149 (MQ=255)
                                                                      tAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTg                                                                    >  1:169505/1‑149 (MQ=255)
                                                                                     gcgcgcTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAg                                                     <  5:212576/149‑1 (MQ=255)
                                                                                     gcgcgcTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAg                                                     <  5:38439/149‑1 (MQ=255)
                                                                                           tGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTcgc                                                                              <  4:158683/118‑1 (MQ=255)
                                                                                           tGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTcgc                                                                              >  3:158683/1‑118 (MQ=255)
                                                                                              cgcTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAgg                                            <  2:203435/149‑1 (MQ=255)
                                                                                                cTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGCGCTGAACAACTGTTCGCGAGTATGTTTCTGACCGAAATCAAGTGCCGGAGGGt                                          <  1:206791/149‑1 (MQ=255)
                                                                                                cTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGt                                          <  2:169505/149‑1 (MQ=255)
                                                                                                   cACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTAt                                        >  5:208309/1‑148 (MQ=255)
                                                                                                   cACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGACAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTAt                                        >  5:248746/1‑148 (MQ=255)
                                                                                                   cACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGACAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTAt                                        >  5:252685/1‑148 (MQ=255)
                                                                                                        aaTTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACt                                  >  2:33721/1‑149 (MQ=255)
                                                                                                          ttGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTc                                <  4:195654/149‑1 (MQ=255)
                                                                                                            gCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTGTGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTcc                               <  3:172865/148‑1 (MQ=255)
                                                                                                              aCTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCaa                            <  7:54631/149‑1 (MQ=255)
                                                                                                               cTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGa                                              >  5:136939/1‑130 (MQ=255)
                                                                                                               cTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGa                                              <  6:136939/130‑1 (MQ=255)
                                                                                                                  cGGGCGCCAGATCAACGAGGCATAGACAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGa                          <  6:252685/147‑1 (MQ=255)
                                                                                                                  cGGGCGCCAGATCAACGAGGCATAGACAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGa                          <  6:248746/147‑1 (MQ=255)
                                                                                                                               aaCGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGACATACTTTCAAGTGc           <  1:33721/149‑1 (MQ=255)
                                                                                                                                       aTAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGACATACTTTCAAGTGCCGACAt     >  7:177762/1‑147 (MQ=255)
                                                                                                                                        tAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGt                                                    >  1:209624/1‑100 (MQ=255)
                                                                                                                                        tAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGa                          <  2:164988/125‑1 (MQ=255)
                                                                                                                                        tAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGa                          >  1:164988/1‑125 (MQ=255)
                                                                                                                                        tAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGACATACTTTCAAGTGCCGACATAAg  <  3:125467/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GCCAAACACCAGGCCGGGCGTTTCCGGGGTACCGCGATGTTCATGGGACCACAGGTACAAGCCACGGTGATAGCCATGTACCCGCGCGCGCTGGCGCTCCACCGAATTGCACTCCGGGCGCCAGATCAACGAGGCATAGGCAAACAGCCACACAGGCCCGCCCTGGTGCCTGGACATGGTTGTCTTCATGGAGCTGAACAACTGTTCGCGAGTATGTTGCTGACCGAAATCAAGTGCCGGAGGGTATGAAACTTCCCAAGACATACTTTCAAGTGCCGACATAAG  >  NC_002947/5944046‑5944330

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: