Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	785,990	C→T	48.6%	A314A (GCG→GCA)	nhaB ←	sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	785,990	0	C	T	48.6%	1.5 / 65.3	74	A314A (GCG→GCA)	nhaB	sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (22/14);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCT  >  minE/785942‑786051                                                 |                                                              cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1019104/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:968877/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:878528/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:726204/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:714785/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:703998/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:518958/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:377684/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:3409151/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:3084568/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:2992367/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:2757496/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:2325073/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1955687/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1553080/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1498562/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1430455/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1387481/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                                           >  1:1167457/1‑69 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                                             >  1:929546/1‑67 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                                             >  1:2070431/1‑67 (MQ=255) cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGACAGCGCAGTCAccagcca                                           >  1:1340766/1‑69 (MQ=255)  aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                                           <  1:2974062/68‑1 (MQ=255)  aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                                           <  1:2766647/68‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:1866549/69‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:2022944/69‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:2034382/69‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:657516/69‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:409368/69‑1 (MQ=255)       aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:1386504/69‑1 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:2885790/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:463654/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     <  1:1390341/68‑1 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:3189941/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:793690/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:1087916/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:1634500/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:1600856/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGAAACGCCAGCGCAGTCACCAGCCAGACGcc                                     >  1:1447334/1‑68 (MQ=255)        atGACTGACAAACCAATCAAGCAAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCAAGACGcc                                     >  1:457354/1‑68 (MQ=255)           actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                                 <  1:565500/69‑1 (MQ=255)           actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                                 <  1:1460705/69‑1 (MQ=255)           actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                                 <  1:2809468/69‑1 (MQ=255)           actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                                 <  1:2940257/69‑1 (MQ=255)           actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                                 <  1:891224/69‑1 (MQ=255)                cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc                            >  1:2407936/1‑69 (MQ=255)                cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc                            >  1:2007468/1‑69 (MQ=255)                cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc                            >  1:1351573/1‑69 (MQ=255)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:1541287/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:556023/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:566585/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:539750/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:3133974/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:1755950/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:1929087/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:2039667/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:2061797/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:21044/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:251461/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:2794961/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:2989402/54‑1 (MQ=25)                   aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                                        <  1:3024253/54‑1 (MQ=25)                          aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACa                     >  1:2281766/1‑66 (MQ=255)                          aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACa                     >  1:3028986/1‑66 (MQ=255)                          aaGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGAACa                     >  1:2284462/1‑66 (MQ=255)                              cAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGa               >  1:581413/1‑68 (MQ=255)                                          gATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCt  <  1:1136200/69‑1 (MQ=255)                                          gATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCt  <  1:1103646/69‑1 (MQ=255)                                          gATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCt  <  1:842842/69‑1 (MQ=255)                                          gATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATACt  <  1:1466114/69‑1 (MQ=255)                                             gCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTa      <  1:1005222/62‑1 (MQ=255)                                             gCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTa      <  1:1679406/62‑1 (MQ=255)                                             gCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTa      <  1:1435351/62‑1 (MQ=255)                                             gCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGAATTTTa      <  1:3643/62‑1 (MQ=255)                                                 |                                                              CAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGACGGATTTTATCCT  >  minE/785942‑786051

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