Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,407,875 G→A 8.1% G207E (GGG→GAG)  PP_3020 → serine/threonine protein phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,407,8750GA 8.1% 138.3 / 2.5 49G207E (GGG→GAG) PP_3020serine/threonine protein phosphatase
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base A (2/2);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

CTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGC  >  NC_002947/3407728‑3408023
                                                                                                                                                    |                                                                                                                                                    
cTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACgggg                                                                                                                                                     <  7:14450/149‑1 (MQ=255)
      tGCGCCAAGTGTGCCAGTGGTCGCGGCGGCGCCGGAAAGAGGGGGACACCCAACCGGTGCCGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTgg                                                                                                                                               <  6:19234/149‑1 (MQ=255)
         gcCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACggggg                                                                                                                                                    <  4:69195/141‑1 (MQ=255)
         gcCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACggggg                                                                                                                                                    >  3:69195/1‑141 (MQ=255)
           cAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCTAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAg                                                                                                                                          >  6:102159/1‑149 (MQ=255)
            aaGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGc                                                                                                                                         >  5:358835/1‑149 (MQ=255)
              gtgtGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCgg                                                                                                                                       <  4:177176/149‑1 (MQ=255)
                 tgCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCa                                                                                                                                    <  4:46790/149‑1 (MQ=255)
                 tgCCAGTGGTCACGTCGGCGCCTGAAAGAGGGTGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGcc                                                                                                                                     <  2:274805/148‑1 (MQ=255)
                   ccAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCAtt                                                                                                                                  >  1:234641/1‑149 (MQ=255)
                             cgTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCTGGCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTgg                                                                                                                         >  1:98652/1‑149 (MQ=255)
                             cgTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGt                                                                                                                        <  8:214495/149‑1 (MQ=255)
                             cgTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGt                                                                                                                        <  6:358835/149‑1 (MQ=255)
                                    gcCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATg                                                                                                                 >  6:182326/1‑149 (MQ=255)
                                     ccTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATgc                                                                                                                <  7:363219/149‑1 (MQ=255)
                                      cTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCa                                                                                                                                    >  6:16926/1‑128 (MQ=255)
                                      cTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCa                                                                                                                                    <  5:16926/128‑1 (MQ=255)
                                         aaaGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCAc                                                                                                            <  8:13961/149‑1 (MQ=255)
                                         aaaGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCAc                                                                                                            >  4:33689/1‑149 (MQ=255)
                                                   acacCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGcaacaa                                                                                                  <  7:39857/149‑1 (MQ=255)
                                                       ccAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGt                                                                                                >  2:371305/1‑147 (MQ=255)
                                                       ccAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCa                                                                                              >  8:32237/1‑149 (MQ=255)
                                                         aaCCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGc                                                                                            <  7:151726/149‑1 (MQ=255)
                                                               tGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCggg                                                                                      <  5:182326/149‑1 (MQ=255)
                                                                    ggCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTg                                                                                                                          >  6:81940/1‑108 (MQ=255)
                                                                    ggCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTg                                                                                                                          <  5:81940/108‑1 (MQ=255)
                                                                     gCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCa                                                                                   <  8:187391/146‑1 (MQ=255)
                                                                     gCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCa                                                                                   >  7:187391/1‑146 (MQ=255)
                                                                         gCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTCGACACGGGGGGCTGGGCCAGCGGCCATTTCACGTTGGTCGACATGCGCACCTTTCACCAAGTCAGCCACGGGCCAGCTgcag                                                                            >  6:272387/1‑149 (MQ=255)
                                                                                    ctgGTCGGCCACACCCCGGTGCTGGAGGCTAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGc                                                                   <  5:102159/148‑1 (MQ=255)
                                                                                    ctgGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGcc                                                                  <  7:32237/149‑1 (MQ=255)
                                                                                       gTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGcccc                                                                <  2:234641/148‑1 (MQ=255)
                                                                                        tCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCggg                                                                                      <  4:149561/125‑1 (MQ=255)
                                                                                        tCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCggg                                                                                      >  3:149561/1‑125 (MQ=255)
                                                                                        tCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCg                                                              >  4:146960/1‑149 (MQ=255)
                                                                                             cacaCCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCaa                                                          >  6:244465/1‑148 (MQ=255)
                                                                                             cacaCCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAAc                                                         <  2:98652/149‑1 (MQ=255)
                                                                                                  cccGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCAc                                                    <  3:146960/149‑1 (MQ=255)
                                                                                                  cccGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCAc                                                    >  2:256928/1‑149 (MQ=255)
                                                                                                     ggTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAg                                                 >  5:34552/1‑149 (MQ=255)
                                                                                                         cTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCAt                                              <  1:371305/148‑1 (MQ=255)
                                                                                                           ggAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGAGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGccccc                                                               <  3:286894/129‑1 (MQ=255)
                                                                                                           ggAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGAGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGccccc                                                               <  3:281849/129‑1 (MQ=255)
                                                                                                           ggAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGAGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGccccc                                                               >  4:281849/1‑129 (MQ=255)
                                                                                                           ggAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGAGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGccccc                                                               >  4:286894/1‑129 (MQ=255)
                                                                                                                 cAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGa                                     <  5:244465/149‑1 (MQ=255)
                                                                                                                  aaGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAg                                    <  7:26720/149‑1 (MQ=255)
                                                                                                                       gctgGTCAATGTGTGGCACCTGGACTCGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACATTTCAGCCCTGGTCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCg                               <  5:272387/149‑1 (MQ=255)
                                                                                                                       gctgGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCCCCTTGCAACCAGTCAGCCCCGGGCCCGAGGCAGCAACACCGCCGCTGCGACAGTCCCCGGCATGCCCCCCGCGGcccc                               >  6:269885/1‑148 (MQ=255)
                                                                                                                                  tgtgGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCAc                     <  3:33689/148‑1 (MQ=255)
                                                                                                                                      gCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCa                 <  6:34552/148‑1 (MQ=255)
                                                                                                                                          cTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGCGGCCCGCTCCCGCCACCCCAGGCTg            >  4:144251/1‑149 (MQ=255)
                                                                                                                                                    ggggCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCAAAAGCAAGACACCCGAGGCCCGCGCCAGCCACCCCCGGCTCCggcctcagc  >  4:50419/1‑142 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGAAGTGCGCCAAGTGTGCCAGTGGTCGCGTCGGCGCCTGAAAGAGGGGGACACCCAACCGGTGCAGGGCCTGCGAGCCCT‑GCTGGTCGGCCACACCCCGGTGCTGGAGGCCAAGCTGCTGGGCAATGTGTGGCACCTGGACACGGGGGGCTGGGCCAGCGGCCATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGC  >  NC_002947/3407728‑3408023

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