Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	ACB143	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
*	ACB143	44,293	0	T	.	92.3%	151.2 / 9.0	40	intergenic (+124/‑367)	oprD_1/czcR_2	Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (9/27);  total (13/27)
Fisher's exact test for biased strand distribution p-value = 2.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTG  >  ACB143/44159‑44407                                                                                                                                       |                                                                                                                   cGATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCttt                                                                                                                  <  7:229287/137‑1 (MQ=255)          aaTTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCGttttttttgcct                                                                                                         >  4:236616/1‑133 (MQ=255)                 tAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGTTAAAAGGTGGATATTTGGCGCCCATTGGGCG‑CAATTTTTTTGCCAAAAAtta                                                                                                 >  2:364604/1‑136 (MQ=255)                   aCTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCaa                                                                                               >  2:147695/1‑137 (MQ=255)                       ggCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAg                                                                                           <  6:84732/137‑1 (MQ=255)                         ccAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                         <  5:187009/137‑1 (MQ=255)                           aGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑GTTTTTTTGCCTAAAATTCAAGTAGGCGt                                                                                       >  2:2793/1‑137 (MQ=255)                           aGTTATCCATTGGAGGTTAATTTTAAAAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTT‑CCTAAAATTCAAGTAGGCGt                                                                                       >  5:221430/1‑137 (MQ=255)                                   aTTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGt                                                                               <  1:246139/137‑1 (MQ=255)                                        cGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgca                                                                          >  4:315186/1‑137 (MQ=255)                                           gTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcg                                                                       <  4:308248/137‑1 (MQ=255)                                               tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATg                                                                   >  5:203257/1‑137 (MQ=255)                                                 ttAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑‑TTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                <  3:27616/137‑1 (MQ=255)                                                  tAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                <  1:157019/137‑1 (MQ=255)                                                  tAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                <  8:243799/137‑1 (MQ=255)                                                  tAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                <  5:260999/137‑1 (MQ=255)                                                   aGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCgg                                                               <  7:19789/137‑1 (MQ=255)                                                    gTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGt                                                              >  5:285850/1‑137 (MQ=255)                                                        aaCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAg                                                          <  7:256374/137‑1 (MQ=255)                                                                ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                         >  8:286251/1‑98 (MQ=255)                                                                ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                         <  7:286251/98‑1 (MQ=255)                                                                   cGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTac                                               <  7:333184/137‑1 (MQ=255)                                                                   cGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTac                                               <  1:227509/137‑1 (MQ=255)                                                                          aCCTTGAGTGTCTTAAGTGCTCCTTTGGTAAAATGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATg                                        <  6:297637/137‑1 (MQ=255)                                                                          aCCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATg                                        >  6:382877/1‑137 (MQ=255)                                                                           ccTTGAGTGGCTTAAGTGCTACTTTGGTAAAAGCTGGATATTTGGCGCCCAATGGCCGCCTATTTTTTTG‑CTAAAATTCAAGTATGCGTTCGAATGTACGCATCGCATGCCGGTGAAGATCATGACTACACCAAAGa                                       >  2:302885/1‑137 (MQ=255)                                                                            cTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                      <  6:380453/137‑1 (MQ=255)                                                                             ttGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                     <  6:285850/137‑1 (MQ=255)                                                                             ttGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                     <  6:221430/137‑1 (MQ=255)                                                                             ttGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAcc                                     <  5:23208/137‑1 (MQ=255)                                                                                  tGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtct                                <  5:382877/137‑1 (MQ=255)                                                                                         aGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTT‑CCTAAAATTGAAGTAGGCGTTCGATTGTACGCATCGCATTCCGGTGAAGTTCATGACTACACCAATGACCACTCTCTTCAAt                         >  3:142949/1‑137 (MQ=255)                                                                                                cTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGTCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCg                  <  1:239978/137‑1 (MQ=255)                                                                                                cTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCg                  <  4:292002/137‑1 (MQ=255)                                                                                                 tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGa                 >  4:144969/1‑137 (MQ=255)                                                                                                 tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGa                 >  2:238639/1‑137 (MQ=255)                                                                                                     gTAAAAGGTGGATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTc             <  7:330846/137‑1 (MQ=255)                                                                                                               tatatTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt   <  1:147695/136‑1 (MQ=255)                                                                                                               gATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt   <  2:139040/137‑1 (MQ=255)                                                                                                               gATATTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCt   <  2:15142/137‑1 (MQ=255)                                                                                                                atatTTGGCGCCCATTGGGCGCC‑TTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTg  <  3:315186/137‑1 (MQ=255)                                                                                                                atatTTGGCGCCCATTGGGCGCC‑GTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTg  <  2:283327/137‑1 (MQ=255)                                                                                                                                       |                                                                                                                   CGATTCTGTAATTCGCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTG  >  ACB143/44159‑44407

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