Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,173,012 T→G 10.8% E14A (GAG→GCG)  xseA ← exodeoxyribonuclease 7 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,173,0120TG10.8% 140.1 / 11.0 56E14A (GAG→GCG) xseAexodeoxyribonuclease 7 large subunit
Reads supporting (aligned to +/- strand):  ref base T (24/23);  new base G (3/3);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.34e-01

GGGCGCCGCTGTCCTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGCCCTGCATC  >  NC_002947/1172866‑1173159
                                                                                                                                                  |                                                                                                                                                   
gggCGCCGCTGTCCTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTcg                                                                                                                                                   >  3:243114/1‑149 (MQ=255)
           tCCTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGccc                                                                                                                                        <  2:272281/149‑1 (MQ=255)
           tCCTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGc                                                                                                                                          <  3:193487/147‑1 (MQ=255)
             cTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGa                                                                                                                                      <  2:372433/149‑1 (MQ=255)
               tGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGt                                                                                                                                    >  2:209989/1‑149 (MQ=255)
                   ggTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGtt                                                                                                                                <  2:340432/149‑1 (MQ=255)
                         gTACATGTGGCCGGACGCCGGGCGGGGGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAg                                                                                                                          <  6:299934/149‑1 (MQ=255)
                             aTGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTc                                                                                                                      >  1:264632/1‑149 (MQ=255)
                             aTGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTc                                                                                                                      >  3:367756/1‑149 (MQ=255)
                             aTGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTc                                                                                                                      >  3:136994/1‑149 (MQ=255)
                              tgtgGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCt                                                                                                                     <  5:193241/149‑1 (MQ=255)
                                  gCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTcc                                                                                                                                             <  2:141300/121‑1 (MQ=255)
                                  gCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTcc                                                                                                                                             >  1:141300/1‑121 (MQ=255)
                                    cGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGGTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGAt                                                                                                                >  5:95452/1‑148 (MQ=255)
                                             ggcgggcgAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCACGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCAtggctgg                                                                                                       >  4:220230/1‑148 (MQ=255)
                                                  gcgAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGTCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGcatcat                                                                                                 <  7:58211/149‑1 (MQ=255)
                                                   cgAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGcatcat                                                                                                 >  5:294910/1‑148 (MQ=255)
                                                    gAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATaa                                                                                               <  4:243114/149‑1 (MQ=255)
                                                        ttGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCAGGCGCGGCCGTTGAGTTGGCCGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGca                                                                                                     >  1:205459/1‑139 (MQ=255)
                                                        ttGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCAGGCGCGGCCGTTGAGTTGGCCGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGca                                                                                                     <  2:205459/139‑1 (MQ=255)
                                                          ggAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATaa                                                                                               >  8:190626/1‑143 (MQ=255)
                                                          ggAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATaa                                                                                               <  7:190626/143‑1 (MQ=255)
                                                            agaTTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATc                                                                                       >  6:335988/1‑149 (MQ=255)
                                                              aTTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCttt                                                                                                                   >  3:127693/1‑119 (MQ=255)
                                                              aTTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCttt                                                                                                                   <  4:127693/119‑1 (MQ=255)
                                                               tttCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGa                                                                                                                 <  5:368502/120‑1 (MQ=255)
                                                               tttCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTATTTGa                                                                                                                 >  6:368502/1‑120 (MQ=255)
                                                                    ccTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATc                                                                                                               >  3:361761/1‑117 (MQ=255)
                                                                    ccTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATc                                                                                                               <  4:361761/117‑1 (MQ=255)
                                                                            ccacACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTtgtctgt                                                                       <  2:264632/149‑1 (MQ=255)
                                                                              acacGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTgtctgtct                                                                     >  4:312224/1‑149 (MQ=255)
                                                                                      cGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGCCGGTAAGGACCACGCGGTCCCGCCCGAGTCTTTCGAAAGGGTATTTGATCATGGGTGGCATCATAAATGACATCGCCCCCGGTTGTATGTATTGCCCCTa                                                             >  6:364334/1‑149 (MQ=255)
                                                                                                      gcagcaCGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGAcc                                              >  6:174796/1‑148 (MQ=255)
                                                                                                             gcgGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATATATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCtt                                      >  4:13051/1‑149 (MQ=255)
                                                                                                                ggCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATATATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTtgtg                                   <  3:13051/149‑1 (MQ=255)
                                                                                                                 gcgcGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTgg                                  <  3:312224/149‑1 (MQ=255)
                                                                                                                   gcgGCCGTTGAGTTGGCTGACTGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGga                                <  5:174796/149‑1 (MQ=255)
                                                                                                                                 ggCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGt                   <  1:209989/148‑1 (MQ=255)
                                                                                                                                 ggCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGt                   <  6:294910/148‑1 (MQ=255)
                                                                                                                                 ggCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGt                   <  5:364334/148‑1 (MQ=255)
                                                                                                                                      aCGGTAAGGACCCCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCTGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCtt              <  3:220230/148‑1 (MQ=255)
                                                                                                                                          tAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAAc                       >  6:312025/1‑135 (MQ=255)
                                                                                                                                          tAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAAc                       <  5:312025/135‑1 (MQ=255)
                                                                                                                                          tAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGcc         >  4:104502/1‑149 (MQ=255)
                                                                                                                                          tAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGcc         >  8:214042/1‑149 (MQ=255)
                                                                                                                                          gctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCa                                                       >  2:208718/4‑103 (MQ=255)
                                                                                                                                          gctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCa                                                       <  1:208718/100‑1 (MQ=255)
                                                                                                                                          gctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAg           <  1:94835/144‑1 (MQ=255)
                                                                                                                                          gctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAg           >  2:94835/4‑147 (MQ=255)
                                                                                                                                           ctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTtg                                     <  7:9413/118‑1 (MQ=255)
                                                                                                                                           ctggACCGCGAGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTtg                                     >  8:9413/3‑120 (MQ=255)
                                                                                                                                             ggACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTc             <  5:332673/142‑1 (MQ=255)
                                                                                                                                             ggACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTc             >  6:332673/1‑142 (MQ=255)
                                                                                                                                              gACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGCCCTg      >  8:37141/1‑148 (MQ=255)
                                                                                                                                               aCCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGCCCTGc     >  6:204454/1‑148 (MQ=255)
                                                                                                                                                 cTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGCCCTGCAtc  >  3:64410/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGGCGCCGCTGTCCTTGAGGGTGAAGTACATGTGGCCGGACGCCGGGCGGGCGAGGTTGGAGATTTCGCCTTCCACCCACACGCTGCGGAAAACGTCTTCCAGCAGCACGCGGGCGCGGCCGTTGAGTTGGCTGACGGTAAGGACCTCGCGGTCCAGCCCGAGTCGTTCGAAAGGGTCTTTGATCATGGCTGGCATCATAAATGACATCGGCCCCGGTTGTCTGTCTTGCCCCTAGTGTCAGCCTTGACCCAATCCTTGTGGGAGACCGCAACCGGTGCCTTCAGCCCTGCATC  >  NC_002947/1172866‑1173159

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