Predicted mutation
evidence seq id position mutation freq annotation gene description
RA pSIJ663_pHT315_araE_xy 5,612 T→A 71.5% intergenic (+881/–) ampR → /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*pSIJ663_pHT315_araE_xy5,6120TA71.5% 29.5 / 9.5 21intergenic (+881/–)ampR/––/–
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base A (9/6);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGGAGCCTATGGAAAAACGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTATGACACAAAAAAGCGCAGCTATTTCAGCTGCGCTTTTTTCGAAACAAACGCATTTGACCAAACAAGCAGACAAAAAAGGCGATGGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACC  >  pSIJ663_pHT315_araE_xy/5476‑5745
                                                                                                                                        |                                                                                                                                     
cgaagccTATGGAAAAACGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTA                                                                                                                                     >  2:263830‑M1/4‑139 (MQ=255)
cgaagccTATGGAAAAACGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTA                                                                                                                                     <  6:93295‑M1/136‑1 (MQ=255)
             aaaaAGCCAGGAAGGCTGACTTTTTAGTGTTAGTTGCCTTTTGATGGCCTTTTGCTCACATGTTGTTTCCTGCGTTTTCGCATGATTCTGTGGATAAAGGTTTTACCGCATTTGATTGATCTGTTATGACACaataatg                                                                                                                        <  8:132058/139‑3 (MQ=255)
                 cGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTATGACACAAAAAAgcg                                                                                                                      >  4:401897/1‑137 (MQ=255)
                                     tACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagt                                                                                                 <  6:223250‑M1/138‑37 (MQ=255)
                                                ggcttttgctggccttttgctcaCATGTTCTTTCCTGCGTTATCCCCTTATTATGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgcctcagccgaacgaccgagcgcagcgagtcagtgagcgagg                                                                                     <  2:31225‑M1/137‑49 (MQ=255)
                                                                ttgctCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaccgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcc                                                                      >  2:235153‑M1/1‑75 (MQ=255)
                                                                        aTGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgca                                                             <  4:413644‑M1/139‑73 (MQ=255)
                                                                                  cTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTATGACACAAAAAAGCGCAGCTATTTCAGCTGCGCTTTTTTCGAAACAAACGCATTTGACCAAACAAGCAGACAAAAAAGGCGa                                                   <  1:263830/139‑1 (MQ=255)
                                                                                  cTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgcccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaacccgcctct                                                   >  8:140298‑M1/1‑57 (MQ=255)
                                                                                                    cTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccga                                  <  5:186671‑M1/138‑100 (MQ=255)
                                                                                                     tgtgGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgagg                                                                                     >  5:314254‑M1/1‑38 (MQ=255)
                                                                                                     tgtgGATAACCGTATTACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgagg                                                                                     <  6:314254‑M1/86‑49 (MQ=255)
                                                                                                       tgGATAACCGTATTACCGCCTTTGAGTGAGCTGTTATGACACAAAAAAGCGCAGCTATTTCAGCTGCGCTTTTTTCGAAACAAACGCATTTGACCAAACAAGCAGACAAAAAAGGCGATGGGATAGATCACCATTTCa                               >  1:110992/1‑138 (MQ=255)
                                                                                                                   ttACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgccca                                                                    >  8:21184‑M1/1‑24 (MQ=255)
                                                                                                                   ttACCGCCTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgccca                                                                    <  7:21184‑M1/89‑66 (MQ=255)
                                                                                                                         ccTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccgattcattaatgcagctggcacga            >  6:178218‑M1/1‑18 (MQ=255)
                                                                                                                         ccTTTGAGTGAGCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccgattcattaatgcagctggcacga            >  1:107955‑M1/1‑18 (MQ=255)
                                                                                                                           tttGAGTGAGCTGATAccgctcgccgcagccgaacgaccgcgcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccgattcattaatgcagctggcacgaca          >  4:247448‑M1/1‑16 (MQ=255)
                                                                                                                                   agCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccgattcattaatgcagctggcacgacaggtttccc  >  5:258800‑M1/1‑8 (MQ=255)
                                                                                                                                   agCTGATAccgctcgccgcagccgaacgaccgagcgcagcgagtcagtgagcgaggaagcggaagagcgcccaatacgcaaaccgcctctccccgcgcgttggccgattcattaatgcagctggcacgacaggtttcac  >  7:147542‑M1/1‑8 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CGGAGCCTATGGAAAAACGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGTTATGACACAAAAAAGCGCAGCTATTTCAGCTGCGCTTTTTTCGAAACAAACGCATTTGACCAAACAAGCAGACAAAAAAGGCGATGGGATAGATCACCATTTCATAAGTATTGATAGGTGCGTTTAAAACACC  >  pSIJ663_pHT315_araE_xy/5476‑5745

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: