Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 44,301 (T)9→8 intergenic (+132/‑359) oprD_1 → / → czcR_2 Porin D/Transcriptional activator protein CzcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB12244,2930T.93.8% 60.6 / ‑2.2 16intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (6/9);  total (7/9)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB122/44164‑44412
                                                                                                                                 |                                                                                                                       
cTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCtttttt                                                                                                                    >  6:161978/1‑135 (MQ=255)
    aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgc                                                                                                                >  5:64486/1‑133 (MQ=255)
     aTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgc                                                                                                                <  5:51013/134‑3 (MQ=255)
          cTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTaaa                                                                                                          <  6:64486/134‑1 (MQ=255)
                                         tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgca                                                                          >  6:93103/1‑135 (MQ=255)
                                         tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgca                                                                          <  6:7517/135‑1 (MQ=255)
                                               ttaaaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  5:159063/132‑1 (MQ=255)
                                                                            tGGCTTTAGTGCTACTTTGGTAAAATGTGGATATTTGGCTCCCATTGTGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGt                                       <  2:196730/135‑1 (MQ=255)
                                                                               cTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtctc                                    >  8:104695/1‑135 (MQ=255)
                                                                                   aGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTc                                 <  5:46782/134‑1 (MQ=255)
                                                                                       cTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAAtgt                            <  5:161978/135‑1 (MQ=255)
                                                                                          cTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTCCTAAAATTAAAGTAGGCGTTCGATTGTTCGCATCGCATGCCGGTGAAGTTCATTACTAAACCAATTACCAGTCTCTTAAAATTGAg                         >  1:62566/1‑135 (MQ=255)
                                                                                            ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTTAGCg                       >  5:60819/1‑135 (MQ=255)
                                                                                             tGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGa                      >  1:69630/1‑135 (MQ=255)
                                                                                                 aaaaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTc                  <  5:166646/135‑1 (MQ=255)
                                                                                                  aaaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACAcc                                                 >  7:166101/1‑103 (MQ=255)
                                                                                                  aaaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACAcc                                                 <  8:166101/103‑1 (MQ=255)
                                                                                                                 gcgcCCATTGGGCGCCTTTTTTTTGCCGAAAATTCAAGTAGGCGTTCGTTTGTACGCATCGCATGCCGGTGAAGTTCATGCCTACACCAATGACCAGTCGCTTCAATGTGAGCGCCCGCGCCGTCATCTGAGCg   >  4:91552/1‑134 (MQ=255)
                                                                                                                  cgcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg  <  3:170186/134‑1 (MQ=255)
                                                                                                                                 |                                                                                                                       
CTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB122/44164‑44412

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: