Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,548 T→G 9.5% E51A (GAG→GCG)  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,5480TG 9.5% 144.4 / 8.2 53E51A (GAG→GCG) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (25/22);  new base G (2/3);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366410‑1366694
                                                                                                                                          |                                                                                                                                                  
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                                                                                                                                >  5:228894/1‑115 (MQ=255)
                                            gCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                <  6:228894/115‑1 (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                                                                                                   >  7:154101/1‑83 (MQ=255)
                                                                                                         tCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGcc                                                                                                   <  8:154101/83‑1 (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)
                                                                                                                                          |                                                                                                                                                  
CGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366410‑1366694

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