Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	402,666	G→C	9.6%	P33A (CCT→GCT)	rlpA ←	minor lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	402,666	0	G	C	9.6%	208.8 / 10.8	73	P33A (CCT→GCT)	rlpA	minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (42/24);  new base C (3/4);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

TCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGA  >  minE/402597‑402734                                                                      |                                                                     tCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                     >  1:834845/1‑71 (MQ=255) tCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                     >  1:2923025/1‑71 (MQ=255)   gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                                    <  1:1165249/70‑1 (MQ=255)   gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTg                                                                   <  1:3126908/71‑1 (MQ=255)   gcgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTg                                                                   <  1:3071031/71‑1 (MQ=255)    cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCt                                                                    <  1:1907768/69‑1 (MQ=255)    cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGc                                                                  >  1:2359280/1‑71 (MQ=255)    cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGc                                                                  >  1:2363582/1‑71 (MQ=255)    cgTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGc                                                                  >  1:2702733/1‑71 (MQ=255)         aGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTac                                                             <  1:1242762/71‑1 (MQ=255)         aGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTa                                                              >  1:3343186/1‑70 (MQ=255)               tCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                       >  1:2320025/1‑71 (MQ=255)               tCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                       >  1:2294638/1‑71 (MQ=255)               tCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                       >  1:2130812/1‑71 (MQ=255)               tCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTAc                                                       >  1:1423523/1‑71 (MQ=255)                 gAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTacac                                                           <  1:2415106/65‑1 (MQ=255)                    aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                     >  1:2894435/1‑52 (MQ=255)                    aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                  >  1:2778022/1‑71 (MQ=255)                    aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                  >  1:3363615/1‑71 (MQ=255)                    aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                  <  1:1218127/71‑1 (MQ=255)                    aCGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                  >  1:1140073/1‑71 (MQ=255)                     cgcgGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTg                                                 <  1:1630515/71‑1 (MQ=255)                        gggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTg                                              <  1:2809636/71‑1 (MQ=255)                         ggTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTg                                              <  1:946039/70‑1 (MQ=255)                          gTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCt                                            <  1:98027/71‑1 (MQ=255)                             cGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGa                                          <  1:1537104/70‑1 (MQ=255)                             cGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGAc                                         >  1:1457264/1‑71 (MQ=255)                                cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACt                                                          >  1:300379/1‑51 (MQ=255)                                 ccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTa                                                        >  1:1685663/1‑52 (MQ=255)                                   gCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatca                                    >  1:594261/1‑70 (MQ=255)                                    cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCatcatc                                  <  1:455983/71‑1 (MQ=255)                                    cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCatcatc                                  >  1:897052/1‑71 (MQ=255)                                    cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCatcatc                                  <  1:3034942/71‑1 (MQ=255)                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCt                                                  >  1:703091/1‑53 (MQ=255)                                         ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTg                                                 >  1:85665/1‑51 (MQ=255)                                         ttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTg                                                 >  1:784820/1‑51 (MQ=255)                                          tCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg                             >  1:2257827/1‑70 (MQ=255)                                           cAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGa                                          >  1:243760/1‑56 (MQ=255)                                             aCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACtt                                                         <  1:2973039/39‑1 (MQ=255)                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCatcatc                                  <  1:447642/61‑1 (MQ=255)                                              cTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                        >  1:3526497/1‑71 (MQ=255)                                                ataGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGc                      >  1:32808/1‑71 (MQ=255)                                                 taGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGc                      <  1:1803103/70‑1 (MQ=255)                                                 taGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGc                      <  1:1408585/70‑1 (MQ=255)                                                  aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCtt                              >  1:496268/1‑61 (MQ=255)                                                  aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcg                    >  1:489078/1‑71 (MQ=255)                                                  aGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcg                    >  1:2631968/1‑71 (MQ=255)                                                   gggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACa                         >  1:1698934/1‑65 (MQ=255)                                                      ccGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAg                >  1:3164089/1‑71 (MQ=255)                                                       cGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAg                >  1:1675023/1‑70 (MQ=255)                                                        gTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACAATCATCGCTTGTACAt                        >  1:2371158/1‑61 (MQ=255)                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCga                                >  1:1146363/1‑49 (MQ=255)                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGc                                >  1:2725271/1‑50 (MQ=255)                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCTGCGAGCATTc           >  1:2888650/1‑71 (MQ=255)                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTc           >  1:2976128/1‑71 (MQ=255)                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTc           >  1:2156797/1‑71 (MQ=255)                                                           aCATACCGAAGGCTGCGGTACACTTACCGTCTGTTGCTg                                           >  1:3372376/1‑39 (MQ=255)                                                             aTACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCt         >  1:1959579/1‑71 (MQ=255)                                                               aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTa                           >  1:2290176/1‑51 (MQ=255)                                                               aCCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTa                           >  1:3666491/1‑51 (MQ=255)                                                                ccGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:2921952/71‑1 (MQ=255)                                                                ccGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:2386934/71‑1 (MQ=255)                                                                ccGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:768977/71‑1 (MQ=255)                                                                ccGCAGCCTGCGGTTTACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:423898/71‑1 (MQ=255)                                                                ccGCAGCCTGCGGTATACTTATCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:2579643/71‑1 (MQ=255)                                                                ccGCAGCCTGCGGTATACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:1066335/71‑1 (MQ=255)                                                                ccGCAGCCTGCGGTATACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc      <  1:1761182/71‑1 (MQ=255)                                                                  gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    <  1:239527/71‑1 (MQ=255)                                                                   cAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgcg   <  1:2385794/71‑1 (MQ=255)                                                                   cAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    >  1:3417188/1‑70 (MQ=255)                                                                   cAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    >  1:3171817/1‑70 (MQ=255)                                                                   cAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    >  1:3083049/1‑70 (MQ=255)                                                                    aGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGa  <  1:2146453/71‑1 (MQ=255)                                                                      gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    <  1:531657/67‑1 (MQ=255)                                                                      gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    <  1:269150/67‑1 (MQ=255)                                                                      gCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc    <  1:2427157/67‑1 (MQ=255)                                                                      |                                                                     TCGCGTTCAGTGGTTCGAAACGCGGGTCCGCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGA  >  minE/402597‑402734

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: ‑