Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,367 A→C 7.7% L107R (CTG→CGG)  PP_2202 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,3670AC 7.7% 152.6 / 4.7 52L107R (CTG→CGG) PP_2202membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base C (2/2);  total (28/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 = 8.55e-01

CACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAG  >  NC_002947/2508229‑2508515
                                                                                                                                          |                                                                                                                                                     
cACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCgg                                                                                                                                             <  1:328598/149‑1 (MQ=255)
 aCGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCggg                                                                                                                                            <  5:225693/149‑1 (MQ=255)
  cgcgAAACCTGAAAGGTGTTTCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGTGTCACCATCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGt                                                                                                                                           <  5:27573/149‑1 (MQ=255)
   gcgAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTc                                                                                                                                          <  1:8884/149‑1 (MQ=255)
     gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTaa                                                                                                                                                 <  1:202711/140‑1 (MQ=255)
     gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTaa                                                                                                                                                 >  2:202711/1‑140 (MQ=255)
     gAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGc                                                                                                                                        >  2:249086/1‑149 (MQ=255)
          cTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATc                                                                                                                                   >  4:247811/1‑149 (MQ=255)
             aaaGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGTATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCAc                                                                                                                                 <  1:85226/148‑1 (MQ=255)
                 gtgtTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGc                                                                                                                             >  8:232551/1‑148 (MQ=255)
                      aCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAg                                                                                                                       >  1:104823/1‑149 (MQ=255)
                      aCCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAg                                                                                                                       >  1:286836/1‑149 (MQ=255)
                        cAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCgg                                                                                                                                     >  5:156643/1‑132 (MQ=255)
                        cAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCgg                                                                                                                                     <  6:156643/133‑2 (MQ=255)
                            cGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGagc                                                                                                                 >  4:212015/1‑149 (MQ=255)
                             gTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGCCATCCAGCTGagcc                                                                                                                >  6:351523/1‑148 (MQ=255)
                                  aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGAt                                                                                                                                    <  7:163384/124‑1 (MQ=255)
                                  aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGAt                                                                                                                                    >  8:163384/1‑124 (MQ=255)
                                  aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCg                                                                                                                                      <  2:223963/122‑1 (MQ=255)
                                  aCCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGACCCGCTTACGCTGATCGCCg                                                                                                                                      >  1:223963/1‑122 (MQ=255)
                                     aGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGAc                                                                                                       <  2:319869/149‑1 (MQ=255)
                                      gCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGAc                                                                                                       >  8:4353/1‑148 (MQ=255)
                                                     cGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGtt                                                                                       <  3:212015/149‑1 (MQ=255)
                                                      gCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTg                                                                                      >  2:79863/1‑149 (MQ=255)
                                                          gCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAg                                                                                  >  3:121517/1‑149 (MQ=255)
                                                                   aGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCtgctgc                                                                         >  7:177944/1‑149 (MQ=255)
                                                                           gCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                   <  7:81815/147‑1 (MQ=255)
                                                                           gCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCAcgc                                                                   >  8:81815/1‑147 (MQ=255)
                                                                             cTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCa                                                               <  2:286836/149‑1 (MQ=255)
                                                                                  ccAGGGTCACCCGCAACCGGAACCGCCAGAATGGCTCAACCATTAACCCGGCGATCACCGTAAGCGGGGCGCCGATCACCGGCATCCAGCTGACCATCATCAGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTg                                                           <  1:249086/149‑1 (MQ=255)
                                                                                        tCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAgcg                                                    >  8:17633/1‑149 (MQ=255)
                                                                                          aCCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCt                                                  >  7:154719/1‑149 (MQ=255)
                                                                                              gCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGaaa                                              <  1:79863/149‑1 (MQ=255)
                                                                                                   aGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGaa                                         >  7:271211/1‑149 (MQ=255)
                                                                                                         cgcCAGAATGGCTAACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgtg                                   >  4:249445/1‑147 (MQ=255)
                                                                                                          gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCgctg                                                                             <  4:198746/106‑1 (MQ=255)
                                                                                                          gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCgctg                                                                             >  3:198746/1‑106 (MQ=255)
                                                                                                          gcCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAgcgc                                  >  4:38614/1‑149 (MQ=255)
                                                                                                           ccAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAg                                                                                  >  8:202512/1‑100 (MQ=255)
                                                                                                            cAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCtt                                >  5:127870/1‑149 (MQ=255)
                                                                                                                aTGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTcg                            <  8:154719/149‑1 (MQ=255)
                                                                                                                   gCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCgcgc                         <  8:271211/149‑1 (MQ=255)
                                                                                                                       aCGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGTGCTTGTCGCGCAg                       <  3:249445/147‑1 (MQ=255)
                                                                                                                       aCGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAg                       <  8:177944/147‑1 (MQ=255)
                                                                                                                               tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGAt              <  4:121517/148‑1 (MQ=255)
                                                                                                                               tCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg             >  1:225881/1‑149 (MQ=255)
                                                                                                                               tCCCGGCGATCAGCGTAAGCGGGGCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATg             <  7:232551/149‑1 (MQ=255)
                                                                                                                                  cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc          >  6:135487/1‑149 (MQ=255)
                                                                                                                                  cGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTAGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcc          >  6:50883/1‑149 (MQ=255)
                                                                                                                                    gCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATggcccg        <  5:135487/149‑1 (MQ=255)
                                                                                                                                        tCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAg                       >  1:3835/1‑130 (MQ=255)
                                                                                                                                        tCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAg                       <  2:3835/130‑1 (MQ=255)
                                                                                                                                          aGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTCGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCcag  <  5:50883/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                     
CACGCGAAACCTGAAAGGTGTTACCAGGCGTGAAACCAGCCCAGGGTAATCATCGCAAGCACCGCATAGCGGCCAGCCTTGGCCAGGGTCACCAGCAACAGGAAGCGCCAGAATGGCTCACGCATTATCCCGGCGATCAGCGTAAGCGGGTCGCCGATCACCGGCATCCAGCTGAGCAGCAGCGACCATTGCCCCCAGCGTTGGTAGCGCTGCTGCGCACGCGCCAGTTGGGCAGCGCTGAAAGGGAACCAGCGCTTGTCGCGCAGGTGCTCGATGGCCCGGCCCAG  >  NC_002947/2508229‑2508515

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: