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.88.2% 113.6 / 12.0 34intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (12/18);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCCAGC  >  ACB122/44168‑44424
                                                                                                                             |                                                                                                                                   
aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgcc                                                                                                                           >  7:285671/1‑133 (MQ=255)
aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgcc                                                                                                                           >  3:145003/1‑133 (MQ=255)
              gCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                            <  4:67037/136‑1 (MQ=255)
              gCCAGATTTCCATTGGCGGGTCCTGTTAGTATACCGAGGATTACTCATTCACCTTGTGTTGCTTAAGTGCTCCTTTGTTAATAGGTGGATATTTGGCGCCGTTTGGGCGCCTTTTTTTTGCCTAATTTTCAAGTAg                                                                                                            <  4:257491/136‑1 (MQ=255)
                          ttGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGt                                                                                                <  5:97465/136‑1 (MQ=255)
                                  tCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                                        <  4:145003/136‑1 (MQ=255)
                                    cTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgc                                                                                       >  1:69226/1‑135 (MQ=255)
                                     tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgcat                                                                                     >  7:94297/1‑136 (MQ=255)
                                       ttAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGc                                                                                   >  3:193433/1‑136 (MQ=255)
                                           tAAAACGAGCATTACGCATTCATCTTGAGTGGCTTAAGTGCTCCTTTGGTAATAGGTGGATATTTGGCGCCCATTGGTCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCgg                                                                                <  1:92685/135‑1 (MQ=255)
                                           tAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                                 <  1:162327/134‑1 (MQ=255)
                                               aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCGCCTTTGGAAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                            <  2:219003/135‑1 (MQ=255)
                                                 gAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                                         <  5:250656/136‑1 (MQ=255)
                                                      ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                    <  3:307699/136‑1 (MQ=255)
                                                      ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                    <  8:285671/136‑1 (MQ=255)
                                                      ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                    <  8:288894/136‑1 (MQ=255)
                                                       tACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                                         >  2:322246/1‑130 (MQ=255)
                                                       tACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                                         <  1:322246/130‑1 (MQ=255)
                                                            aTTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGGCTTAAATTGAAGTATGCGTTGGCTTGTACGCCTCGCATGGCGGGGAAGTTGATGACTacac                                                              >  5:6200/1‑136 (MQ=255)
                                                                     gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                       <  5:198384/134‑1 (MQ=255)
                                                                      aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                       <  2:215881/133‑1 (MQ=255)
                                                                         ggCTAAAGTGCTCCTTTGGTAAAAGGTGGATATTAGGCTACCATTCCTCGCCATTTTTATGCCTAAAAATAAAGTAGCAGTTACATTATACGCAAAGCATGACCGTGAAGAAGAAGACTAAACCAAAGAACaata                                                  >  8:99446/1‑132 (MQ=255)
                                                                           cTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTTGCGCCCATTGGGCGCCTTTTTTTTTCCGAAAAATAAAGTAGGCGTTCGATTTTAAGCCTCGCAAGCCGGGGAAGTTAATGACGACCCCAAAGACCacgcgag                                               >  8:185164/1‑130 (MQ=255)
                                                                              aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTc                                             >  8:336734/1‑135 (MQ=255)
                                                                              aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGACCATTGGGCGCCTTTTTTTTTCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTc                                             <  8:285912/135‑1 (MQ=255)
                                                                                     ccTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGCCTAAAAATAAATTAGGCGTTGGATTGTACGCCACGCATGCCGGTGAACTTAATTACTACCCCCATGACCAGTCTCTTCAATTTGAg                                     >  6:200027/1‑136 (MQ=255)
                                                                                        ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGc                                                                                                                           <  5:154145/48‑1 (MQ=255)
                                                                                        ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGc                                                                                                                           >  6:154145/1‑48 (MQ=255)
                                                                                        ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTAGATTTTAAGCATCGAATGCCGGTTAAGTTCATGACTACACCAATTACCAGTCTCTTAAATGTGAGCGa                                  >  3:192717/1‑136 (MQ=255)
                                                                                                       tatTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt                    <  2:69226/135‑1 (MQ=255)
                                                                                                        atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATc                     >  1:268113/1‑133 (MQ=255)
                                                                                                        atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt                    >  3:116201/1‑134 (MQ=255)
                                                                                                             gcgcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg              <  3:279386/135‑1 (MQ=255)
                                                                                                                cccATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGGTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCGTCAATGTGAGCGCCCTCGCCGTCATCGGAGCGGCACg          >  4:45378/1‑136 (MQ=255)
                                                                                                                 ccATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGa         <  7:11958/136‑1 (MQ=255)
                                                                                                                        gcgcCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCCAGc  <  3:161836/136‑1 (MQ=255)
                                                                                                                             |                                                                                                                                   
AATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCCAGC  >  ACB122/44168‑44424

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: