Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 675,591 A→T 11.6% A36A (GCT→GCA)  PP_0577 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947675,5910AT11.6% 103.2 / 3.9 41A36A (GCT→GCA) PP_0577hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base T (3/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCA  >  NC_002947/675444‑675737
                                                                                                                                                   |                                                                                                                                                  
gCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGcagcca                                                                                                                                                    <  1:121153/148‑1 (MQ=255)
 ccGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGcagccagc                                                                                                                                                  <  6:44754/149‑1 (MQ=255)
         ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAgccagcca                                                                                                                                                <  4:261824/143‑1 (MQ=255)
         ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAgccagcca                                                                                                                                                >  3:261824/1‑143 (MQ=255)
         ccagctccagctGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCaa                                                                                                                                          <  1:87234/149‑1 (MQ=255)
                    tGCGGTAGTGCTCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCGTGGGCGGGGGCCGCTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCACCCAGGGCAATGGGCAGCAgg                                                                                                                               >  8:270866/1‑149 (MQ=255)
                            tGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGca                                                                                                                       >  1:319274/1‑149 (MQ=255)
                              cccATCATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCagccagcc                                                                                                                                                 <  5:317032/121‑1 (MQ=255)
                              cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCagccagcc                                                                                                                                                 >  6:317032/1‑121 (MQ=255)
                              cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCACGcac                                                                                                                      >  8:201410/1‑148 (MQ=255)
                              cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGACAATGGGCCGCAGGGCCAGGc                                                                                                                        >  2:282577/1‑146 (MQ=255)
                              cccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGACAATGGGCAGCAGGGCCAGGc                                                                                                                        <  1:282577/146‑1 (MQ=255)
                               ccATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCacaa                                                                                                                    <  2:204644/149‑1 (MQ=255)
                                           gCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACgg                                                                                                        >  2:43043/1‑149 (MQ=255)
                                                 gcAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACg                                                                                                              >  6:64923/1‑137 (MQ=255)
                                                 gcAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACg                                                                                                              <  5:64923/137‑1 (MQ=255)
                                                    gggTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTAttt                                                                                               <  7:267283/149‑1 (MQ=255)
                                                                   tCGAAGGTGGGGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGcggcgg                                                                                <  2:319274/149‑1 (MQ=255)
                                                                         ccggcccccgccCAGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGAt                                                                              >  5:255694/8‑145 (MQ=255)
                                                                         ccggcccccgccCAGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGAt                                                                              <  6:255694/138‑1 (MQ=255)
                                                                            gtggCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCg                                                                         <  5:295791/147‑1 (MQ=255)
                                                                            gtggCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCg                                                                         >  6:295791/1‑147 (MQ=255)
                                                                               gccgccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                           >  3:182268/1‑126 (MQ=255)
                                                                               gccgccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                           <  4:182268/126‑1 (MQ=255)
                                                                                      gggCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCgggg                                                                      <  7:8683/140‑1 (MQ=255)
                                                                                      gggCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCgggg                                                                      >  8:8683/1‑140 (MQ=255)
                                                                                             ggCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGc                                                      <  1:43043/149‑1 (MQ=255)
                                                                                                cGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCACACCCAGTGCAgg                                                   >  5:131308/1‑149 (MQ=255)
                                                                                                         cGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCa                                          <  6:131308/149‑1 (MQ=255)
                                                                                                                 tgctgGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCAGGc                                  >  7:146682/1‑149 (MQ=255)
                                                                                                                    tgGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGaa                               <  8:146682/149‑1 (MQ=255)
                                                                                                                      gCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGGATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAgg                             >  7:204425/1‑149 (MQ=255)
                                                                                                                             aCGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtggggtg                      >  5:222394/1‑149 (MQ=255)
                                                                                                                                    gctgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACgg               >  8:291468/1‑149 (MQ=255)
                                                                                                                                     ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGAc                 <  3:218099/145‑1 (MQ=255)
                                                                                                                                     ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGAc                 >  4:218099/2‑146 (MQ=255)
                                                                                                                                     ctgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                        <  8:74834/139‑1 (MQ=255)
                                                                                                                                     ctgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                        >  7:74834/1‑139 (MQ=255)
                                                                                                                                      tgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTg             >  7:60211/1‑149 (MQ=255)
                                                                                                                                       gcCCTGTCTGCCAGCCAGGGCAATGGGCATCAGGGCCAGGCACAACCGCCGGTACGGGATGTTTTTCTCCCGGCGCGGGGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            <  7:270866/149‑1 (MQ=255)
                                                                                                                                        cccTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTggg           >  6:235026/1‑149 (MQ=255)
                                                                                                                                        cccTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            <  8:60211/148‑1 (MQ=255)
                                                                                                                                         ccTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGc          <  6:222394/149‑1 (MQ=255)
                                                                                                                                             gCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            <  6:67307/143‑1 (MQ=255)
                                                                                                                                             gCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            <  2:297350/143‑1 (MQ=255)
                                                                                                                                             gCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            >  5:67307/1‑143 (MQ=255)
                                                                                                                                             gCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg            >  1:297350/1‑143 (MQ=255)
                                                                                                                                               agccagccAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTc    >  6:137562/1‑149 (MQ=255)
                                                                                                                                                 ccagccagGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtggg                         <  7:159634/126‑1 (MQ=255)
                                                                                                                                                 ccagccagGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtggg                         >  8:159634/1‑126 (MQ=255)
                                                                                                                                                 ccagccagGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCa  <  7:201410/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTCCA  >  NC_002947/675444‑675737

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