Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 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
*ACB14344,2930T.93.3% 181.9 / 3.4 44intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (12/29);  total (14/30)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

ACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB143/44177‑44412
                                                                                                                    |                                                                                                                       
aCTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCa                                                                                                     <  3:232679/136‑1 (MQ=255)
 cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAAtta                                                                                                      >  6:36056/1‑133 (MQ=255)
    ggCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAg                                                                                                   <  3:273758/134‑1 (MQ=255)
     gCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAg                                                                                                <  1:158851/136‑1 (MQ=255)
       cAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                              <  5:166467/136‑1 (MQ=255)
       cAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                              <  2:290274/136‑1 (MQ=255)
               ccttGGCGGGTACTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGAtt                                                                                      <  4:178485/134‑1 (MQ=255)
               cATTGGCGGGTCCTGTTAGTAAAACGAGCATTACTCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGAtt                                                                                      <  6:288897/136‑1 (MQ=255)
                     cGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCG‑CTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcg                                                                                 >  1:186734/1‑135 (MQ=255)
                         tCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                            <  4:294369/136‑1 (MQ=255)
                         tCCTGTTAGTAAAACGAGCATTACGCAATCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                            <  4:112645/136‑1 (MQ=255)
                            tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgcat                                                                         >  1:196978/1‑136 (MQ=255)
                                 gTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTTCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  8:182920/135‑1 (MQ=255)
                                 gTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  5:45114/135‑1 (MQ=255)
                                 gTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  5:199707/135‑1 (MQ=255)
                                     aaCGATCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                <  4:187885/136‑1 (MQ=255)
                                     aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                <  3:214326/136‑1 (MQ=255)
                                     aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                <  3:292367/136‑1 (MQ=255)
                                             ttACGCATTCACCTTGCGTGGCTTAAGTTCTCATTTGGTAAAAGGTGGAGATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                        <  7:55872/136‑1 (MQ=255)
                                             ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                        <  3:160382/136‑1 (MQ=255)
                                             ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                        <  7:114536/136‑1 (MQ=255)
                                             ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                        <  7:208167/136‑1 (MQ=255)
                                                 gCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTac                                                    <  2:114848/136‑1 (MQ=255)
                                                       aCCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCaa                                               <  3:239830/135‑1 (MQ=255)
                                                        ccTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAAt                                              >  7:103690/1‑135 (MQ=255)
                                                          ttGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     <  7:223844/110‑1 (MQ=255)
                                                          ttGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                     >  8:223844/1‑110 (MQ=255)
                                                            gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                          <  5:103827/135‑1 (MQ=255)
                                                            gAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                           <  3:180679/134‑1 (MQ=255)
                                                            gAGTGGCGTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                          <  7:120243/135‑1 (MQ=255)
                                                                ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATTACCAGtct                                     >  5:212795/1‑136 (MQ=255)
                                                                 gCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtct                                     <  8:186877/135‑1 (MQ=255)
                                                                    tAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCtt                                  >  8:215980/1‑135 (MQ=255)
                                                                              tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCAAGCCGGTGAAGTTCATGATTAAACCAATGACCAGTCTCTTCATTGTTAGCg                       >  2:281704/1‑136 (MQ=255)
                                                                               ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTAAACCAATGACCAGTCTCTTCAATGTGAGCGa                      >  4:193437/1‑136 (MQ=255)
                                                                                   tAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTc                  <  2:92895/136‑1 (MQ=255)
                                                                                     aaaGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGc                >  8:225801/1‑136 (MQ=255)
                                                                                          tGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCa                                                                                                     <  7:57614/46‑1 (MQ=255)
                                                                                          tGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCa                                                                                                     >  8:57614/1‑46 (MQ=255)
                                                                                          tGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAg                                                                                                <  2:13362/51‑1 (MQ=255)
                                                                                          tGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAg                                                                                                >  1:13362/1‑51 (MQ=255)
                                                                                               atTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTGATGACTAAACCAATGAGCAGTCTCTTCAATGTGGGCGCCCACGGCGTCATCt        >  4:229222/1‑134 (MQ=255)
                                                                                               atTTGGCGCCCATTGGGCG‑CGTTTTTTTTGGCTAAAATTCAAGTAGGCGTTCGATTGTACGCATAGCATGCCGGTGAAGTTCATGACTAACCCAATGACCAGTCTCTTCAATGTGAGAGACCTCGCCGTCATcg        >  7:245523/1‑133 (MQ=255)
                                                                                                    gcgcCCATTGGGCG‑CTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg  <  2:186734/135‑1 (MQ=255)
                                                                                                                    |                                                                                                                       
ACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGG  >  ACB143/44177‑44412

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