Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	402,650	C→T	45.0%	G38D (GGC→GAC)	rlpA ←	minor lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	402,650	0	C	T	45.0%	16.6 / 46.3	40	G38D (GGC→GAC)	rlpA	minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (16/6);  new base T (5/13);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGC  >  minE/402581‑402717                                                                      |                                                                    tAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGcc                                                                    >  1:722658/1‑71 (MQ=255)   aTCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGTCGt                                                                  <  1:2586617/71‑1 (MQ=255)      cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGTCGTTa                                                                <  1:2625281/70‑1 (MQ=255)      cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGTCGTTa                                                                <  1:3544360/70‑1 (MQ=255)                  cgcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                   >  1:3029002/1‑71 (MQ=255)                  cgcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                   >  1:3232101/1‑71 (MQ=255)                   gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTg                                                  >  1:4007554/1‑71 (MQ=255)                    cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGTCGTTACATACCGCAGGCTGc                                                 <  1:4203825/71‑1 (MQ=255)                        cAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGc                                                                     >  1:1162057/1‑47 (MQ=255)                            ggTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                         >  1:1912257/1‑71 (MQ=255)                                 gAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTAcc                                     <  1:1034449/70‑1 (MQ=255)                                           tCCGCCCCGCTAATTTCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTg                          <  1:871381/71‑1 (MQ=255)                                                 ccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcat                     >  1:1561237/1‑70 (MQ=255)                                                  cGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatca                   <  1:3772803/71‑1 (MQ=255)                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGc                                                 >  1:2504215/1‑37 (MQ=255)                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc               >  1:1746810/1‑71 (MQ=255)                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc               >  1:3060584/1‑71 (MQ=255)                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc               >  1:3050164/1‑71 (MQ=255)                                                       aTTTCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTa                                             >  1:1942150/1‑40 (MQ=255)                                                       aTTTCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTa                                             >  1:2314519/1‑40 (MQ=255)                                                       aTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCt              >  1:4121672/1‑71 (MQ=255)                                                        tttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatca                   <  1:1104995/65‑1 (MQ=255)                                                        tttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTAACcatca                   <  1:1301271/65‑1 (MQ=255)                                                         ttCAACTATATGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg            >  1:2625656/1‑71 (MQ=255)                                                         ttCAACTATAGGGTCTTTACATACCGCAGGCTGCGGTACACTTACCGTc                                  <  1:1757768/49‑1 (MQ=39)                                                         ttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTc                                  <  1:2253177/49‑1 (MQ=255)                                                         ttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTc                                  <  1:3878587/49‑1 (MQ=255)                                                         ttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg            >  1:3339098/1‑71 (MQ=255)                                                         ttCAACTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg            >  1:1344348/1‑71 (MQ=255)                                                         ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGa                         >  1:1136820/1‑58 (MQ=255)                                                            aaCTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTa          <  1:3681672/70‑1 (MQ=255)                                                             aCTATAGGGTCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACa        <  1:3282075/71‑1 (MQ=255)                                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt       >  1:2840464/1‑71 (MQ=255)                                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt       >  1:2639155/1‑71 (MQ=255)                                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt       >  1:1745293/1‑71 (MQ=255)                                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt       >  1:1047966/1‑71 (MQ=255)                                                                  aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGc     >  1:4046990/1‑69 (MQ=255)                                                                   tggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc  <  1:3819734/70‑1 (MQ=255)                                                                   gggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc  <  1:2723260/71‑1 (MQ=255)                                                                   gggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc  <  1:3578787/71‑1 (MQ=255)                                                                   gggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc  <  1:1264180/71‑1 (MQ=255)                                                                    ggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgc  <  1:996320/70‑1 (MQ=255)                                                                      |                                                                    TAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGC  >  minE/402581‑402717

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