Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,695,944 A→T 7.7% L75Q (CTG→CAG)  ligD ← DNA ligase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,695,9440AT 7.7% 149.0 / 3.1 52L75Q (CTG→CAG) ligDDNA ligase D
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base T (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

CATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGC  >  NC_002947/3695813‑3696088
                                                                                                                                      |                                                                                                                                                
cATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTggg                                                                                                                                  <  8:169514/148‑1 (MQ=255)
   gcgcGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCAACGCCAGGCGGCGCACCTTGGGGTc                                                                                                                               >  3:155613/1‑148 (MQ=255)
           ccTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGgcggcg                                                                                                                                           >  1:150582/1‑128 (MQ=255)
           ccTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGgcggcg                                                                                                                                           <  2:150582/128‑1 (MQ=255)
                  gCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTg                                                                                                               <  6:165424/149‑1 (MQ=255)
                    tCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTg                                                                                                                                    <  4:115350/126‑1 (MQ=255)
                    tCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTg                                                                                                                                    >  3:115350/1‑126 (MQ=255)
                    tCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCACCGACGGGCCTTTggg                                                                                                             >  5:103538/1‑149 (MQ=255)
                      ccAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGa                                                                                                           >  2:25887/1‑149 (MQ=255)
                         gATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATCACAGGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGa                                                                                                           >  1:94957/1‑149 (MQ=255)
                         gATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCg                                                                                                        <  1:25887/149‑1 (MQ=255)
                               gcgGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGc                                                                                                  <  3:273266/149‑1 (MQ=255)
                                     ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCGCCGCCTGGCGGTGCatgt                                                                                                                                      <  7:2963/107‑4 (MQ=255)
                                     ccAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCGCCGCCTGGCGGTGCatgt                                                                                                                                      >  8:2963/1‑104 (MQ=255)
                                      cAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAg                                                                                           >  8:117980/1‑149 (MQ=255)
                                       aGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAgg                                                                                          <  3:170037/149‑1 (MQ=255)
                                         aCGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGt                                                                                        <  4:200000/149‑1 (MQ=255)
                                             tCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTggg                                                                                                                                  <  5:62590/103‑1 (MQ=255)
                                             tCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTggg                                                                                                                                  >  6:62590/1‑103 (MQ=255)
                                              cACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTg                                                                                       <  7:217904/145‑1 (MQ=255)
                                              cACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTg                                                                                       >  8:217904/1‑145 (MQ=255)
                                               aCGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCATCTTTTCAGGGTGCCAt                                                                                   <  6:242516/148‑1 (MQ=255)
                                                    cccTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCt                                                                             <  5:276689/149‑1 (MQ=255)
                                                        gCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTc                                                                            >  6:118107/1‑146 (MQ=255)
                                                        gCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTc                                                                            <  5:118107/146‑1 (MQ=255)
                                                          aCCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGc                                                                       <  7:117980/149‑1 (MQ=255)
                                                            cATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGAACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCg                                                                      >  2:19551/1‑148 (MQ=255)
                                                              tAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGAGCCTTTGGGGATCGCCCAGCttttg                                                                                             >  4:29020/1‑122 (MQ=255)
                                                                    ccTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTa                                                             <  1:19551/149‑1 (MQ=255)
                                                                         cGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGc                                                         <  6:190766/148‑1 (MQ=255)
                                                                               atGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGAc                                                  >  2:126950/1‑149 (MQ=255)
                                                                               atGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGAc                                                  >  5:130943/1‑149 (MQ=255)
                                                                                 ggCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTg                                                >  1:264048/1‑149 (MQ=255)
                                                                                     cTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcgt                                            >  5:142887/1‑149 (MQ=255)
                                                                                     cTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCACGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcgt                                            >  1:221687/1‑149 (MQ=255)
                                                                                     cTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAAGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGcgt                                            >  1:225757/1‑149 (MQ=255)
                                                                                                 ggCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCt                                   <  2:94957/149‑1 (MQ=255)
                                                                                                     tAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGAt                               <  6:103538/149‑1 (MQ=255)
                                                                                                       gTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGc                             <  4:193656/149‑1 (MQ=255)
                                                                                                          cAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAg                           >  1:351/1‑148 (MQ=255)
                                                                                                          cAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGATTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAg                           >  1:64221/1‑148 (MQ=255)
                                                                                                           aGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGaa                         <  4:155613/149‑1 (MQ=255)
                                                                                                           aGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGa                          >  3:159696/1‑148 (MQ=255)
                                                                                                              gggTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTg                      <  6:142887/149‑1 (MQ=255)
                                                                                                                gTGATCCTCCACATGCCCCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTg                      >  4:59959/1‑147 (MQ=255)
                                                                                                                gTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGc                                              <  3:21708/123‑1 (MQ=255)
                                                                                                                gTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGc                                              >  4:21708/1‑123 (MQ=255)
                                                                                                                gTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTg                      >  2:224585/1‑147 (MQ=255)
                                                                                                                        ccACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATgg            >  7:205597/1‑149 (MQ=255)
                                                                                                                          aggtGCGCCGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGc          >  5:142301/4‑149 (MQ=255)
                                                                                                                               cACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGtt     <  3:59959/149‑1 (MQ=255)
                                                                                                                                  cGCCTGGCGGTGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGc  <  6:142301/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                
CATGCGCGTCGCCTTCAGGCTCCCAGATACCGCGGTCCCAGACGATCACGTCCCCTGCACCATAGTGGCCTTCCGGGATATGGCCCTCGAAATC‑‑‑GGCGTAGTCCAGCGGGTGATCCTCCACATGCACCGCCAGGCGGCGCACCTTGGGGTCCAACGACGGGCCTTTGGGGATCGCCCAGCTTTTCAGGGTGCCATCGAGCTCCAGGCGGAAGTCGTAGTGCAAATGACTGGCGTCGTGTTTCTGGATGCAGAACTGCAGCGCATGGGCCCGTTTGC  >  NC_002947/3695813‑3696088

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