Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 120,975 Δ1 bp 100% intergenic (+272/+100) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440120,9750T.94.7% 239.1 / 2.4 57intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (33/21);  total (33/24)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCAT  >  KT2440/120853‑121104
                                                                                                                          |                                                                                                                                 
gAATGGCTGGTGAAGCCCTTCGATGCTCCTATGGTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGgttgttttgccc                                                                                                                        >  7:148926/1‑130 (MQ=255)
 aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCg                                                                                                                      >  7:64974/1‑134 (MQ=255)
 aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCg                                                                                                                      >  7:318130/1‑134 (MQ=255)
 aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTCGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCg                                                                                                                      >  3:271850/1‑134 (MQ=255)
    ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCGCCCCCTACGAAGTCCGCGGTGTTGTTTGCCCGGAAt                                                                                                                  >  1:300526/1‑135 (MQ=255)
      cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAttt                                                                                                                >  1:254705/1‑135 (MQ=255)
                        gCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                              <  5:303459/135‑1 (MQ=255)
                            cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCAAACGAACTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgc                                                                                           >  2:36845/1‑134 (MQ=255)
                               tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       >  3:173265/1‑135 (MQ=255)
                                gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  8:64974/134‑1 (MQ=255)
                                gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  3:217451/134‑1 (MQ=255)
                                gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  5:274206/134‑1 (MQ=255)
                                       tAAGTGTTGGTAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAt                                                                                                                  <  2:92672/100‑1 (MQ=255)
                                       tAAGTGTTGGTAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAt                                                                                                                  >  1:92672/1‑100 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCTGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                               >  4:81933/1‑135 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  5:301231/127‑1 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       >  6:301231/1‑127 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                >  1:178311/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                >  5:127567/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                >  4:157705/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                               >  3:173678/1‑135 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGGTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATgg                                                                                >  6:221635/1‑133 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGGTTTGTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                               >  3:49871/1‑135 (MQ=255)
                                         aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCTCGAATGCGATGTCGa                                                                             >  7:227012/1‑135 (MQ=255)
                                              tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       >  8:185507/1‑120 (MQ=255)
                                              tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  7:185507/120‑1 (MQ=255)
                                              tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATaca                                                                        >  5:323838/1‑135 (MQ=255)
                                                      atCTGGACCCCGCCCTAGCGGTCCTGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCa                                                                <  7:236476/135‑1 (MQ=255)
                                                       gcgggACCCCGCCCTAGCGGTCCGGACTTGGCCCCCGAACTCCCCCCTCCCCATCCGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                               <  4:85014/132‑1 (MQ=255)
                                                        cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGTGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAtt                                                              >  5:307482/1‑135 (MQ=255)
                                                            cccccGCCCTCGCGGGCCGGACGTGGACCCCGAACTCCCCCCTCCCCCTACGACGGCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCGCGCGAATGCGATGTCGAATCCACCGCCGCATTcgcg                                                          <  4:93509/134‑1 (MQ=255)
                                                             ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCgg                                                         <  4:173265/135‑1 (MQ=255)
                                                                 gCCCTAGCGGTCCGGACTTGGACCCCGGACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                     >  3:50467/1‑135 (MQ=255)
                                                                 gCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                     <  1:36845/135‑1 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGATTTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTacac                                                   >  8:148272/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                   >  7:336355/1‑135 (MQ=255)
                                                                    cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGGTGGGGAATACACCGCCGCATTCGCGGGTAAAcc                                                  >  5:197189/1‑135 (MQ=39)
                                                                          gTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCggg                                              >  5:62895/1‑133 (MQ=255)
                                                                            ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAg                                          <  6:307482/135‑1 (MQ=255)
                                                                             cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATAAACCTCCGCATTCGCGGGTAAACCCGGGTCAgg                                         >  3:12111/1‑135 (MQ=255)
                                                                              ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGGTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGACTACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                        <  6:4271/135‑1 (MQ=255)
                                                                              ggACTTGGACCCCGAACTTCCCCCACCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                        <  5:239475/135‑1 (MQ=255)
                                                                              ggACATGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                        <  1:244778/135‑1 (MQ=255)
                                                                                aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                     >  7:29476/1‑88 (MQ=255)
                                                                                aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                     <  8:29476/88‑1 (MQ=255)
                                                                                       ccccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTAGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCCAACCa                                >  4:283413/1‑134 (MQ=255)
                                                                                        cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  1:57682/78‑1 (MQ=255)
                                                                                        cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       >  2:57682/1‑78 (MQ=255)
                                                                                        cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCAc                               >  4:299205/1‑134 (MQ=255)
                                                                                             aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       <  2:329287/73‑1 (MQ=255)
                                                                                             aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                       >  1:329287/1‑73 (MQ=255)
                                                                                             aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACg                         >  3:109801/1‑135 (MQ=255)
                                                                                                ccccccACCCCATACGAGGGCGGGGTTTTTTTTTGCCCGTATTTTGCGGGGAGGGGTTTTTCGCGCTATTGCGATGTCTAAACCCCCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAcc                      <  2:254705/134‑1 (MQ=255)
                                                                                                 cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCg                     >  5:164050/1‑135 (MQ=255)
                                                                                                   cccTCCCCATCCGAAGGCGGGGGTTTTTTTTGCCCGTATTTGGCGGGGTGCGGGTTTTCGCGCTAATGCTTGGTCTAATACACCGCCGCATGCGCGGTTAACCCGGGGTCAGGCGAAAACCACTAAACGACCGGc                   <  6:241386/135‑1 (MQ=255)
                                                                                                   ccatCCCAATACTAAGCCGGGGTTTTTTTTTGCGCGTATTTTGGCGGGCGCGGGTTTTCTCGCTACTGCTATGTCTAATCCACCGCCGCATTCGCGGTTAAACCGGGGTCAGGCGAAAACCACTAAACGACCGGc                   <  6:62895/132‑1 (MQ=255)
                                                                                                        cccATCCGAAGTCCGGGGTTTTTTTTGCCCGGATTTGGCCGGGAGCGGGTTTTCTCGCGAATGCTTTGTCTAATACACCGCCGCATTCGGGGGTAAACCGGGGTCAGGCGCAAACCACTACACGACCGGCTTGTc              <  6:127567/135‑1 (MQ=255)
                                                                                                                aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGtt      <  3:299205/135‑1 (MQ=255)
                                                                                                                      cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccat  >  8:213007/2‑134 (MQ=255)
                                                                                                                          |                                                                                                                                 
GAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCAT  >  KT2440/120853‑121104

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: