Predicted mutation
evidence seq id position mutation annotation gene description
RA KT2440 120,983 (T)9→8 intergenic (+280/+92) 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.98.0% 216.7 / ‑1.9 49intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (34/14);  minor base G (0/1);  total (34/15)
Fisher's exact test for biased strand distribution p-value = 3.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

TCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  KT2440/120872‑121105
                                                                                                       |                                                                                                                                  
tCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTACCCCCTCCCCATACGAAGTCCTGGGTTTTTTTTGCACTGAATTTGCCGGGAGCGGGt                                                                                                    >  6:151126/1‑135 (MQ=255)
          tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          >  5:93803/1‑135 (MQ=255)
          tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          <  5:207722/135‑1 (MQ=255)
            tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  1:92071/1‑135 (MQ=255)
             gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  4:92415/134‑1 (MQ=255)
             gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  8:280887/134‑1 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGCCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGGGATGt                                                                                 >  5:199561/1‑134 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGCACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGagggc                                                                                >  6:46135/1‑131 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCGGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  1:242310/1‑134 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  6:130257/1‑135 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  1:170510/1‑134 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  6:304542/1‑134 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 <  8:37711/134‑1 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  2:260304/1‑135 (MQ=255)
                    tAAGTGTTGGCAGATATCTGGACCCCGACCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGACCGGAATTTGCCGGGAGCGGGATGTCTCGCGAATGCGATGt                                                                                 >  8:61477/1‑134 (MQ=255)
                      aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  1:265504/1‑135 (MQ=255)
                      aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACATCCCCCTCCCCATACGAAGACCGGGGTTTTTTTTGACCGGTATTTGCCTGGAGCGGGTTTTCTAGCGAATGAGATGTCGa                                                                              >  7:134701/1‑135 (MQ=255)
                         gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                            >  3:212311/1‑134 (MQ=255)
                                     cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAtt                                                               >  6:129214/1‑135 (MQ=255)
                                     cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                <  3:211313/134‑1 (MQ=255)
                                              gCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                      <  2:260044/135‑1 (MQ=255)
                                                ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGTGCGGGTTTTCTCGCGAATGCGATGTCGATTACACCGCCGCCTTCGCGGGTCAAc                                                    >  2:129504/1‑135 (MQ=255)
                                                ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                    >  7:218989/1‑135 (MQ=255)
                                                ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                    >  5:89544/1‑135 (MQ=255)
                                                ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAAATTGACGGGAGCGGGTTTTCTCGCCAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                    >  1:262377/1‑135 (MQ=255)
                                                ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCAGACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTGTTCTCGCGAATGCGATGTCGAATACACCGCCGCAGTCGCGGGTAAAc                                                    >  8:221041/1‑135 (MQ=255)
                                                         ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAg                                           <  2:170510/135‑1 (MQ=255)
                                                          cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAgg                                          >  2:61840/1‑135 (MQ=255)
                                                          cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTAGAATACACCGCCGCATTCGCGGGTAAGACCGGGTCAgg                                          >  5:150709/1‑135 (MQ=255)
                                                          cGGACTTGGACCCCGAACTTCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACACCGCCCCCTTCGCGGGTAAACCCGGGTCAGgc                                         >  5:243209/1‑135 (MQ=255)
                                                           ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  7:11147/135‑1 (MQ=255)
                                                           ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         >  8:194333/1‑135 (MQ=255)
                                                           ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  1:223931/135‑1 (MQ=255)
                                                           ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  1:228548/135‑1 (MQ=255)
                                                            gACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         >  3:82849/1‑134 (MQ=255)
                                                             aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      <  1:232433/88‑1 (MQ=255)
                                                             aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      >  2:232433/1‑88 (MQ=255)
                                                                    ccccGAACTTCCCCCTCCCCATACCAAGTCCGGGGTTTTTTTTGCCCGGAATTCGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACCCCGCCGCAGTCGCGGGTAAACCCGGGTCGGGCGCAACCCa                                 >  5:144259/1‑134 (MQ=255)
                                                                     cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTACACCCGGGTCAGGCGCAcaccac                                >  7:160183/1‑134 (MQ=255)
                                                                     cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCAc                                >  5:257650/1‑134 (MQ=255)
                                                                     cccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCAc                                >  3:225593/1‑134 (MQ=255)
                                                                        gAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaaa                            <  8:218989/135‑1 (MQ=255)
                                                                           cgtCCCCCTCCCCATCCCATGCCGGGGGTTTTTTATGCCCGGAATTTGCCGGTAGCGGGTTTTCGCGCTAATGCTATGTCGAATGCACCGCCTCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAAc                           <  4:218232/131‑1 (MQ=255)
                                                                                cccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                    <  7:194333/135‑1 (MQ=255)
                                                                                 ccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCCGGCGCAAACCACTAAACGACCGGcg                   >  6:196136/1‑134 (MQ=255)
                                                                                       cATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTACACCCGGGTCAGGCGCAAACCACTACACGACCGGCTTGTCtt             >  4:173433/1‑135 (MQ=255)
                                                                                            gAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               >  8:6346/1‑128 (MQ=255)
                                                                                            gAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  7:6346/128‑1 (MQ=255)
                                                                                             aaGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGtt       >  6:250832/1‑135 (MQ=255)
                                                                                                  ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACAAAACGACCGGCTTGTCTTCCAGTtccat   <  6:164896/133‑1 (MQ=255)
                                                                                                   cggggTTTTTTTTGCCCGTAATTGGCCGGGTGGGGGTTTTCGCGCGAATGCTTTGTCTAATACACCGCGGCATTCGGGGGTAAACCCGGGTCAGGCGAAAACCACAAAACGACCGGCTTGTCTTCCAGTtccatc  <  4:72942/134‑1 (MQ=255)
                                                                                                   cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc  >  7:126076/2‑135 (MQ=255)
                                                                                                   cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccat   >  3:173825/2‑134 (MQ=255)
                                                                                                   cggggTTTTTTTTGCCCGCAATTTGCCGGGAGGGGGTTTTCTGGCGATTGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGAAAACAACTACACGACCGGCTTGTCTTCCAGTtccatc  <  8:229766/134‑1 (MQ=255)
                                                                                                       |                                                                                                                                  
TCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  KT2440/120872‑121105

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: