Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,551 A→T 9.4% V50D (GTC→GAC)  PP_1191 ← S4 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,366,5510AT 9.4% 141.0 / 3.4 53V50D (GTC→GAC) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base T (2/3);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACAT  >  NC_002947/1366410‑1366699
                                                                                                                                             |                                                                                                                                                    
cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctg                                                                                                                                               <  4:254328/149‑1 (MQ=255)
 gcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctgc                                                                                                                                              >  3:125430/1‑149 (MQ=255)
 gcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGACCATGTTCATGATGGCGCGGGCATCTTCAATGTCCATGCAGGCTGGCTGCTTCAGCAGCAGGTTCACCGGCTCCAGGATTTCGGCAAGGGCCACCGGCAGCAGCACGACCCgctgc                                                                                                                                              >  4:246892/1‑149 (MQ=255)
   aCGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctgctg                                                                                                                                            >  3:11333/1‑149 (MQ=255)
     gCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCGCGTCGAgctgctgcc                                                                                                                                          <  7:242412/149‑2 (MQ=255)
          cgcgGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCaa                                                                                                                                       <  7:163893/147‑1 (MQ=255)
           gcgGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                     <  5:116464/148‑1 (MQ=255)
             ggTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTg                                                                                                                                  <  4:125430/149‑1 (MQ=255)
              gTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGa                                                                                                                                 <  1:28851/149‑1 (MQ=255)
                   ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctg                                                                                                                                               <  2:190871/130‑1 (MQ=255)
                   ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTCGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctg                                                                                                                                               >  1:190871/1‑130 (MQ=255)
                                gCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGacca                                                                                                               >  1:178275/1‑149 (MQ=255)
                                    aTGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAc                                                                                                                                                     <  2:68027/107‑1 (MQ=255)
                                    aTGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAc                                                                                                                                                     >  1:68027/1‑107 (MQ=255)
                                            gCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCCACCTTGAACTGCGGCTGCTCGACCCCCCCGCCATCGa                                                                                                   >  4:32673/1‑149 (MQ=255)
                                            gCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                     <  6:228894/115‑1 (MQ=255)
                                            gCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                     >  5:228894/1‑115 (MQ=255)
                                            gCTGGCGCGGGCATCTTCACTGTCAATGCCGGCTGGATGGTTCAGCAGCAGGGTCACCGGCTCAAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAACTTTAACTGCGGCTGCTCGACCACAACGCCATCGa                                                                                                   >  3:59272/1‑149 (MQ=255)
                                             cTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAc                                                                                                  >  3:173465/1‑149 (MQ=255)
                                                gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGTGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                                 <  3:246892/147‑1 (MQ=255)
                                                  gcgGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTa                                                                                             >  6:130259/1‑149 (MQ=255)
                                                          cTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCCACCGTTACCCAGcca                                                                                     >  1:238503/1‑148 (MQ=255)
                                                                   ccATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACGGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAAt                                                                             <  3:228429/148‑1 (MQ=255)
                                                                   ccATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCATCCTTCAATg                                                                            >  8:163389/1‑149 (MQ=255)
                                                                    cATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATg                                                                            <  5:130259/148‑1 (MQ=255)
                                                                        ccGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACt                                                                        <  4:11333/148‑1 (MQ=255)
                                                                               gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                                 <  4:173465/148‑1 (MQ=255)
                                                                                       agcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggcg                                                        >  4:221712/1‑149 (MQ=255)
                                                                                       agcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggcg                                                        >  6:156348/1‑149 (MQ=255)
                                                                                         cagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                       >  1:55199/1‑148 (MQ=255)
                                                                                          agcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCTCCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgc                                                            <  8:136591/142‑1 (MQ=255)
                                                                                          agcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgc                                                            >  7:136591/1‑142 (MQ=255)
                                                                                          agcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                     >  1:326758/1‑149 (MQ=255)
                                                                                              gggTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGagcag                                                 >  8:283292/1‑149 (MQ=255)
                                                                                                 tCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGcc                                               >  8:31234/1‑148 (MQ=255)
                                                                                                 tCACCGGCTCAAGTTTTTCGGCACGGGCGCCCGGCAGCAGCTCGACTCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCTTCCTCGCGGCGGGAGCAGCCg                                              <  4:59272/149‑1 (MQ=255)
                                                                                                   aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                            <  2:178275/149‑1 (MQ=255)
                                                                                                         tCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGcc                                                                                                        <  8:154101/83‑1 (MQ=255)
                                                                                                         tCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGcc                                                                                                        >  7:154101/1‑83 (MQ=255)
                                                                                                           cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGa                                    <  5:156348/149‑1 (MQ=255)
                                                                                                           cAGGGGTGCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcc                                     <  3:32673/148‑2 (MQ=255)
                                                                                                             gggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATg                                  >  7:77052/1‑149 (MQ=255)
                                                                                                                      cACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                         >  5:126590/1‑149 (MQ=255)
                                                                                                                               ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGATGCGttt                          <  7:134834/139‑1 (MQ=255)
                                                                                                                               ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGATGCGttt                          >  8:134834/1‑139 (MQ=255)
                                                                                                                                gagcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgct                                          >  4:33336/3‑122 (MQ=255)
                                                                                                                                gagcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgct                                          <  3:33336/120‑1 (MQ=255)
                                                                                                                                  gcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGa                                    <  5:34115/126‑1 (MQ=255)
                                                                                                                                  gcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCACCCGAGCTGCTCGa                                    >  6:34115/1‑126 (MQ=255)
                                                                                                                                      cagcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGa                       <  5:238087/135‑1 (MQ=255)
                                                                                                                                      cagcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCCAACAGCTCCCTGAGGCGTTTGGa                       >  6:238087/1‑135 (MQ=255)
                                                                                                                                       agcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTc        <  7:283292/149‑1 (MQ=255)
                                                                                                                                         cTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg       >  6:261990/1‑148 (MQ=255)
                                                                                                                                             aCGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACAt  <  5:261990/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
CGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACAT  >  NC_002947/1366410‑1366699

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