Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 5,298,739 T→G 17.8% D780A (GAC→GCC)  recC ← RecBCD enzyme subunit RecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24405,298,7390TG17.8% 90.1 / 2.1 52D780A (GAC→GCC) recCRecBCD enzyme subunit RecC
Reads supporting (aligned to +/- strand):  ref base T (12/29);  new base G (9/2);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 3.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAGCG  >  KT2440/5298608‑5298871
                                                                                                                                   |                                                                                                                                     
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGGCCAGCGGCGCGTC‑gg                                                                                                                                   >  5:86381/1‑135 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑gg                                                                                                                                   <  5:287372/135‑1 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑g                                                                                                                                    >  4:16040/1‑134 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑g                                                                                                                                    <  3:16040/134‑1 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCTGGCCCGGCGGCGCGTC‑gg                                                                                                                                   >  6:65124/1‑135 (MQ=255)
     tGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCaaa                                                                                                                              >  3:330973/1‑135 (MQ=255)
          aTAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGGGGGCCGGCGGGGCGGC‑GGGCCAATggcc                                                                                                                         >  3:274627/1‑132 (MQ=255)
            aGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCTGCGCGTC‑GGCCCAAAGCCAGc                                                                                                                       >  2:346450/1‑135 (MQ=255)
                 aaCAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCAGGCCCGGCTGCGCGGG‑GGCCCAATgcccgccccccc                                                                                                                 >  8:296169/1‑132 (MQ=255)
                    aGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCAGCTGCGCGTG‑GGCCAAAAGCCAcccccccccg                                                                                                               >  6:247588/1‑129 (MQ=255)
                     gTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGa                                                                                                              <  6:425113/135‑1 (MQ=255)
                         gCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGTGGCCCTGCTGCGCGGC‑GGCCAAAAGCCAGCCCGCCGCGCtgtg                                                                                                          >  4:303381/1‑135 (MQ=255)
                         gCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCGGGGCCGGCGGGGCGGG‑GGCCCAAATCCAGCCCCCCGCGCtgt                                                                                                           >  7:7139/1‑134 (MQ=255)
                           ggCCCTTCTGGAAGTACCGCGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGGCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATc                                                                                                       <  7:331236/135‑1 (MQ=255)
                              ccTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGGGGCCCAGCTGCGCGGC‑GGCCCAATGCCAGCCCCCCGCGCTGTGATcgc                                                                                                     >  2:329288/1‑135 (MQ=255)
                                   tGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCCGCTGCGCGGCGGGCCAAAT‑CCCGCCCCCCGCCCTGTGATCGGGCCACt                                                                                               >  8:301966/1‑136 (MQ=255)
                                     gAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCGGCTGCGCGGC‑GGCCAAATGCCCGCCCCCCCCCATGTGATCCCGCCACtc                                                                                              >  6:111142/1‑134 (MQ=255)
                                      aaGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAAc                                                                                                <  5:183303/132‑1 (MQ=255)
                                      aaGGACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGGCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAAc                                                                                                <  5:65124/132‑1 (MQ=255)
                                       aggaCCGAGGGCTGATCGGCTGCAGCGGGTGCTCCTGGGTCAGGGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  8:399046/133‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  3:322658/135‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  6:151548/135‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCTGCGCGTC‑CGCCCAAATCCACCCCCCCACGCTGTGATAGCGCAACAGGcc                                                                                           >  4:185998/1‑135 (MQ=255)
                                         tACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCg                                                                                          <  6:257656/135‑1 (MQ=255)
                                             gAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGC‑GGCCCAATGCCCGCCCGCAGCGATGTGGTGGCGCAACTTGCCgctac                                                                                      >  4:322658/1‑131 (MQ=255)
                                                        ggCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCGGCGCGTC‑GGCCCAATGCCCGCCCGCCGCGATGTGATCGCGCCACGGGCCGATCAGCACCCCGGcc                                                                           >  2:308279/1‑135 (MQ=255)
                                                                cGGGTGCTCCTGGGGCAGTGCGTGCAACAGTTGCTCGCCCTGTCTCGCGGGCTGGCCCTGCTGCGCTTC‑GGCCAAATGCCAGCCCGCAGCGAGGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCg                                                                  <  8:192137/136‑1 (MQ=255)
                                                                       tCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGGGCGGC‑GGCCCCATGGCCGCCCGCCGCGATGGGCTAGCGCAACAGGGCGATCAGCCCCGCGGCCGGGGGGTCGcgtttg                                                            >  4:99327/1‑130 (MQ=255)
                                                                          tGGGGCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTc                                                          <  6:415377/134‑1 (MQ=255)
                                                                               cAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGAt                                                    <  5:22149/135‑1 (MQ=255)
                                                                                 gtgcgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACt                                                 <  8:7139/136‑1 (MQ=255)
                                                                                 gtgcgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCC‑CCGCGCTGTGATCGCGCAACTGGCCGCTCAGCCCCGCGGCCGGGCGGGCGCTATTTTCGCGGGTACt                                                 >  2:172649/1‑135 (MQ=255)
                                                                                    cgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCAGCGCGTC‑GGCCCAATGCCCGCCCGCCGCGCTGTGATAGCGCAACTGGCCGATCAGCCCCGCGGCCGGGCGGTCGCTATTGTCGCGGATACTgcg                                              >  8:311682/1‑136 (MQ=255)
                                                                                         aaCAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGa                                         <  7:432982/136‑1 (MQ=255)
                                                                                            aGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGAcc                                       <  5:400829/135‑1 (MQ=255)
                                                                                            aGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGAcc                                       <  4:418092/135‑1 (MQ=255)
                                                                                            aGTTGCGCGGCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGC‑GGCGAAATGCCAGCCGGCAGCGATGTGTTCGCGCTACTGGCCGATCAGCACCGAGGCCGGGCGGTCGCTATTGTCGGGGATACTGCGCCCGTccc                                      <  8:196086/136‑1 (MQ=255)
                                                                                               tGCTCGCCCGGTCTCGCGGGCGGGGCCTGCGGCGCGTC‑GGGCCAATGCCAGCCCGCCGCGCTGTGGTAGCGCGACGGGCCGCTGCGCACCCCGGCCGGGGGGTGGCGATTGTGGCGGGTGCTGGGCCCGCCCCAg                                    >  2:329733/1‑135 (MQ=255)
                                                                                                  tCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCt                                  <  6:293001/134‑1 (MQ=255)
                                                                                                        ggTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCGGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATAc                           <  6:86381/135‑1 (MQ=255)
                                                                                                            tcGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGtt                       <  1:161968/135‑1 (MQ=255)
                                                                                                               cgGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGt                    <  4:185172/135‑1 (MQ=255)
                                                                                                                   cTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  7:296169/136‑1 (MQ=255)
                                                                                                                    tGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  3:185998/135‑1 (MQ=255)
                                                                                                                    tGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  5:247588/135‑1 (MQ=255)
                                                                                                                       cccTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCCACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGCCCCAGCTGACATACAGTTGGTCACGTGca            >  4:220231/1‑134 (MQ=255)
                                                                                                                         ctgctgCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGAc         <  7:178714/136‑1 (MQ=255)
                                                                                                                          tgctgcGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGAc         <  4:81388/135‑1 (MQ=255)
                                                                                                                                gcgTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAgc   <  1:346450/135‑1 (MQ=255)
                                                                                                                                gcgTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTAATGGCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTa     <  4:3045/133‑1 (MQ=255)
                                                                                                                                gcgTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTA‑TGTCGCGGATACTGCGCCCGACCCAGCTGTCATACAGTTGGTCACGTGCTGACAGTAgcg  <  4:243971/135‑1 (MQ=255)
                                                                                                                                 cgTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAgcg  <  5:99548/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
ATTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTC‑GGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAGCG  >  KT2440/5298608‑5298871

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑