Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,432,242 A→C 25.9% E36A (GAG→GCG)  gpFI → Putative prophage major tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,432,2420AC25.9% 62.9 / 4.1 59E36A (GAG→GCG) gpFIPutative prophage major tail sheath protein
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base C (7/9);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCAAGGCGGTG  >  KT2440/3432108‑3432375
                                                                                                                                      |                                                                                                                                     
cgcgTCCCTTGTCTTATCTGGAGCCCAATTGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGGCCTCCTCGATCATTGGCGTGGGCGACACCTTCACCGa                                                                                                                                       <  1:337352/135‑1 (MQ=255)
cgcgGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGa                                                                                                                                       <  2:173590/135‑1 (MQ=255)
   ggCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACCCCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCGCGATCATTTGCTTGGTCGCCACCTTCCCCGCGGGc                                                                                                                                   >  8:101079/1‑136 (MQ=255)
     cccTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCCATGCGCTGCCGTCCGCCGCGATCCTTGGCTTGGTCGCCACCTTCACCCCGGGCgg                                                                                                                                 >  8:114696/1‑136 (MQ=255)
        ttGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACAACGTCACCGACGGCGGTAc                                                                                                                              >  8:365025/1‑136 (MQ=255)
                tCTGGAGCCCACATGAATGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACg                                                                                                                       <  8:104770/135‑1 (MQ=255)
                 cTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCGACg                                                                                                                       <  1:347626/134‑1 (MQ=255)
                  tGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACCCCGGGGCGCGCAGCATTGCGCTGCCGGCCTCCTCGCTCATTGGCTTGGTCGACCCCGTCACCCAGGGCGGGACCCCCCAGGcc                                                                                                                    >  7:61346/1‑136 (MQ=255)
                  tGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGcc                                                                                                                    >  8:353700/1‑136 (MQ=255)
                  tGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCAATTCGCGGCCGTCCACCGCGATCATTGGCTTGGTCCCCACCATCACCGACGGCGGTACACCCACGGc                                                                                                                     >  5:163580/1‑135 (MQ=255)
                   gccccccACATGACTGGATTCTTTCACGGCGTTACCGAAACGAACGTCGCCACCGGGGCGCGCAGCATTGCGCTGCCGGCCGCCTCGCTCATTGGCTTGGGCGACACCTTCACCCAGGGCGGTACCGCCACGGcc                                                                                                                    <  4:47928/131‑1 (MQ=255)
                       cccACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCCCGGCCAAGt                                                                                                                >  4:96028/1‑135 (MQ=255)
                       cccACATGACTGGATTCTTTAACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGAGCTGCCGTGCTCCTCGCTCATTGGGTTTGTCGACACCATAACCCCGGGCCGTAACGCCACCGCCAACTa                                                                                                               >  8:12488/1‑136 (MQ=255)
                       cccACATGACTGGATGCTTGCCCGGCGTTACCGTAACGACCGTCGACACCGGGGCGCGCAGCATTGCGCGGCCGTCCTCCTCGATCATTGGCTTGGGCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTa                                                                                                               <  8:323697/136‑1 (MQ=255)
                        cccaaTGACTGGATTCTTTCACGGCGTTACCGTAACCACCGTCGACACCGGGGCGCGAAGCTTGGCGCGGCCGTCCTCCGCGATCATTGGCTTGGGCGACACCTGCCCCGCGGGCGGTACCGCCACGGCCAAGTa                                                                                                               <  6:85459/131‑1 (MQ=255)
                        ccACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTTCGCTGGCGGCCTGCGCGATCATTGGGTTTGTGGACACCTTCACCGAGGGCGGGTCCCCCCCCGCCaacg                                                                                                                >  6:353470/1‑132 (MQ=255)
                        ccACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTTCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGCCACCTTCACCCCGGGCGGTTCCGCCACCGCCAAGt                                                                                                                >  3:401556/1‑134 (MQ=255)
                        ccACATGACTGGAGTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGCCCTCCTCGATCATGGGCTGGGGCGACACCTGCCCCGAGGGCGGTACCGCCACGGCCAAGTa                                                                                                               <  4:227853/135‑1 (MQ=255)
                        ccAAATGAGTGGATTCTTACAGGGCGTTCCCTAAACGACCGTCGACACCGGGGCGCGCAGATTTGCGCGGCCGTCCGCCTCGATCATTGGCTGGGTCGACACCTTCACCGCGGGCGGAACCCCCACGGCCAAGTa                                                                                                               <  4:201763/135‑1 (MQ=255)
                           caTGCCTGGATTCTTTCACGGCGTTACCGAAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACaa                                                                                                            <  4:85049/135‑1 (MQ=255)
                            aTGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACg                                                                                                          <  7:365025/136‑1 (MQ=255)
                               aCTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCCGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACCCCTTCCCCGCAGGCGGTACCGCCACCGCCAACGACAACGccc                                                                                                       >  8:84759/1‑136 (MQ=255)
                                   gATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGCGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGt                                                                                                   <  7:353700/136‑1 (MQ=255)
                                   gATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTa                                                                                                     <  7:230394/134‑1 (MQ=255)
                                    aTTCTTCCCCGGCGTCCCCGTACCCACGGTCGACCCCGGGGCGCGCAGCATTGCGCGGCCGTCCGCCTCGATCTTTGGCTTGGTCGACACCTTCCCCGCGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGt                                                                                                   <  5:266492/135‑1 (MQ=255)
                                      tCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGc                                                                                                 <  2:140658/135‑1 (MQ=255)
                                      tCTGTCACGGCGTTACCTTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGc                                                                                                 <  8:439918/135‑1 (MQ=255)
                                       cTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCCTTGCGCTGGCCTGCACCGCGCTAATTGGGGTGGTCGCCCCCTTCACCCCGGGCGGTACCCCCCCGGCCCAGTACAACGCCCAAGTggt                                                                                                >  2:142173/1‑133 (MQ=255)
                                                             aacgTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTGCTCCTCGATAATTGGCTTGGTCGCCACCTTCACCGAGGGCGGTACCCCCACGGCCAAGTACAACGACCTAGTGCTGATCACAAACGAGCGCGAGGCg                                                                          >  5:296725/1‑135 (MQ=255)
                                                             aacgTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGCTCATTGGCGTGGTCGACACCTTCACCGAGGGCCGGCCCGCCACGGCCAAGTACAACCACCTAGTGCGGATCACCAACGAGCGCGAGGCg                                                                          >  1:257274/1‑135 (MQ=255)
                                                             aacgTCGACACCGGGGCGCGCAGCAATGCGCTGCCGTCCTACGCGATCAATGGCTTGGGCGCCACCCTCACCGCGGGCGGTAACGCCACGGCCAACTACAACGACCTAGTGCTGATCAACAACGAGGGCCAGGCg                                                                          >  4:313053/1‑135 (MQ=255)
                                                                  cGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGGCGACACCTTACCCGAGGGCGGTACCGCCCCGGCCAAGTACAAAGACCTAGTGCTGATCACCAACGCGCGCGCGGCGGTcgcc                                                                    >  7:51230/1‑136 (MQ=255)
                                                                            gcgcgcAGCTTTGCGCGGCCGTCCTCCTCGATCATGGGCTTGGTCCACACCTTCCCCGCGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGt                                                           <  2:341046/135‑1 (MQ=255)
                                                                             cgcgcAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCCCCGAGGGCGGTTCCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGt                                                           >  3:410598/1‑134 (MQ=255)
                                                                               cgccgcATTGCGCGGCCGTCCTCCTCGATCATGGGCGGGGCCCACCCCTTCCCCGCGGGCGGTACCCCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGaaa                                                       <  7:77036/136‑1 (MQ=255)
                                                                                              cgtgcctccgCGATCATTGGCTTGGCCGACACCTTCACCGAGGGCGTTACCCCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCAcc                                         <  5:111373/131‑1 (MQ=255)
                                                                                              ccGGCCGCCGCGATCATTGGCGTGGGCGACACCTTCCCCGTGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCAcc                                         <  6:247309/135‑1 (MQ=255)
                                                                                               cGTCCTCCTCGATCATGGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCAcc                                         <  2:332866/134‑1 (MQ=255)
                                                                                                 tcctccTCGATCATTGGCTTGGTCGACACCTTCACCGa                                                                                                                                       <  4:188849/38‑1 (MQ=255)
                                                                                                 tcctccTCGATCATTGGCTTGGTCGACACCTTCACCGa                                                                                                                                       >  3:188849/1‑38 (MQ=255)
                                                                                                 tcctccTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAgg                                     >  7:29079/1‑136 (MQ=255)
                                                                                                 tcctccTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGGGAAACCTCGGCGATCACCAAgg                                     >  7:187046/1‑136 (MQ=255)
                                                                                                   ctcctcGATCATGGGCTTGGGCGACACCTTCACCGCGGGCGGTCCCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGc                                    <  1:38389/135‑1 (MQ=255)
                                                                                                     ccgcGATCATTGGCTGGGCCCCCACCTCCCCCGCGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGt                                  <  5:442067/132‑1 (MQ=255)
                                                                                                      cgcGATCTTTGGCTTGGTCGACCCCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGc                                <  7:233717/134‑1 (MQ=255)
                                                                                                       tcGATCATTGGCTTGGTCGACACCTTCACCGCGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCa                               <  8:337572/136‑1 (MQ=255)
                                                                                                         gctCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCa                               <  7:362510/132‑1 (MQ=255)
                                                                                                          agcATTGGCTGGGGCGACACCTGCACCGAGGGCGGCACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCa                               <  6:161738/131‑1 (MQ=255)
                                                                                                          aTCCTTGGCTGGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCa                               <  4:94713/133‑1 (MQ=255)
                                                                                                          aTCATTGGCTGGGACCACACATCCCCCGGGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCa                               <  5:353470/133‑1 (MQ=255)
                                                                                                              ttGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCa                         >  4:318227/1‑135 (MQ=255)
                                                                                                                  cTTGGTCGACACCTTCACCGAGGCCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCt                      >  4:375146/1‑134 (MQ=255)
                                                                                                                   ttGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTaca                   >  7:119106/1‑136 (MQ=255)
                                                                                                                        cGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGCTCACCCAGGCGTGCAAGGCCATCAACACCcgc              >  7:383613/1‑136 (MQ=255)
                                                                                                                        cGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGc                                    <  4:331974/114‑1 (MQ=255)
                                                                                                                        cGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGc                                    >  3:331974/1‑114 (MQ=255)
                                                                                                                          acacCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGCCCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCcgcg             >  4:338272/1‑135 (MQ=255)
                                                                                                                          acacCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTacac                  <  2:350546/130‑1 (MQ=255)
                                                                                                                          acacCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTacac                  >  1:350546/1‑130 (MQ=255)
                                                                                                                              cTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCaa         >  3:391130/1‑135 (MQ=255)
                                                                                                                                 ccccGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCAAGGc      <  8:354117/133‑1 (MQ=255)
                                                                                                                                    cGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCAATGCgg    >  6:253878/1‑134 (MQ=255)
                                                                                                                                    cGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCAAGGCGGTg  >  8:322183/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                     
CGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGGTACCGCCACGGCCAAGTACAACGACCTAGTGCTGATCACCAACGAGCGCGAGGCGGTCGCCGCGTTCGGTGAAACCTCGGCGATCACCAAGGCGTGCAAGGCCATCTACACCCGCGCCAAGGCGGTG  >  KT2440/3432108‑3432375

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