Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 407,069 A→G 15.4% intergenic (+39/‑164) glnE → / → rfaF Bifunctional glutamine synthetase adenylyltransferase/adenylyl‑removing enzyme/ADP‑heptose‑‑LPS heptosyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440407,0690AG15.4% 60.1 / 5.6 26intergenic (+39/‑164)glnE/rfaFBifunctional glutamine synthetase adenylyltransferase/adenylyl‑removing enzyme/ADP‑heptose‑‑LPS heptosyltransferase 2
Reads supporting (aligned to +/- strand):  ref base A (5/17);  new base G (2/2);  total (7/19)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGT  >  KT2440/406937‑407095
                                                                                                                                    |                          
ccTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAgg                          <  2:227929/135‑1 (MQ=255)
ccTGCAGAAGGAGGCCGGGGTGAGCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAgg                          <  2:76880/135‑1 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACGGGTACt                      >  4:104947/1‑134 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACGGGTACTg                     >  4:109419/1‑135 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACt                      >  6:80589/1‑134 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  2:243386/135‑1 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  1:241212/135‑1 (MQ=255)
     aGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGACGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  5:134578/135‑1 (MQ=255)
       aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACGGGTACTGca                   <  3:104947/135‑1 (MQ=255)
       aaGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACGGGTACTGca                   <  3:109419/135‑1 (MQ=255)
          gaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGACCGCCGTGAAGTGCGGCGGATCTGGGGTGAGTGGGGTTTAAGCGGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAg                <  6:134408/135‑1 (MQ=255)
          gaggCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAg                <  2:144482/135‑1 (MQ=255)
                ggggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa         <  7:135387/136‑1 (MQ=255)
                ggggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGGTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTg          <  4:127651/135‑1 (MQ=255)
                 gggTGTTGGTTGCGGGACAGTTTGCCGATGTTCGCGGTGAAGTGCGGCGGATGTGTGGTGAGTTGGGTTTAAGCTGATAGCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTTCTGCAGAGGCCGTGa         <  4:148009/135‑1 (MQ=255)
                 gggTGATCGATGCGGCACAGTTTGCCGATGAACTCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa         <  5:196806/135‑1 (MQ=255)
                 gggTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGa         <  5:80589/135‑1 (MQ=255)
                        cGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAgtgt  <  2:23947/135‑1 (MQ=255)
                                        gCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     >  2:78079/1‑100 (MQ=255)
                                        gCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  1:78079/100‑1 (MQ=255)
                                                  gCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     >  5:84554/1‑90 (MQ=255)
                                                  gCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  6:84554/90‑1 (MQ=255)
                                                                                               aCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     <  5:189056/45‑1 (MQ=255)
                                                                                               aCCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTg                     >  6:189056/1‑45 (MQ=255)
                                                                                                           tATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTg          >  5:57780/1‑44 (MQ=255)
                                                                                                           tATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTg          <  6:57780/44‑1 (MQ=255)
                                                                                                                                    |                          
CCTGCAGAAGGAGGCCGGGGTGATCGATGCGGCACAGTTTGCCGATGAACGCCGTGAAGTGCGGCGGATCTGGGGTGAGTTGGGTTTAAGCTGATACCGTGTTGGCTTATTCGCGGGCAAGCCCGCTCCCACAGGTACTGCACAGGCCCTGAGCAGTGT  >  KT2440/406937‑407095

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