Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,022 C→A 7.4% R655L (CGT→CTT)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,0220CA 7.4% 171.5 / 3.3 54R655L (CGT→CTT) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/28);  new base A (2/2);  total (24/30)
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.65e-01

AAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGC  >  NC_002947/5538879‑5539166
                                                                                                                                               |                                                                                                                                                
aaGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCg                                                                                                                                              <  3:38697/148‑1 (MQ=255)
  gTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTg                                                                                                                                           >  6:41691/1‑149 (MQ=255)
  gTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTg                                                                                                                                           <  4:129576/149‑1 (MQ=255)
  gTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGACCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAACATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCGTCTGCACAAGGCGttc                                                                                                                                           >  1:130320/1‑148 (MQ=255)
    aTCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCt                                                                                                                                         >  4:181167/1‑149 (MQ=255)
     tCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCtc                                                                                                                                        >  8:175941/1‑149 (MQ=255)
     tCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCtc                                                                                                                                        <  5:41691/149‑1 (MQ=255)
     tCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGACCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAACATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCGTCTGCACAAGGCGttctt                                                                                                                                         <  2:130320/148‑4 (MQ=255)
        gTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGGCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCAGCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCtcgtc                                                                                                                                     <  2:180375/149‑1 (MQ=255)
             cGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGctc                                                                                                                                 <  2:33899/148‑1 (MQ=255)
                  ggTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACg                                                                                                                           <  5:147562/149‑1 (MQ=255)
                   gTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACg                                                                                                                           <  8:143616/148‑1 (MQ=255)
                      cAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGTCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtct                                                                                                                       >  6:284720/1‑149 (MQ=255)
                       aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtct                                                                                                                       <  2:210883/148‑1 (MQ=255)
                       aGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGAtct                                                                                                                       >  1:210883/1‑148 (MQ=255)
                                 gTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCATTATATTACAGACCTTCGGGTCGTGGTATCTCAATAAATACTTGTGAAGACTTCGTTTCTCGTCGCTCTTGACGATCTCCACCCCTTCg                                                                                                            <  2:234524/149‑1 (MQ=255)
                                     gTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGctc                                                                                                         <  3:181167/148‑1 (MQ=255)
                                        cTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCtt                                                                                                       <  1:16024/147‑1 (MQ=255)
                                        cTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCtt                                                                                                       >  2:16024/1‑147 (MQ=255)
                                           gcccggcccgATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAg                                                                                                   <  7:116373/148‑1 (MQ=255)
                                                          gTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAc                                                                                              <  2:239462/138‑1 (MQ=255)
                                                          gTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCAc                                                                                              >  1:239462/1‑138 (MQ=255)
                                                           tcggctcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCgg                                                                                  >  1:276584/1‑149 (MQ=255)
                                                                tcggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCt                                                                             >  5:187879/1‑149 (MQ=255)
                                                                  ggtccATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCtt                                                                            <  5:284720/145‑1 (MQ=255)
                                                                  ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGTGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCtt                                                                            <  4:272456/148‑1 (MQ=255)
                                                                  ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCt                                                                                                               >  3:146783/1‑113 (MQ=255)
                                                                  ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCt                                                                                                               <  4:146783/113‑1 (MQ=255)
                                                                  ggcCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCtt                                                                            <  7:122970/148‑1 (MQ=255)
                                                                     ccATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCCCCCCTTCGCTCTTGTAGGTCACCTTGCGCCGCGGGTTCTTGCCt                                                                        >  6:76771/1‑149 (MQ=255)
                                                                     ccATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                        >  5:212479/1‑149 (MQ=255)
                                                                      cATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                        >  8:188820/1‑148 (MQ=255)
                                                                       attcGCTGCAGCGCTTCACGCAGCATCGGCCAGCCGTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCTCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGa                                                                      <  6:80367/146‑1 (MQ=255)
                                                                       aTACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCt                                                                        >  7:102321/1‑147 (MQ=255)
                                                                         aCGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGt                                                                    <  8:102321/149‑1 (MQ=255)
                                                                          cGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTg                                                                   <  7:175941/149‑1 (MQ=255)
                                                                                    tgaagcGCTGCAGCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggcc                                                         <  6:47091/144‑1 (MQ=255)
                                                                                          gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTc                                                                              >  7:195958/1‑122 (MQ=255)
                                                                                          gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTc                                                                              <  8:195958/122‑1 (MQ=255)
                                                                                          gcagcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACAtggccgtggcc                                                   <  1:73398/149‑1 (MQ=255)
                                                                                            agcaTCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGC‑GGTTCTTTCCTGAGTGGTAAAaggccgtggccgcg                                                >  7:28882/1‑147 (MQ=255)
                                                                                                tCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGc                                              >  3:174528/1‑148 (MQ=255)
                                                                                                                 tCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGaa                            >  8:41444/1‑149 (MQ=255)
                                                                                                                       tAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTg                      >  2:23972/1‑149 (MQ=255)
                                                                                                                            gagcaaCGCCGTCTGCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCt                 >  5:237784/5‑149 (MQ=255)
                                                                                                                                aaCGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTc             <  7:188820/149‑1 (MQ=255)
                                                                                                                                 aCGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCg            <  1:23972/149‑1 (MQ=255)
                                                                                                                                 aCGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCg            <  4:174528/149‑1 (MQ=255)
                                                                                                                                    ccTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGgcg         >  4:200923/1‑149 (MQ=255)
                                                                                                                                       tctgCACAAGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCgg      <  6:237784/147‑1 (MQ=255)
                                                                                                                                         tgCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCAAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGaa    <  6:212479/149‑1 (MQ=255)
                                                                                                                                           cAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGc  >  7:121262/1‑149 (MQ=255)
                                                                                                                                           cAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGc  <  2:249834/149‑1 (MQ=255)
                                                                                                                                           cAGACGGCGTTGCTCGTCGCTCTTGACAATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGc  >  1:170479/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
AAGTATCGGTCTGCGGCTGGTGCAGCGGGATCAGTTGGTGCTTGCCCGGCCCGATCAGGTCGGCTCGGCCCATACGCTGCAGCGCTTCACGCAGCATCGGCCAGCCCTTCGGGTCGTGGTAGCGCAAGAACGCCTTGTGCAGACGGCGTTGCTCGTCGCTCTTGACGATCTCCACCCCTTCGCTCTTGTAGGTCACCTTGCGCAGCGGGTTCTTGCCTGAGTGGTACATGGCCGTGGCCGAGGCCATTGGCGACGGGTAGAACGCCTGTACCTGGTCGGCGCGGAAGC  >  NC_002947/5538879‑5539166

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