Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,114,823 C→T 7.5% R115W (CGG→TGG)  PP_1883 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,114,8230CT 7.5% 150.2 / 2.0 53R115W (CGG→TGG) PP_1883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base T (2/2);  total (32/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 = 8.95e-01

CTGGCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGACAAGGTTTTC  >  NC_002947/2114683‑2114969
                                                                                                                                              |                                                                                                                                                  
cTGGCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCg                                                                                                                                              <  1:209091/148‑1 (MQ=255)
  ggCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACTGGCCCGCTg                                                                                                                                           >  8:290693/1‑148 (MQ=255)
  ggCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACTGGCCCGCTg                                                                                                                                           >  6:28081/1‑148 (MQ=255)
  ggCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTg                                                                                                                                           >  2:158384/1‑148 (MQ=255)
   gCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTg                                                                                                                                           <  5:118120/147‑1 (MQ=255)
                 cTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCATTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCtt                                                                                                                            <  7:78194/148‑1 (MQ=255)
                 cTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCATTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGATGACCGGCCCGCTGAAGGCATATGGGCtt                                                                                                                            >  8:78194/1‑148 (MQ=255)
                   gATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGc                                                                                                                         >  2:161403/1‑149 (MQ=255)
                       gATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGt                                                                                                                     >  7:100561/1‑149 (MQ=255)
                        aTAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGtt                                                                                                                    <  6:313866/149‑1 (MQ=255)
                             gCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCTGTTGccaa                                                                                                               >  8:95637/1‑147 (MQ=255)
                             gCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCa                                                                                                               <  3:236967/149‑1 (MQ=255)
                               aCTAGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACt                                                                                                             >  1:298595/1‑149 (MQ=255)
                                    tCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCaa                                                                                                        <  3:187242/149‑1 (MQ=255)
                                          ccTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATg                                                                                                    >  2:282946/1‑147 (MQ=255)
                                          ccTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATg                                                                                                    <  1:282946/147‑1 (MQ=255)
                                          ccTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGc                                                                                                   <  2:298595/148‑1 (MQ=255)
                                           cTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAg                                                                                                                                        <  5:102026/110‑1 (MQ=255)
                                           cTTGCGCGGTGGCCACGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAg                                                                                                                                        >  6:102026/1‑110 (MQ=255)
                                            ttGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACtt                                                                                                            <  3:303063/137‑1 (MQ=255)
                                            ttGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACtt                                                                                                            >  4:303063/1‑137 (MQ=255)
                                                       ccAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGa                                                                                     >  6:193298/1‑149 (MQ=255)
                                                            tgtTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGAc                                                                                <  1:161403/149‑1 (MQ=255)
                                                            tgtTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGAc                                                                                <  5:193298/149‑1 (MQ=255)
                                                                  aTTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAAc                                                                          >  4:301727/1‑149 (MQ=255)
                                                                      gTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGAGCGTGAGTCCAAGCGTCGTTTCTTGTACGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGTTTGCCCAGTTAAAAGATGCCCATGAGGGGGTGATCGAGGTCACCCTg                                                                       >  1:53344/1‑149 (MQ=255)
                                                                         cACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTc                                                                    <  1:158384/148‑1 (MQ=255)
                                                                          aCGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCaa                                                                  >  4:144357/1‑149 (MQ=255)
                                                                                 tGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTATCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCACAATCCCGaa                                                           >  6:174750/1‑149 (MQ=255)
                                                                                   gACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGaa                                                           <  8:100561/147‑1 (MQ=255)
                                                                                         gCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATAtt                                                   >  5:294765/1‑149 (MQ=255)
                                                                                         gCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATAtt                                                   >  7:15250/1‑149 (MQ=255)
                                                                                            ggCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGaa                                                           >  1:274608/1‑138 (MQ=255)
                                                                                            ggCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGaa                                                           <  2:274608/138‑1 (MQ=255)
                                                                                             gCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACTGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTg                                               <  5:28081/149‑1 (MQ=255)
                                                                                             gCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACTGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTg                                               <  7:290693/149‑1 (MQ=255)
                                                                                                      tgGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAgcg                                      >  1:166782/1‑149 (MQ=255)
                                                                                                             cTGAGTCCAAGCTTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGTCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTAt                                <  7:95637/149‑1 (MQ=255)
                                                                                                             cTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTAt                                >  2:321328/1‑149 (MQ=255)
                                                                                                             cTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTAt                                >  8:194349/1‑149 (MQ=255)
                                                                                                               gAGTCCAAGCGTCGTTTCTTCTTCGTTTACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTAAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATgc                              >  8:94706/1‑149 (MQ=255)
                                                                                                                aGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATgcc                             >  5:67905/1‑149 (MQ=255)
                                                                                                                              tttCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCa                   >  4:154654/1‑146 (MQ=255)
                                                                                                                              tttCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCa                   >  3:46489/1‑146 (MQ=255)
                                                                                                                              tttCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCa                   >  3:306723/1‑146 (MQ=255)
                                                                                                                              tttCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAg                <  6:294765/149‑1 (MQ=255)
                                                                                                                                tcttcttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGAt              >  6:7514/1‑149 (MQ=255)
                                                                                                                                 cttcttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGAt              >  8:8730/1‑148 (MQ=255)
                                                                                                                                 cttcttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGAt              <  4:306723/148‑1 (MQ=255)
                                                                                                                                 cttcttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGAGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGAt              >  2:166967/1‑148 (MQ=255)
                                                                                                                                  ttcttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGa            >  4:214972/1‑149 (MQ=255)
                                                                                                                                     ttcGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGAc           >  3:24940/1‑147 (MQ=255)
                                                                                                                                       cGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGACAAg        <  7:94706/148‑1 (MQ=255)
                                                                                                                                              cGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGACAAGGTTTTc  <  6:293433/147‑1 (MQ=255)
                                                                                                                                              cGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGACAAGGTTTTc  >  5:293433/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTGGCCGCGGGGAGATCCTGATCGATAACGCACTCGTCACTGCCTTGCGCGGTGGCCAGGTGTTCTATTCGTTCACGCGTCTGGACAACGCCGGCAATCCTGTGGGGCCTGAGTCCAAGCGTCGTTTCTTCTTCGTTGACCGGCCCGCTGAAGGCATATGGGCTTTGCCGGTTGCCCACTTTCAAGATGCCCATGAGGGGGTGATCGACGTCAACCTGCTCAATCCCGAAGGCATATTGCTGGCAACAGCGCCCTATGCCGCCATGGCCAAAGATGACAAGGTTTTC  >  NC_002947/2114683‑2114969

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: