Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,200,558 C→T 7.8% L72L (CTG→TTG)  gspG → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,200,5580CT 7.8% 145.7 / 2.2 51L72L (CTG→TTG) gspGprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (2/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.00e-01

ATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCATGCCCG  >  NC_002947/1200414‑1200705
                                                                                                                                                |                                                                                                                                                   
aTGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTgc                                                                                                                                                  <  5:226479/148‑1 (MQ=255)
 tggtggTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCATCGTGCTGGGCAACCATGAGAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTgcgc                                                                                                                                                <  4:234646/149‑1 (MQ=255)
  ggtggtGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTgc                                                                                                                                                  <  4:14/146‑1 (MQ=255)
           cTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGCCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTCGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGc                                                                                                                                       >  5:127353/1‑148 (MQ=255)
                  atcGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGaa                                                                                                                                 >  7:217986/1‑147 (MQ=255)
                  atcGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACa                                                                                                                               <  7:117203/149‑1 (MQ=255)
                    cGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAgg                                                                                                                             >  1:78834/1‑149 (MQ=255)
                     ggCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAggg                                                                                                                            >  8:10082/1‑149 (MQ=255)
                           ctgATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTgg                                                                                                                       <  3:116388/148‑1 (MQ=255)
                            tgATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGAATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGcc                                                                                                                     >  3:183260/1‑149 (MQ=255)
                                   cGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTgg                                                                                                              >  7:95312/1‑149 (MQ=255)
                                    gTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGt                                                                                                             <  7:10082/149‑1 (MQ=255)
                                      ctggttgcCCAGCACGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTg                                                                                                            >  6:54694/7‑148 (MQ=255)
                                      ctggttgcCCAGCACGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTg                                                                                                            >  6:49810/7‑148 (MQ=255)
                                        tAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGaa                                                                                                          >  5:82543/1‑148 (MQ=255)
                                             cccAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGcc                                                                                                    <  6:82543/149‑1 (MQ=255)
                                              ccAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCg                                                                                                   >  8:195488/1‑149 (MQ=255)
                                              ccAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCg                                                                                                   <  8:263549/149‑1 (MQ=255)
                                                             aCCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGc                                                                                                                           <  8:9436/110‑1 (MQ=255)
                                                             aCCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGc                                                                                                                           >  7:9436/1‑110 (MQ=255)
                                                              ccAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGc                                                                                   >  5:52444/1‑149 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACTTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGa                                                                                           >  5:172725/1‑138 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACTTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGa                                                                                           <  4:275990/138‑1 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACTTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGa                                                                                           >  3:275990/1‑138 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACTTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGa                                                                                           <  6:172725/138‑1 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTc                                                                                <  5:49810/149‑1 (MQ=255)
                                                                 ggACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTc                                                                                <  5:54694/149‑1 (MQ=255)
                                                                  gACAAGGCCATGAAGCAGAAGGTGATGGCCGACGTGGCCACCGTGGAGCTGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCg                                                                               <  6:127353/149‑1 (MQ=255)
                                                                   aCAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCgg                                                                              <  2:78834/149‑1 (MQ=255)
                                                                        gCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTg                                                                         >  6:308905/1‑149 (MQ=255)
                                                                               aGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGt                                                                  <  8:95312/149‑1 (MQ=255)
                                                                               aGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGt                                                                  <  6:52444/149‑1 (MQ=255)
                                                                                      ggTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGc                                                            >  8:27012/1‑148 (MQ=255)
                                                                                          aTGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGgcgc                                                                     >  8:154410/1‑135 (MQ=255)
                                                                                          aTGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGgcgc                                                                     <  7:154410/135‑1 (MQ=255)
                                                                                           tGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTg                                                      <  1:257835/149‑1 (MQ=255)
                                                                                           tGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTg                                                      >  2:257835/1‑149 (MQ=255)
                                                                                              ccGACCTGGCCACCATGGAGCAGGAGCTGGACATGTATAGCCTGGACAACCTGCGGTTCCCGGGCAAGGAACAGGGCATGGCCGAGCTGGTGAAAAAGCCGGCTCAGGAGCAGCTGACTCGGGCCTGGCGCAGTGATGGCTCCGTGCGt                                                   >  6:71808/1‑149 (MQ=255)
                                                                                                      gCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGcc                                           <  7:195488/149‑1 (MQ=255)
                                                                                                                               tGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTAt                  <  7:27012/149‑1 (MQ=255)
                                                                                                                                gTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTg                               <  8:165094/135‑1 (MQ=255)
                                                                                                                                gTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTg                               >  7:165094/1‑135 (MQ=255)
                                                                                                                                   tCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCt                                              >  1:60149/1‑117 (MQ=255)
                                                                                                                                   tCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCt                                              <  2:60149/117‑1 (MQ=255)
                                                                                                                                   tCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAg               >  5:73545/1‑148 (MQ=255)
                                                                                                                                      ccTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTAcc           >  4:308421/1‑149 (MQ=255)
                                                                                                                                       cTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCg          >  3:101843/1‑149 (MQ=255)
                                                                                                                                         ggACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCa        >  6:30005/1‑149 (MQ=255)
                                                                                                                                          gACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCAccc                       <  3:85758/133‑1 (MQ=255)
                                                                                                                                          gACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCAccc                       >  4:85758/1‑133 (MQ=255)
                                                                                                                                              aCCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCATGccc   >  2:307171/1‑149 (MQ=255)
                                                                                                                                               ccTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCATGCCCg  >  4:212961/1‑149 (MQ=255)
                                                                                                                                               ccTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCATGCCCg  >  8:35049/1‑149 (MQ=255)
                                                                                                                                                cTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGAc                                    <  8:191280/114‑1 (MQ=255)
                                                                                                                                                cTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGAc                                    >  7:191280/1‑114 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATGGTGGTGATCTTTATCATCGGCCTGCTGATTGCCGTGGTAGCGCCCAGCGTGCTGGGCAACCAGGACAAGGCCATGAAGCAGAAGGTGATGGCCGACCTGGCCACCCTGGAGCAGGCGCTGGACATGTATCGCCTGGACAACCTGCGCTTCCCCAGCAACGAACAGGGCCTGGCCGCGCTGGTGAAAAAGCCGGCTCAGGAGCCGCTGCCTCGGGCCTGGCGCAGTGATGGCTACGTGCGTCGCCTGCCGGAAGACCCGTGGGGCACCCCTTATCAGTACCGCATGCCCG  >  NC_002947/1200414‑1200705

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