Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 5,641,230 (C)7→8 coding (349/1320 nt) PP_4950 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,641,2231.C90.9% 186.0 / 11.2 55coding (342/1320 nt)PP_4950hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (20/30);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGCGACGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑CCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGGCGCCTTCGGCTGG  >  NC_002947/5641092‑5641355
                                                                                                                                     |                                                                                                                                      
ggcgACGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑cccc                                                                                                                                   >  1:661231/1‑136 (MQ=255)
    aCGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑cccccccc                                                                                                                              <  1:1093265/136‑2 (MQ=255)
    aCGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑cccccccc                                                                                                                              >  1:570568/1‑135 (MQ=255)
    aCGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑cccccccc                                                                                                                              >  1:646089/1‑135 (MQ=255)
         cAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGccc                                                                                                                          <  1:485110/136‑1 (MQ=255)
          aaCTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCg                                                                                                                         >  2:1176043/1‑136 (MQ=255)
          aaCTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCg                                                                                                                         >  1:721341/1‑136 (MQ=255)
              ggATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACt                                                                                                                     <  2:721341/136‑1 (MQ=255)
               gATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGTCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑CCCCCCCTGCCCGCACTGa                                                                                                                   <  2:541290/136‑1 (MQ=255)
               gATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCATACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTg                                                                                                                    <  1:679770/136‑1 (MQ=255)
                aTCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGa                                                                                                                   <  1:158288/136‑1 (MQ=255)
                aTCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGa                                                                                                                   <  1:29264/136‑1 (MQ=255)
                 tCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGaa                                                                                                                  <  1:1059732/136‑1 (MQ=255)
                  cGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAAc                                                                                                                 >  1:411308/1‑136 (MQ=255)
                     cAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGaa                                                                                                              <  1:1014625/136‑1 (MQ=255)
                     cAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGaa                                                                                                              <  1:979542/136‑1 (MQ=255)
                     cAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGaa                                                                                                              <  1:1145694/136‑1 (MQ=255)
                        cgcCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑CCCCCCCTGCCCGCACTGAACGAAGTAc                                                                                                          <  1:605742/136‑1 (MQ=255)
                         gcCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTAc                                                                                                          <  1:258673/136‑1 (MQ=255)
                         gcCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTAc                                                                                                          <  1:510508/136‑1 (MQ=255)
                           cTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTg                                                                                                        <  2:863064/136‑1 (MQ=255)
                                cGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCt                                                                                                   <  1:659297/136‑1 (MQ=255)
                                 gTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTg                                                                                                  <  1:382673/136‑1 (MQ=255)
                                 cgtGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTg                                                                                                  <  1:218079/134‑1 (MQ=255)
                                      gcTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCgc                                                                                             >  2:529401/1‑135 (MQ=255)
                                      gcTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACGGACCGAAGTACTGGCCCTGCTCgg                                                                                             >  2:504423/1‑135 (MQ=255)
                                        tCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGc                                                                                           <  1:705275/136‑1 (MQ=255)
                                             cagcTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGaa                                                                                      >  1:650846/1‑136 (MQ=255)
                                                      cAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTgac                                                                             <  1:335873/136‑1 (MQ=255)
                                                         aCCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTgacctg                                                                          <  1:630366/136‑1 (MQ=255)
                                                         aCCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTgacctg                                                                          <  1:1276470/136‑1 (MQ=255)
                                                               ccGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTgacctggacctg                                                                    <  2:322858/136‑1 (MQ=255)
                                                               ccGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTgacctggacctg                                                                    <  2:650846/136‑1 (MQ=255)
                                                                cGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCGCGAGCGTGAAGCAGGTgacctggacctg                                                                    >  2:60317/1‑136 (MQ=255)
                                                                cGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGa                                                                                                                   >  2:1353034/1‑88 (MQ=255)
                                                                cGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGa                                                                                                                   <  1:1353034/88‑1 (MQ=255)
                                                                cGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGGGAAGCAGGTGacctggacctgg                                                                   >  2:502918/1‑136 (MQ=255)
                                                                       ggCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCAt                                                            <  1:875298/136‑1 (MQ=255)
                                                                       ggCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGGGACCTGGACCTGGTCGGCAt                                                            <  1:956350/136‑1 (MQ=255)
                                                                                  cTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGATC‑‑CCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGccgccg                                                 >  1:636496/1‑136 (MQ=255)
                                                                                        tGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTg                                                                                                        >  2:106975/1‑75 (MQ=255)
                                                                                        tGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTg                                                                                                        <  1:106975/75‑1 (MQ=255)
                                                                                             cTGGACGAAAACGCCTGGCACAGCCTTAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGc                                      >  2:844461/1‑136 (MQ=255)
                                                                                             cTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTTCTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGTCCTGGTCGGCATCTACGCCGCCGGCCCGGTATGc                                      >  1:270952/1‑136 (MQ=255)
                                                                                             cTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGCCCTGGTCGGCCTCTACGCCGCCGGCCCGGTATGc                                      >  1:1047859/1‑136 (MQ=255)
                                                                                             cTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGc                                      >  2:719371/1‑136 (MQ=255)
                                                                                             cTGGAAGAAAACGCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCACGCACTGAACGAAGTACTGCCCAAGCACGAGCCAGAACCAGCTGACCAGGACCAGGTCGGCATCAACACCGCCGGCCCCGTATGc                                      >  2:541094/1‑136 (MQ=255)
                                                                                                         gCCTGGCACAGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGcc                          <  2:1145747/136‑1 (MQ=255)
                                                                                                                acaGCCATAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGt                   >  1:438175/1‑136 (MQ=255)
                                                                                                                      aTAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcg             <  2:122512/136‑1 (MQ=255)
                                                                                                                      aTAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcg             <  2:646089/136‑1 (MQ=255)
                                                                                                                      aTAGCCTGCAGGAA‑‑CCACCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTAGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCTGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcgc            >  1:871181/1‑136 (MQ=255)
                                                                                                                       tAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcgc            >  1:1120327/1‑136 (MQ=255)
                                                                                                                       tAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcgc            >  1:1233747/1‑136 (MQ=255)
                                                                                                                       tAGCCTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGgcgc            <  2:626985/136‑1 (MQ=255)
                                                                                                                           cTGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGGCGCCTTc        <  2:770970/136‑1 (MQ=255)
                                                                                                                            tGCAGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGGCGCCTTCg       >  2:452745/1‑136 (MQ=255)
                                                                                                                               aGGAACCCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCGTCGCCCGCTCGTTCGGCGCCTTCGGcg    >  1:616132/1‑135 (MQ=255)
                                                                                                                                 gAA‑‑CCCCCCACTGCCCGCACTGAACGAAGTACTGGCACTGCTCGAGCGTGAAGCAGGTGACCTGGACATGGACGGCAACTACGCCGCCGGCCCGGTATGCCGCGGCATCGCCAGCTCCTTCGGCGCCTACGGCTgg  >  2:146853/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                      
GGCGACGTGCAACTGGATCGCCAGCGCCTGAACGTTGCGCTCGAGCAGCTGCGCCAGACCTTGCCGCTGCTGGCGCTCGACCCTTACCTGCGCCTGGACGAAAACGCCTGGCACAGCCATAGCCTGCAGGAA‑‑CCCCCCCTGCCCGCACTGAACGAAGTACTGGCCCTGCTCGAGCGTGAAGCAGGTGACCTGGACCTGGTCGGCATCTACGCCGCCGGCCCGGTATGCCGCGGCTTCGCCAGCTCGTTCGGCGCCTTCGGCTGG  >  NC_002947/5641092‑5641355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: