Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,882,969 C→T 7.7% G75G (GGG→GGA PP_5156 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,882,9690CT 7.7% 152.3 / 2.1 52G75G (GGG→GGAPP_5156hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/24);  new base T (2/2);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GGTAATGCTGGTTTCGGTCTGTGCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTGGGCGGCTCCTGGCTCTT  >  NC_002947/5882822‑5883112
                                                                                                                                                   |                                                                                                                                                
ggTAATGCTGGTTTCGGTCTGTGCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGcc                                                                                                                                                 >  7:1221/1‑149 (MQ=255)
                gTCTGTGCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTg                                                                                                                                 <  3:206121/149‑1 (MQ=255)
                   tgtgCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGt                                                                                                                               >  1:37672/1‑148 (MQ=255)
                   tgtgCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGt                                                                                                                               >  4:115732/1‑148 (MQ=255)
                        cTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCg                                                                                                                         >  7:100277/1‑149 (MQ=255)
                         ttGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGttt                                                                                                                                  >  4:189504/1‑139 (MQ=255)
                            ccGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAgg                                                                                                                     >  2:146485/1‑149 (MQ=255)
                             cGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGt                                                                                                                    <  2:275411/149‑1 (MQ=255)
                               aaGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGttt                                                                                                                  >  6:72730/1‑149 (MQ=255)
                                 gTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGa                                                                                                                <  1:87470/149‑1 (MQ=255)
                                 gTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGa                                                                                                                <  7:58065/149‑1 (MQ=255)
                                    cagcaCTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGc                                                                                                             >  3:152189/1‑149 (MQ=255)
                                      gcaCTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGccc                                                                                                           <  8:100277/149‑1 (MQ=255)
                                            gCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATAccc                                                                                                     >  6:93782/1‑149 (MQ=255)
                                             cGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCgggg                                                                                                                                     >  6:73672/1‑116 (MQ=255)
                                             cGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCgggg                                                                                                                                     <  5:73672/116‑1 (MQ=255)
                                             cGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCt                                                                                                    <  5:136855/149‑1 (MQ=255)
                                             cGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCt                                                                                                    <  5:93782/149‑1 (MQ=255)
                                             cGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCAACAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCACATACCCt                                                                                                    <  6:184143/149‑1 (MQ=255)
                                                tCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTg                                                                                                                                 >  6:370167/1‑117 (MQ=255)
                                                tCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTg                                                                                                                                 <  5:370167/117‑1 (MQ=255)
                                                   cTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATAc                                                                                              >  4:319717/1‑149 (MQ=255)
                                                    ttCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATAcc                                                                                             >  5:278613/1‑149 (MQ=255)
                                                       ggCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGaa                                                                                          <  5:72730/149‑1 (MQ=255)
                                                       ggCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGGCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGaa                                                                                          <  4:263876/149‑1 (MQ=255)
                                                         cAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACt                                                                                        <  5:160194/149‑1 (MQ=255)
                                                              cGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAc                                                                                   <  4:152189/148‑1 (MQ=255)
                                                               gACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGTCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAc                                                                                   <  2:352582/147‑1 (MQ=255)
                                                               gACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGTCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAc                                                                                   <  2:347052/147‑1 (MQ=255)
                                                               gACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGTCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAc                                                                                   >  1:352582/1‑147 (MQ=255)
                                                               gACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGTCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAc                                                                                   >  1:347052/1‑147 (MQ=255)
                                                                  gCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTg                                                                               >  4:146997/1‑148 (MQ=255)
                                                                        ctctGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCt                                                                        <  3:115732/149‑1 (MQ=255)
                                                                            ggCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTc                                                                    >  5:158434/1‑149 (MQ=255)
                                                                              cATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGa                                                                  <  6:158434/149‑1 (MQ=255)
                                                                                    gCTGGTCGGCAGCGGCATGTCAGCCCCCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACACTCGATACCGAACATGACGGTGAAGGTCTTGTCGACCGTgg                                                            >  5:165973/1‑149 (MQ=255)
                                                                                        gTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCAc                                                         <  3:319717/148‑1 (MQ=255)
                                                                                                           ccACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGAt                                      >  8:287722/1‑148 (MQ=255)
                                                                                                            cACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGa                                    <  8:1221/149‑1 (MQ=255)
                                                                                                            cACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGa                                    <  1:146485/149‑1 (MQ=255)
                                                                                                               cGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGa                                 >  1:109886/1‑149 (MQ=255)
                                                                                                                     tttGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCg                           >  6:58771/1‑149 (MQ=255)
                                                                                                                           gACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGcctcct                     >  3:292132/1‑149 (MQ=255)
                                                                                                                           gACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGcctcc                      <  7:287722/148‑1 (MQ=255)
                                                                                                                             cATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTg                   >  4:127670/1‑149 (MQ=255)
                                                                                                                             cATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTg                   <  2:109886/149‑1 (MQ=255)
                                                                                                                               tCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTggg                 <  3:146997/149‑1 (MQ=255)
                                                                                                                                  accaGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTc                                                                         >  4:345115/1‑90 (MQ=255)
                                                                                                                                  accaGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTc                                                                         <  3:345115/90‑1 (MQ=255)
                                                                                                                                  accaGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAg                                  <  4:155898/129‑1 (MQ=255)
                                                                                                                                  accaGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAg                                  >  3:155898/1‑129 (MQ=255)
                                                                                                                                  accaGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTGggcgg              >  8:293820/1‑149 (MQ=255)
                                                                                                                                                gtgCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGAGAGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTGGGCGGCTCCTGGGTCtt  >  1:286004/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGTAATGCTGGTTTCGGTCTGTGCCTTGCCGAAGTGCAGCACTTGCGGTCCCTTCGGCAGTGCGACGCCGGCCTCTGGCATCAGGCTGGTCGGCAGCGGCATGTCAGCCACCGGCTCTTTGTCGACATCCACCAGCAGATGGTGGTGCCCGGTGTGCGGGGTTTGGTCACCCGCAGGTTTGAGGCCCATACCCTCGATACCGAACTTGACGGTGAAGGTCTTGTCGACCGTGGCACCATCGGCGGGGGAAACGATGAAGACCTTCGCCTCCTTGGGCGGCTCCTGGCTCTT  >  NC_002947/5882822‑5883112

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: