Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,143,968 C→T 7.4% R118H (CGC→CAC)  PP_4528 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,143,9680CT 7.4% 154.8 / 2.7 54R118H (CGC→CAC) PP_4528hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GCTCGGCAGCGTTCAGAAGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGAATCCGCGCTCGAATAGCAGCTTTTCA  >  NC_002947/5143820‑5144111
                                                                                                                                                    |                                                                                                                                               
gCTCGGCAGCGTTCAGAAGCTACCGATTATCGGGTTTTGTACGCTGCACTTGCCAGACTTAAAACGACCTCAAGCATTCTGTGCCTTCGCAGAAACCGGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGgc                                                                                                                                                 <  5:249976/149‑1 (MQ=255)
gCTCGGCAGCGTTCAGAAGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGgc                                                                                                                                                 <  6:201658/149‑1 (MQ=255)
           ttCAGAAGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAACCGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTc                                                                                                                                       >  3:119110/1‑148 (MQ=255)
            tCAGAAGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCg                                                                                                                                      >  4:12079/1‑148 (MQ=255)
                 aGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCGATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCg                                                                                                                                <  6:232947/149‑1 (MQ=255)
                   cTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACATCAAGCAATCTGTGACTTCGCAGAAACCCGTTCGAACGTCGTTTAGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTAGAAGCCggg                                                                                                                              >  1:293119/1‑149 (MQ=255)
                    tACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCggg                                                                                                                              >  6:304792/1‑148 (MQ=255)
                      ccGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCgggg                                                                                                                             <  4:275027/147‑1 (MQ=255)
                      ccGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGAt                                                                                                                           >  1:152296/1‑149 (MQ=255)
                       cGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGAtg                                                                                                                          <  7:216086/149‑1 (MQ=255)
                             tCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATgctgct                                                                                                                     <  5:304792/148‑1 (MQ=255)
                               gggTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTc                                                                                                                  >  6:140154/1‑149 (MQ=255)
                                     tGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATc                                                                                                              >  2:41403/1‑147 (MQ=255)
                                     tGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATc                                                                                                              <  1:41403/147‑1 (MQ=255)
                                       aaCGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGTGCgg                                                                                                                                             <  6:37714/114‑1 (MQ=255)
                                       aaCGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGTGCgg                                                                                                                                             <  6:42630/114‑1 (MQ=255)
                                       aaCGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGTGCgg                                                                                                                                             >  5:42630/1‑114 (MQ=255)
                                       aaCGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGTGCgg                                                                                                                                             >  5:37714/1‑114 (MQ=255)
                                            tGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCggg                                                                                                                              <  3:3490/124‑1 (MQ=255)
                                            tGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCggg                                                                                                                              >  4:3490/1‑124 (MQ=255)
                                                 ttGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAgcg                                                                                                >  8:148078/1‑149 (MQ=255)
                                                 ttGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAgcg                                                                                                >  8:222746/1‑149 (MQ=255)
                                                   gCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAgcgc                                                                                               >  7:140909/1‑148 (MQ=255)
                                                          ttAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCcgcg                                                                                       >  7:54311/1‑149 (MQ=255)
                                                          ttAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCcgcg                                                                                       >  7:291291/1‑149 (MQ=255)
                                                              aaCGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCg                                                                                   <  6:23993/149‑1 (MQ=255)
                                                                cGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGt                                                                                 <  8:40954/149‑1 (MQ=255)
                                                                    cTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAgg                                                                             <  2:152296/149‑1 (MQ=255)
                                                                    cTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAg                                                                              >  2:19160/1‑148 (MQ=255)
                                                                       aaGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGa                                                                            <  7:222746/147‑1 (MQ=255)
                                                                         gCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGaa                                                                        <  5:140154/149‑1 (MQ=255)
                                                                         gCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGaa                                                                        <  2:135247/149‑1 (MQ=255)
                                                                                gtgCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCg                                                                 >  5:147988/1‑149 (MQ=255)
                                                                                 tgCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGttt                                                                    <  7:148078/145‑1 (MQ=255)
                                                                                    cTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCAt                                                             <  8:140909/149‑1 (MQ=255)
                                                                                    cTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCAt                                                             >  4:33371/1‑149 (MQ=255)
                                                                                    cTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCAt                                                             >  2:257442/1‑149 (MQ=255)
                                                                                    cTTCGCAGAAACCCGTTCGAACGTCCTTTCGTCCAGCGCATCCTCCTGCTCATACAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATCCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCAt                                                             >  1:42039/1‑149 (MQ=255)
                                                                                     ttCGCAGAAACCCGTTCGAACGACGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCAtt                                                            >  3:209584/1‑149 (MQ=255)
                                                                                            aaaCCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGccgccg                                                     >  2:217853/1‑149 (MQ=255)
                                                                                             aaCCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGt                                                    <  6:147988/149‑1 (MQ=255)
                                                                                             aaCCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGt                                                    <  3:12079/149‑1 (MQ=255)
                                                                                                   ttCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCtt                                              <  8:54311/149‑1 (MQ=255)
                                                                                                      gAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCa                                            <  3:33371/148‑1 (MQ=255)
                                                                                                      gAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCa                                            <  8:291291/148‑1 (MQ=255)
                                                                                                         cgtcgtTTCGTCCATCGAATCCTGCTTCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGTGATGCTGCTGTCGTTCAGGCTCATCAAGCGGGAGCCGCGTTCGTTCTGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTc                                        <  2:293119/149‑1 (MQ=255)
                                                                                                          gtcgtTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACg                                                                          <  6:224923/114‑1 (MQ=255)
                                                                                                          gtcgtTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACg                                                                          >  5:224923/1‑114 (MQ=255)
                                                                                                                     cAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGa                            <  1:127063/149‑1 (MQ=255)
                                                                                                                               tcctGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAg                                                                              >  7:220628/1‑89 (MQ=255)
                                                                                                                               tcctGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAg                                                                              <  8:220628/89‑1 (MQ=255)
                                                                                                                               tcctGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGAATCCGCGCTc                  >  5:59236/1‑149 (MQ=255)
                                                                                                                                      cATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTgcc                                                         <  6:284810/103‑1 (MQ=255)
                                                                                                                                      cATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTgcc                                                         >  5:284810/1‑103 (MQ=255)
                                                                                                                                               ccTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGAATCCGCGCTCGAATAGCAGCTTTTCa  <  1:257442/149‑1 (MQ=255)
                                                                                                                                                cTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGAATCCGCGCTCGAATAGCAGCTTTTCa  >  4:25778/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCTCGGCAGCGTTCAGAAGCTACCGATTATCGGGTTTTGAACGCTGCACTTGCCAGCCTTAAAACGACCTCAAGCAATCTGTGCCTTCGCAGAAACCCGTTCGAACGTCGTTTCGTCCAGCGCATCCTCCTGCTCATCCAGCACCTGGCGCGGGTGCTCGAAGCCGGGGATGCTGCTGTCGATCAGGCTCATCAAGCGCGAGCCGCGTTCGGTCAGGACGAAGTTTTCGCCATTGCCGCCGTCCTCTTCAGGTCGGCATTCAATGAATCCGCGCTCGAATAGCAGCTTTTCA  >  NC_002947/5143820‑5144111

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