Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 44,293 Δ1 bp 100% intergenic (+124/‑367) 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.95.2% 172.9 / ‑1.1 41intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (8/31);  total (10/31)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

TTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB122/44162‑44412
                                                                                                                                   |                                                                                                                       
ttCTGTAAATCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAGGGTGGATATTTGGCGCCCATTGGGCGCCtttt                                                                                                                      <  1:124573/135‑1 (MQ=255)
      aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgc                                                                                                                >  8:345307/1‑133 (MQ=255)
      aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttttttttgc                                                                                                                >  3:295867/1‑133 (MQ=255)
                cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAAtt                                                                                                       >  3:212461/1‑133 (MQ=255)
                cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCa                                                                                                     <  7:345307/135‑1 (MQ=255)
                    gCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAg                                                                                                   <  8:96120/133‑1 (MQ=255)
                       aGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                              <  8:34531/135‑1 (MQ=255)
                               aTTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGAtt                                                                                      <  7:397400/135‑1 (MQ=255)
                                      tgTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGTCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGACTAAAATTCAAGTAGGCGTTCGATTGTAcgca                                                                               <  3:96570/134‑1 (MQ=255)
                                         ccTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGGCTTAAATTAAAGTAGGCGTTCGATTTTAagcctc                                                                             >  8:181273/1‑134 (MQ=255)
                                              tAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGc                                                                       <  6:355372/135‑1 (MQ=255)
                                                 tAAAACGAGCATTACGCATTCACCTTGAGTTGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  4:423767/134‑1 (MQ=255)
                                                  aaaaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCgg                                                                    <  6:45384/134‑1 (MQ=255)
                                                    aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTg                                                                  <  7:55671/134‑1 (MQ=255)
                                                     aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                <  7:181273/135‑1 (MQ=255)
                                                      cGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAg                                                               <  8:129545/135‑1 (MQ=255)
                                                       gAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCGTTTTTTTTTCCTAAAATTCAAGTAGGCGTTCGATTGTAAGCATCGCATTCCGGTGAAGt                                                              >  2:160859/1‑135 (MQ=255)
                                                        aGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                             <  7:226458/135‑1 (MQ=255)
                                                        aGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGtt                                                             <  1:6051/135‑1 (MQ=255)
                                                          cATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCa                                                           <  8:69352/135‑1 (MQ=255)
                                                           aTTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCAt                                                          <  1:160761/135‑1 (MQ=255)
                                                            cagcGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATg                                                         <  3:133716/132‑1 (MQ=255)
                                                             tACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                        <  6:14400/135‑1 (MQ=255)
                                                                  aTTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAAGCATCGGATGCCGGTGAAGTTCAGGACTacc                                                   >  7:445722/1‑134 (MQ=255)
                                                                            aggggCTTAAGTGCTCCTATGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                          <  2:452429/131‑1 (MQ=255)
                                                                            aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                          <  6:59422/134‑1 (MQ=255)
                                                                            aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                           <  8:153835/133‑1 (MQ=255)
                                                                            aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                           <  5:77778/133‑1 (MQ=255)
                                                                            aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                           <  1:170110/133‑1 (MQ=255)
                                                                             gTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTATGCGTTAGATTGTACGCATAGCATGCCGGTGAAGTTCATGACTAAACCAATGAcct                                         >  4:133716/1‑133 (MQ=255)
                                                                                    aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCGAAAAATCAACGAGGCGTTCCAGTGTACCCAACCCCAGCCGGGGAAGGTCAAGACGACACCAATGCCCAGGCTCtt                                  >  7:101115/1‑134 (MQ=255)
                                                                                        gCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAAtg                             <  6:97655/135‑1 (MQ=255)
                                                                                        gCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAATTAGGCGGTCGATTTTACGCATCGCATGCCGGGGAAGTTAATTACTACACCAATGACCAGTCTCTTAAAtt                             >  6:466177/1‑134 (MQ=255)
                                                                                             tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGc                        <  6:172330/135‑1 (MQ=255)
                                                                                                ggTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCg                       <  6:51541/133‑1 (MQ=255)
                                                                                                ggTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCg                       <  5:378593/133‑1 (MQ=255)
                                                                                                   aaaaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTc                  <  7:16082/135‑1 (MQ=255)
                                                                                                    aaaGGTGGATATTTGGCGCCCATTCGGCGCCTTTTTTTTGGCAAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGGGAGCGCCCTCg                 >  3:177559/1‑135 (MQ=255)
                                                                                                     aaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGc                <  6:213764/135‑1 (MQ=255)
                                                                                                              atTTGGCGCCCATTGGGCGCCGTTTTTTTTCCTAAAAATCAAGTAGGCGGTCGATTGTAAGCATAGCATGCCGGTGAAGTTCATGACTACACCAAAGACCACTCTCTTCAATGTGAGGCGCCTCGCCGTCaaca        >  2:150642/1‑131 (MQ=255)
                                                                                                              atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATc         >  3:50038/1‑133 (MQ=255)
                                                                                                              atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTg       <  5:458522/135‑1 (MQ=255)
                                                                                                              atTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTg       <  1:142079/135‑1 (MQ=255)
                                                                                                                    cgcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg  <  1:150642/134‑1 (MQ=255)
                                                                                                                                   |                                                                                                                       
TTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB122/44162‑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: