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

ATTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAGCG  >  KT2440/5298608‑5298871
                                                                                                                                   |                                                                                                                                     
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGGCCAGCGGCGCGTCgg                                                                                                                                   >  5:86381/1‑135 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCgg                                                                                                                                   <  5:287372/135‑1 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCg                                                                                                                                    >  4:16040/1‑134 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCg                                                                                                                                    <  3:16040/134‑1 (MQ=255)
aTTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCTGGCCCGGCGGCGCGTCgg                                                                                                                                   >  6:65124/1‑135 (MQ=255)
     tGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCaaa                                                                                                                              >  3:330973/1‑135 (MQ=255)
          aTAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGGGGGCCGGCGGGGCGGCGGGCCAATggcc                                                                                                                         >  3:274627/1‑132 (MQ=255)
            aGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCTGCGCGTCGGCCCAAAGCCAGc                                                                                                                       >  2:346450/1‑135 (MQ=255)
                 aaCAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCAGGCCCGGCTGCGCGGGGGCCCAATgcccgccccccc                                                                                                                 >  8:296169/1‑132 (MQ=255)
                    aGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCAGCTGCGCGTGGGCCAAAAGCCAcccccccccg                                                                                                               >  6:247588/1‑129 (MQ=255)
                     gTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGa                                                                                                              <  6:425113/135‑1 (MQ=255)
                         gCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGTGGCCCTGCTGCGCGGCGGCCAAAAGCCAGCCCGCCGCGCtgtg                                                                                                          >  4:303381/1‑135 (MQ=255)
                         gCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCCGGCGGGGCCGGCGGGGCGGGGGCCCAAATCCAGCCCCCCGCGCtgt                                                                                                           >  7:7139/1‑134 (MQ=255)
                           ggCCCTTCTGGAAGTACCGCGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGGCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATc                                                                                                       <  7:331236/135‑1 (MQ=255)
                              ccTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGGGGGCCCAGCTGCGCGGCGGCCCAATGCCAGCCCCCCGCGCTGTGATcgc                                                                                                     >  2:329288/1‑135 (MQ=255)
                                   tGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCCGCTGCGCGGCGGGCCAAATCCCGCCCCCCGCCCTGTGATCGGGCCACt                                                                                               >  8:301966/1‑136 (MQ=255)
                                     gAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCGGGCCCGGCTGCGCGGCGGCCAAATGCCCGCCCCCCCCCATGTGATCCCGCCACtc                                                                                              >  6:111142/1‑134 (MQ=255)
                                      aaGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAAc                                                                                                <  5:183303/132‑1 (MQ=255)
                                      aaGGACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGGCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAAc                                                                                                <  5:65124/132‑1 (MQ=255)
                                       aggaCCGAGGGCTGATCGGCTGCAGCGGGTGCTCCTGGGTCAGGGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  8:399046/133‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  3:322658/135‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGcc                                                                                           <  6:151548/135‑1 (MQ=255)
                                        gTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCTGCGCGTCCGCCCAAATCCACCCCCCCACGCTGTGATAGCGCAACAGGcc                                                                                           >  4:185998/1‑135 (MQ=255)
                                         tACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCg                                                                                          <  6:257656/135‑1 (MQ=255)
                                             gAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGCGGCCCAATGCCCGCCCGCAGCGATGTGGTGGCGCAACTTGCCgctac                                                                                      >  4:322658/1‑131 (MQ=255)
                                                        ggCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCAGCGGCGCGTCGGCCCAATGCCCGCCCGCCGCGATGTGATCGCGCCACGGGCCGATCAGCACCCCGGcc                                                                           >  2:308279/1‑135 (MQ=255)
                                                                cGGGTGCTCCTGGGGCAGTGCGTGCAACAGTTGCTCGCCCTGTCTCGCGGGCTGGCCCTGCTGCGCTTCGGCCAAATGCCAGCCCGCAGCGAGGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCg                                                                  <  8:192137/136‑1 (MQ=255)
                                                                       tCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGGGCGGCGGCCCCATGGCCGCCCGCCGCGATGGGCTAGCGCAACAGGGCGATCAGCCCCGCGGCCGGGGGGTCGcgtttg                                                            >  4:99327/1‑130 (MQ=255)
                                                                          tGGGGCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTc                                                          <  6:415377/134‑1 (MQ=255)
                                                                               cAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGAt                                                    <  5:22149/135‑1 (MQ=255)
                                                                                 gtgcgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACt                                                 <  8:7139/136‑1 (MQ=255)
                                                                                 gtgcgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCCCGCGCTGTGATCGCGCAACTGGCCGCTCAGCCCCGCGGCCGGGCGGGCGCTATTTTCGCGGGTACt                                                 >  2:172649/1‑135 (MQ=255)
                                                                                    cgtgcAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCAGCGCGTCGGCCCAATGCCCGCCCGCCGCGCTGTGATAGCGCAACTGGCCGATCAGCCCCGCGGCCGGGCGGTCGCTATTGTCGCGGATACTgcg                                              >  8:311682/1‑136 (MQ=255)
                                                                                         aaCAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGa                                         <  7:432982/136‑1 (MQ=255)
                                                                                            aGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGAcc                                       <  5:400829/135‑1 (MQ=255)
                                                                                            aGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGAcc                                       <  4:418092/135‑1 (MQ=255)
                                                                                            aGTTGCGCGGCCGGTCTCGCGGGCTGGCCCTGCTGCGCGGCGGCGAAATGCCAGCCGGCAGCGATGTGTTCGCGCTACTGGCCGATCAGCACCGAGGCCGGGCGGTCGCTATTGTCGGGGATACTGCGCCCGTccc                                      <  8:196086/136‑1 (MQ=255)
                                                                                               tGCTCGCCCGGTCTCGCGGGCGGGGCCTGCGGCGCGTCGGGCCAATGCCAGCCCGCCGCGCTGTGGTAGCGCGACGGGCCGCTGCGCACCCCGGCCGGGGGGTGGCGATTGTGGCGGGTGCTGGGCCCGCCCCAg                                    >  2:329733/1‑135 (MQ=255)
                                                                                                  tCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCt                                  <  6:293001/134‑1 (MQ=255)
                                                                                                        ggTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCGGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATAc                           <  6:86381/135‑1 (MQ=255)
                                                                                                            tcGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGtt                       <  1:161968/135‑1 (MQ=255)
                                                                                                               cgGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGt                    <  4:185172/135‑1 (MQ=255)
                                                                                                                   cTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  7:296169/136‑1 (MQ=255)
                                                                                                                    tGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  3:185998/135‑1 (MQ=255)
                                                                                                                    tGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGt               <  5:247588/135‑1 (MQ=255)
                                                                                                                       cccTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCCACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGCCCCAGCTGACATACAGTTGGTCACGTGca            >  4:220231/1‑134 (MQ=255)
                                                                                                                         ctgctgCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGAc         <  7:178714/136‑1 (MQ=255)
                                                                                                                          tgctgcGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGAc         <  4:81388/135‑1 (MQ=255)
                                                                                                                                gcgTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAgc   <  1:346450/135‑1 (MQ=255)
                                                                                                                                gcgTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTAATGGCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTa     <  4:3045/133‑1 (MQ=255)
                                                                                                                                gcgTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATGTCGCGGATACTGCGCCCGACCCAGCTGTCATACAGTTGGTCACGTGCTGACAGTAgcg  <  4:243971/135‑1 (MQ=255)
                                                                                                                                 cgTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAgcg  <  5:99548/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                     
ATTCATGGGCATAGCTGAACAGTGGGCTGCCCTTCTGGAAGTACCGAGGGCTGAACGGCTGCAGCGGGTGCTCCTGGGTCAGTGCGTGCAACAGTTGCTCGCCCGGTCTCGCGGGCTGGCCCTGCTGCGCGTCGGCCAAATGCCAGCCCGCAGCGATGTGATCGCGCAACTGGCCGATCAGCACCGAGGCCGGGCGTTCGCTATTGTCGCGGATACTGCGCCCGACCCAGCTGACATACAGTTGGTCACGTGCTGACAGTAGCG  >  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: