Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,113,719 G→T 7.8% Q204K (CAG→AAG)  PP_1882 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,113,7190GT 7.8% 152.1 / 3.1 51Q204K (CAG→AAG) PP_1882membrane protein
Reads supporting (aligned to +/- strand):  ref base G (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 = 5.20e-01

CCTTGCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGCC  >  NC_002947/2113572‑2113864
                                                                                                                                                   |                                                                                                                                                 
ccTTGCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGAt                                                                                                                                                 >  2:433214/1‑150 (MQ=255)
  ttGCACACCCACAACATCCCTTTCACCAGGGCACATTGCAGGCGCCTTCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTgcatc                                                                                                                                               <  1:458029/150‑5 (MQ=255)
    gCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAgg                                                                                                                                             <  1:723795/150‑1 (MQ=255)
      tCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGgc                                                                                                                                           >  1:154769/1‑150 (MQ=255)
            acaacaACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTcgc                                                                                                                                     <  1:196716/150‑1 (MQ=255)
                caacCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCaa                                                                                                                                 <  2:674285/150‑1 (MQ=255)
                 aacCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                <  2:790264/150‑1 (MQ=255)
                 aacCCGTACACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                <  2:154769/150‑1 (MQ=255)
                          aCCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCgg                                                                                                                       >  2:546095/1‑150 (MQ=255)
                          aCCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCgg                                                                                                                       >  2:386018/1‑150 (MQ=255)
                           ccAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTTATCAGGgcgacaacg                                                                                                                                    >  2:419635/1‑130 (MQ=255)
                           ccAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTTATCAGGgcgacaacg                                                                                                                                    <  1:419635/136‑7 (MQ=255)
                                gCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTAtcg                                                                                                                 >  2:656155/1‑150 (MQ=255)
                                gCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTAtcg                                                                                                                 >  2:563058/1‑150 (MQ=255)
                                        cAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTc                                                                                                         >  2:710376/1‑150 (MQ=255)
                                         aGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                <  1:838026/126‑1 (MQ=255)
                                         aGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAAt                                                                                                                                >  2:838026/1‑126 (MQ=255)
                                            cgcCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAtt                                                                                                     <  1:623788/150‑1 (MQ=255)
                                             gcCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAttt                                                                                                    >  1:373935/1‑150 (MQ=255)
                                                       cAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAt                                                                                                      >  1:143190/1‑138 (MQ=255)
                                                       cAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGAt                                                                                                      <  2:143190/138‑1 (MQ=255)
                                                          cTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                       <  1:563058/150‑1 (MQ=255)
                                                           tCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGTTATCGTCGCGTTCGATTTCCAGCAGGATctcc                                                                                      >  1:15128/1‑150 (MQ=255)
                                                                 gTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGc                                                                                <  2:373935/150‑1 (MQ=255)
                                                                             atcgatGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGa                                                                                                                            <  2:448662/94‑1 (MQ=255)
                                                                             atcgatGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGa                                                                                                                            >  1:448662/1‑94 (MQ=255)
                                                                                  tgtgACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCt                                                               <  1:386018/150‑1 (MQ=255)
                                                                                      aCCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGc                                                           >  1:234469/1‑150 (MQ=255)
                                                                                        ccAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGcat                                                         >  1:335480/1‑150 (MQ=255)
                                                                                               gtgtTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGGAGGATCTTCTCACGCTCGTCGCTGATCGAGCTGCGCATCATGTaa                                                  <  1:412501/150‑1 (MQ=255)
                                                                                                     aTCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAg                                            >  2:9242/1‑150 (MQ=255)
                                                                                                      tCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCTCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                     >  1:619940/1‑140 (MQ=255)
                                                                                                      tCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                     <  2:619940/140‑1 (MQ=255)
                                                                                                               tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCTCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGaa                                       <  2:174089/145‑1 (MQ=255)
                                                                                                               tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGaa                                       >  1:174089/1‑145 (MQ=255)
                                                                                                               tCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGCGGCTGCTCACGCTCGTCGATGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGtt                                  >  1:438393/1‑150 (MQ=255)
                                                                                                                   aTAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                       >  1:702693/1‑93 (MQ=255)
                                                                                                                   aTAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGAtctc                                                                                       <  2:702693/93‑1 (MQ=255)
                                                                                                                    tAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                     >  1:475465/1‑126 (MQ=255)
                                                                                                                    tAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATg                                                     <  2:475465/126‑1 (MQ=255)
                                                                                                                        gggCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGa                         >  2:799123/1‑150 (MQ=255)
                                                                                                                        gggCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCCATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGa                         <  1:9242/150‑1 (MQ=255)
                                                                                                                            aaCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCTTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGa                                                             >  1:125401/1‑110 (MQ=255)
                                                                                                                            aaCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGa                                                             <  2:125401/110‑1 (MQ=255)
                                                                                                                             aCGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCAGCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCg                    >  1:830588/1‑150 (MQ=255)
                                                                                                                                        gTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGa                                                                      >  2:290185/1‑89 (MQ=255)
                                                                                                                                        gTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGa                                                                      <  1:290185/89‑1 (MQ=255)
                                                                                                                                           agcGCCCTTATCAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGa                                             >  1:95207/3‑111 (MQ=255)
                                                                                                                                           agcGCCCTTATCAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGa                                             <  2:95207/109‑1 (MQ=255)
                                                                                                                                             cGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTgg    <  1:546095/150‑1 (MQ=255)
                                                                                                                                               cccTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGcc  <  2:335480/150‑1 (MQ=255)
                                                                                                                                               cccTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGcc  >  1:306211/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTTGCTCACGCACAACAACCCGTTCACCAGGGCGTATTCCAGGCGCCAGCTGGCCAGCTCATGGGTTGGCAGCACGATCGATGTGACCCAGTCCGTGTTGATCAGGTTCATCGCATAAAGGGCAACGGGCTGGCCGTTGTCGCCCTGATAAGGGCGCTCGCGCAATTCGATGCGGCTATCGTCGCGTTCGATTTCCAGCAGGATCTCCTCACGCTCGTCGCTGATCGAGCTGAGCATCATGTAAGACGAGTTGAACGGTTGGTACTGGAGGTCGATCCTGTGCCCCGTGGCC  >  NC_002947/2113572‑2113864

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