Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,883,615 C→G 11.3% A197G (GCC→GGC) ‡ PP_3428 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,883,6150CG11.3% 155.1 / 11.4 53A197G (GCC→GGC) ‡PP_3428TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base G (3/3);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

CGCCATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGC  >  NC_002947/3883469‑3883761
                                                                                                                                                  |                                                                                                                                                  
cGCCATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCt                                                                                                                                                 <  1:904283/150‑1 (MQ=255)
       ggCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCCAGTCTGGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcga                                                                                                                                          <  1:798032/150‑1 (MQ=255)
       ggCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCagcga                                                                                                                                          >  1:249110/1‑150 (MQ=255)
           tctATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAgcgc                                                                                                                                      >  1:429984/1‑150 (MQ=255)
               tttCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTCCAGGGCTGCagcgg                                                                                                                                          >  1:453992/1‑141 (MQ=255)
               tttCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTCCAGGGCTGCagcgg                                                                                                                                          <  2:453992/142‑2 (MQ=255)
                       cccTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGGGCTGCagcggc                                                                                                                                         >  2:471679/1‑133 (MQ=255)
                       cccTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGGGCTGCagcggc                                                                                                                                         <  1:471679/135‑3 (MQ=255)
                               tCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAAc                                                                                                                  >  2:1186865/1‑150 (MQ=255)
                                  gCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTc                                                                                                               >  2:768589/1‑150 (MQ=255)
                                    ccTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGa                                                                                                             <  2:249110/150‑1 (MQ=255)
                                     cTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAt                                                                                                            <  2:429984/150‑1 (MQ=255)
                                              tACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGccag                                                                                                       >  1:1301251/1‑144 (MQ=255)
                                              tACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGccag                                                                                                       <  2:1301251/146‑3 (MQ=255)
                                                 cAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACt                                                                                                >  1:67311/1‑150 (MQ=255)
                                                  aGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTc                                                                                               >  1:1227274/1‑150 (MQ=255)
                                                  aGCAAGCCAGGCAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCCGCAGCGAGCGCTCGCGGTTGCCCCCGCCAACCTCGATGCCCTGGCGACTc                                                                                               <  2:454762/150‑1 (MQ=255)
                                                    cAAGCCCTGGAGTAACAGCCTGAAGATGCGGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAAAGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAGGCCCTGGCGACTCAg                                                                                             <  1:358838/150‑1 (MQ=255)
                                                       gCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCt                                                                                                       >  2:150488/1‑137 (MQ=255)
                                                       gCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCt                                                                                                       <  1:150488/137‑1 (MQ=255)
                                                          cTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGc                                                                                       <  1:494026/150‑1 (MQ=255)
                                                          cTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGc                                                                                       <  1:768589/150‑1 (MQ=255)
                                                          cTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGc                                                                                       >  2:494026/1‑150 (MQ=255)
                                                             gAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTg                                                                                    >  1:1173341/1‑150 (MQ=255)
                                                                cAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTc                                                                                 <  2:1227274/150‑1 (MQ=255)
                                                                  aCAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGa                                                                               >  2:1319200/1‑150 (MQ=255)
                                                                   cAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAg                                                                              >  2:1228235/1‑150 (MQ=255)
                                                                     gCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAgcg                                                                            >  2:639269/1‑150 (MQ=255)
                                                                       cTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGcc                                                                          >  2:686886/1‑150 (MQ=255)
                                                                         gAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACCCGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAg                                                                        <  2:1314848/150‑1 (MQ=255)
                                                                               gCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGccac                                                                  >  1:1290868/1‑150 (MQ=255)
                                                                                     gCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCacgac                                                            >  1:475034/1‑150 (MQ=255)
                                                                                       ccAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGc                                                          <  2:67311/150‑1 (MQ=255)
                                                                                               aCCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGcc                                                  <  2:1337395/150‑1 (MQ=255)
                                                                                                 cTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAg                                                >  1:565413/1‑150 (MQ=255)
                                                                                                               ccTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGACGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGc                                  >  1:643459/1‑150 (MQ=255)
                                                                                                                 tGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGccc                                >  2:371564/1‑150 (MQ=255)
                                                                                                                  gCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAt                                                                                                            <  2:878009/73‑1 (MQ=255)
                                                                                                                  gCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGAt                                                                                                            >  1:878009/1‑73 (MQ=255)
                                                                                                                    aCGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTgg                             <  1:1319200/150‑1 (MQ=255)
                                                                                                                      gAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGa                                                                               <  2:1219339/98‑1 (MQ=255)
                                                                                                                      gAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGa                                                                               >  1:1219339/1‑98 (MQ=255)
                                                                                                                      gAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGgcg                           >  2:911710/1‑150 (MQ=255)
                                                                                                                          aCGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGa                       <  1:795946/150‑1 (MQ=255)
                                                                                                                            gACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGc                     >  2:34743/1‑150 (MQ=255)
                                                                                                                            gACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGc                     >  1:922256/1‑150 (MQ=255)
                                                                                                                                  acgcacgcCGCTGCAGCCCTGCAGCGAGCGCTcga                                                                                                                                  >  1:858601/1‑34 (MQ=255)
                                                                                                                                  acgcacgcCGCTGCAGCCCTGCAGCGAGCGCTcga                                                                                                                                  <  2:858601/35‑2 (MQ=255)
                                                                                                                                    gcacgcCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTc             >  1:684531/1‑150 (MQ=255)
                                                                                                                                       gctcgcTGCAGGGCTGCAGCG‑GCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAg                                                                              <  2:11531/78‑1 (MQ=37)
                                                                                                                                       gctcgcTGCAGGGCTGCAGCG‑GCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAg                                                                              >  1:11531/4‑81 (MQ=37)
                                                                                                                                         ccgcTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCt        <  1:34743/150‑1 (MQ=255)
                                                                                                                                               cAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGc  <  2:684531/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCCATCGGCCTCTATTTCGGCCCCCTGCAGTCAGCCCTGGAGCAATACCAGCAAGCCCTGGAGCAACAGCCTGAAGATGCCGACGCCCAAGTGCACCTGGCGCGCCTGTACCTGCACGAAAACGACTACACGCACGCCGCTGCAGCCCTGCAGCGAGCGCTCGCGGTTGCCCCCGGCAACCTCGATGCCCTGGCGACTCAGGTACGCCTGCTCGAGCGCCAGGGCCACCACGACGCCTCCCGCCAGGTGCTGGGCAAAGCCCTGGCGCTGAGCCCCGACTCGGCCTTTCTGC  >  NC_002947/3883469‑3883761

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑