Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,211,889 C→A 7.4% S309* (TCA→TAA) ‡ PP_2815 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,211,8890CA 7.4% 168.3 / 3.6 54S309* (TCA→TAA) ‡PP_2815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/27);  new base A (2/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GGACCTGAACCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCGTTGCACGCCTTCC  >  NC_002947/3211743‑3212032
                                                                                                                                                  |                                                                                                                                               
ggACCTGAACCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCATTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCAc                                                                                                                                               >  7:23642/1‑149 (MQ=255)
  aCCTGAACCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTAGCGGGTTCAccc                                                                                                                                             <  3:308813/149‑1 (MQ=255)
   ccTGAACCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACccg                                                                                                                                            >  4:350602/1‑149 (MQ=255)
       aaCCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACccgcc                                                                                                                                          <  2:71189/147‑1 (MQ=255)
       aaCCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGt                                                                                                                                        >  4:225440/1‑149 (MQ=255)
        aCCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGGAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGtt                                                                                                                                       <  8:23642/149‑1 (MQ=255)
         ccTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTGAACccgccaaac                                                                                                                                      <  5:62243/149‑5 (MQ=255)
             ctgcAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCAc                                                                                                                                  <  3:37491/149‑1 (MQ=255)
             ctgcAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTGAACccgccaaacggc                                                                                                                                   >  2:44543/1‑141 (MQ=255)
             ctgcAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTGAACccgccaaacggc                                                                                                                                   <  1:44543/148‑8 (MQ=255)
                 aGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCCGGCATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGCGTTCACCCGCCGTTGCCAcgatc                                                                                                                              >  2:46728/1‑149 (MQ=255)
                  ggACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGCGTTCACCCGCCGTTGCCAcgatcg                                                                                                                             <  7:60131/149‑1 (MQ=255)
                   gACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCAcgatcga                                                                                                                            >  5:158318/1‑149 (MQ=255)
                        cGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAgg                                                                                                                       >  7:310131/1‑149 (MQ=255)
                                aaTTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATg                                                                                                               <  5:274284/149‑1 (MQ=255)
                                       gAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGc                                                                                                        >  6:29389/1‑149 (MQ=255)
                                          gATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATTAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCtga                                                                                                     <  8:271725/149‑1 (MQ=255)
                                          gATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCtga                                                                                                     <  2:201329/149‑1 (MQ=255)
                                          gATACGCGCGGCAATAGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCtga                                                                                                     <  8:22103/149‑1 (MQ=255)
                                           aTACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCtgat                                                                                                    <  6:158318/149‑1 (MQ=255)
                                            tACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCtgatg                                                                                                   <  3:350602/149‑1 (MQ=255)
                                                    gCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGa                                                                                           <  3:225440/149‑1 (MQ=255)
                                                           gTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAgg                                                                                            >  6:69055/1‑141 (MQ=255)
                                                           gTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAgg                                                                                            <  5:69055/141‑1 (MQ=255)
                                                           gTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCcgc                                                                                    >  4:233173/1‑149 (MQ=255)
                                                           gTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGAGGTTGATGGGCCGTTTGGCGGCTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCcgc                                                                                    >  8:100406/1‑149 (MQ=255)
                                                               ggcgggcgCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCcgcgcg                                                                                 <  7:3352/148‑1 (MQ=255)
                                                                gcgggcgCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGc                                                                               >  2:317745/1‑149 (MQ=255)
                                                                  gggcgCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCa                                                                             <  6:140020/149‑1 (MQ=255)
                                                                     cgcgGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGc                                                                          >  3:245272/1‑149 (MQ=255)
                                                                          ggATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTc                                                                     <  2:99566/149‑1 (MQ=255)
                                                                             tGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGGCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTc                                                                  <  3:297770/149‑1 (MQ=255)
                                                                                        cGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGtt                                                       <  4:140500/149‑1 (MQ=255)
                                                                                         ggCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTg                                                      >  2:339530/1‑149 (MQ=255)
                                                                                             cTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGgc                                                  >  4:111698/1‑149 (MQ=255)
                                                                                                      ggAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGa                                         >  5:32622/1‑149 (MQ=255)
                                                                                                      ggAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGa                                         >  5:223161/1‑149 (MQ=255)
                                                                                                          gCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGa                                     <  6:223161/149‑1 (MQ=255)
                                                                                                            cTTGACCACCCAGCGGTTGTGGGGCCGTTTGGCGGGTTCACCGGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGc                                   <  1:429167/149‑1 (MQ=255)
                                                                                                            cTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCtt                                                                >  6:48991/1‑120 (MQ=255)
                                                                                                            cTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCtt                                                                <  5:48991/120‑1 (MQ=255)
                                                                                                              tGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAgg                                                                           <  1:377817/107‑1 (MQ=255)
                                                                                                              tGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAgg                                                                           >  2:377817/1‑107 (MQ=255)
                                                                                                               gACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTgg                                            >  3:210982/1‑137 (MQ=255)
                                                                                                               gACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTgg                                            <  4:210982/137‑1 (MQ=255)
                                                                                                                    cccAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTg                           >  6:341046/1‑149 (MQ=255)
                                                                                                                        gCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGGCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTc                       <  8:218027/149‑1 (MQ=255)
                                                                                                                         cGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCa                      <  6:264657/149‑1 (MQ=255)
                                                                                                                          ggTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAg                     <  7:100406/149‑1 (MQ=255)
                                                                                                                                tGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCg               >  6:61008/1‑149 (MQ=255)
                                                                                                                                  ggCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCGtt             <  3:233173/149‑1 (MQ=255)
                                                                                                                                     gcaacggCGGGTGAACCCGTCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACgg                             >  1:206334/6‑130 (MQ=255)
                                                                                                                                     cGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCGTTGCa          >  5:418818/1‑149 (MQ=255)
                                                                                                                                             gggTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCGTTGCACGCCTTcc  >  2:227644/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGACCTGAACCTGCTGCAGGACGTCGCCCGGCAATTTCCGAAGATACGCGCGGCAATCGGTTCGGCGGGCGCGGGGATGAATGGCTTCCGGCGCTCGCACCTGGAAGCCTTGACCACCCAGCGGTTGATGGGCCGTTTGGCGGGTTCACCCGCCGTTGCCACGATCGATCAGGTGAGGATGGTGTCGCTGATGACCCAGGATGCCCGCGCGGCCAGGCAGTTCGTCTTGAGCACGTTGGGGCGCTTGGCGACCGAGCCCACGGTGCTTCAGCATTCGTTGCACGCCTTCC  >  NC_002947/3211743‑3212032

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