Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,690,488 A→G 6.8% K257E (AAG→GAG)  PP_3255 → non‑homologous end joining protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,690,4880AG 6.8% 172.6 / 2.6 59K257E (AAG→GAG) PP_3255non‑homologous end joining protein
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base G (2/2);  total (30/29)
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.90e-01

AGGATTATCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTCGACTGCT  >  NC_002947/3690345‑3690634
                                                                                                                                               |                                                                                                                                                  
aGGATTATCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAgc                                                                                                                                                 <  3:56352/147‑1 (MQ=255)
     tATCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCt                                                                                                                                          >  1:166819/1‑149 (MQ=255)
       tCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTgg                                                                                                                                        >  4:6889/1‑149 (MQ=255)
       tCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTg                                                                                                                                         <  2:225199/148‑1 (MQ=255)
         aTGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGt                                                                                                                                      >  4:123907/1‑149 (MQ=255)
          tGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTg                                                                                                                                     >  5:225789/1‑149 (MQ=255)
             cGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTggg                                                                                                                                       >  1:289156/1‑144 (MQ=255)
             cGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTggg                                                                                                                                       <  2:289156/144‑1 (MQ=255)
              gCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCaa                                                                                                                                 <  1:137767/149‑1 (MQ=255)
              gCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCaa                                                                                                                                 >  2:255071/1‑149 (MQ=255)
              gCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATAGAGGCGGTGGAGAAGGGAGAAACAGGGGCGGCCGAGAAAGGTCCGGACATCCTCGACCTGACCTAACTGCTCAAGCGCAGCCTGGGTGGaaa                                                                                                                                 >  8:312305/1‑149 (MQ=255)
                 ttCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGa                                                                                                                               >  3:199064/1‑148 (MQ=255)
                     gCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCg                                                                                                                          <  6:225789/149‑1 (MQ=255)
                          aCGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTg                                                                                                                       >  7:205380/1‑147 (MQ=255)
                            gATCATGGACCTGGTGGAAGAGAAGGCCTGCAAGGGCAAGATCGAGCCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTgaaga                                                                                                                   <  6:237271/149‑1 (MQ=255)
                              tCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTgaag                                                                                                                    >  7:171102/1‑146 (MQ=255)
                                   gACCTTGTGGAAGAGAAGTCCAGCAAGGGCAAGATCGAGGCGTTGGAGAATGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGAcgc                                                                                                             <  7:312305/148‑1 (MQ=255)
                                             aaGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGTCCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCa                                                                                                  <  1:255071/149‑1 (MQ=255)
                                              agagAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCAt                                                                                                 <  3:6889/149‑1 (MQ=255)
                                                agaAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGTGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATg                                                                                                >  5:187631/1‑148 (MQ=255)
                                                    ggCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGaaa                                                                                                         <  6:232823/135‑1 (MQ=255)
                                                    ggCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGaaa                                                                                                         >  5:232823/1‑135 (MQ=255)
                                                            aGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATg                                                                                   <  2:166819/149‑1 (MQ=255)
                                                              ggCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGcc                                                                                 <  8:205380/149‑1 (MQ=255)
                                                               gCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCa                                                                                <  3:123907/149‑1 (MQ=255)
                                                                 aaGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCAtt                                                                              <  2:120068/149‑1 (MQ=255)
                                                                   gATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGt                                                                                                                                      <  7:86462/91‑1 (MQ=255)
                                                                   gATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGt                                                                                                                                      >  8:86462/1‑91 (MQ=255)
                                                                   gATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGc                                                                                                                                   >  1:1245/1‑94 (MQ=255)
                                                                   gATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGc                                                                                                                                   <  2:1245/94‑1 (MQ=255)
                                                                        aGGCGGTGGAGAAGGGTGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCa                                                                       <  6:187631/149‑1 (MQ=255)
                                                                         ggCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCaa                                                                      <  8:171102/149‑1 (MQ=255)
                                                                               ggAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgc                                                                >  3:50582/1‑149 (MQ=255)
                                                                                 agaAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgcgc                                                              <  4:183861/149‑1 (MQ=255)
                                                                                              aGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACccagc                                                  >  7:289159/1‑148 (MQ=255)
                                                                                                        gagaAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCAc                                       >  4:346026/1‑149 (MQ=255)
                                                                                                        gagaAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCAc                                       <  3:346026/149‑1 (MQ=255)
                                                                                                           aaaGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGACGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCaa                                    <  4:199064/149‑1 (MQ=255)
                                                                                                                     ccctCCTCGACCTGACCGAACTGCTCGAGCGCAGACTGGGTGGCAAGCAACCGGTGAAGAAGACGCGAAAGGCGTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgcgc                                                              >  5:301878/4‑113 (MQ=255)
                                                                                                                     aCATCCTCGCCCTGACCGAACTGCTCGAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgcgc                                                              <  6:301878/113‑1 (MQ=255)
                                                                                                                     aCATCCTCGACCTGACCGAACTGCTCGAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgcgc                                                              >  5:301875/1‑113 (MQ=255)
                                                                                                                     aCATCCTCGACCTGACCGAACTGCTCGAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTgcgc                                                              <  6:301875/113‑1 (MQ=255)
                                                                                                                          cTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCgg                         >  1:141711/1‑145 (MQ=255)
                                                                                                                          cTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCgg                         <  2:141711/145‑1 (MQ=255)
                                                                                                                          cTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGgcgc                     >  2:224839/1‑149 (MQ=255)
                                                                                                                           tCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCt                                                            >  5:194983/1‑109 (MQ=255)
                                                                                                                           tCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCt                                                            <  6:194983/109‑1 (MQ=255)
                                                                                                                              aCCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCt                 >  6:136632/1‑149 (MQ=255)
                                                                                                                                   aCCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCAc            >  7:135629/1‑149 (MQ=255)
                                                                                                                                      gAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTc         >  3:113519/1‑149 (MQ=255)
                                                                                                                                       aaCTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGa                                                                            <  3:220952/81‑1 (MQ=255)
                                                                                                                                       aaCTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGa                                                                            >  4:220952/1‑81 (MQ=255)
                                                                                                                                          tGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGgcgc                     >  4:46941/1‑133 (MQ=255)
                                                                                                                                          tGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGgcgc                     <  3:46941/133‑1 (MQ=255)
                                                                                                                                            cTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTCGActg    <  4:50582/148‑1 (MQ=255)
                                                                                                                                             tCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCgg                   <  3:86348/132‑1 (MQ=255)
                                                                                                                                             tCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCgg                   >  4:86348/1‑132 (MQ=255)
                                                                                                                                             tCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTCGActgct  >  2:355924/1‑149 (MQ=255)
                                                                                                                                             tCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTCGActgct  >  8:253570/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
AGGATTATCATGACGCCTTCCGCCAGACGATCATGGACCTGGTGGAAGAGAAGGCCAGCAAGGGCAAGATCGAGGCGGTGGAGAAGGGCGACACAGGGGCGGCCGAGAAAGGTGCGGACATCCTCGACCTGACCGAACTGCTCAAGCGCAGCCTGGGTGGCAAGAAACCGGTGAAGAAGACGCGAAAGGCTTCATGAACACGGCCAATGCCATTGAGTTCAAGACTGCGCCTGCAACCCAGCCAGGTAATCACCAAACCCCTCGCGGGCGCGGCTGTCACGTCGACTGCT  >  NC_002947/3690345‑3690634

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