Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 675,587 A→T 12.7% C38S (TGC→AGC)  PP_0577 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947675,5870AT12.7% 109.9 / 2.7 43C38S (TGC→AGC) PP_0577hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base T (3/3);  total (22/21)
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

AACTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTC  >  NC_002947/675439‑675735
                                                                                                                                                    |                                                                                                                                                    
aaCTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGca                                                                                                                                                      <  3:16295/149‑1 (MQ=255)
     gCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGcagcca                                                                                                                                                  <  1:121153/148‑1 (MQ=255)
     gCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCAGGGcagcaa                                                                                                                                                  <  7:74502/148‑3 (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                                                                                         <  4:182268/126‑1 (MQ=255)
                                                                                    gccgccTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTcc                                                                                         >  3:182268/1‑126 (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               >  4:218099/2‑146 (MQ=255)
                                                                                                                                          ttgcCCTGGCTGGCTGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGAc               <  3:218099/145‑1 (MQ=255)
                                                                                                                                          ctgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                      >  7:74834/1‑139 (MQ=255)
                                                                                                                                          ctgcCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAggtgggg                      <  8:74834/139‑1 (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          <  2:297350/143‑1 (MQ=255)
                                                                                                                                                  gCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTgg          <  6:67307/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)
                                                                                                                                                    |                                                                                                                                                    
AACTGGCCGTTGCGCCAGCTCCAGCTGCGGTAGTGCCCATAATCACGAGCGGGCGCAGGGTAAGCACTGGCGTCGAAGGTGGTGGCCGCCTGGGCGGGGGCCGGTGGTAACGGGCGGCTGCTGGCCACATACGGGTTGCTGCCCTGGCAGCCAGCCAGGGCAATGGGCAGCAGGGCCAGGCACAACAGACGGTACGGCATGTATTTCTCCCGGCGCGGCGGATTAGCGGGGGCGGCAAACCCAGTGCAGGTAGCGGCCAAGCCCGGCGAAGGTGGGGTGGCGACGGTGGGCCAGTTC  >  NC_002947/675439‑675735

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