Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	785,972	A→G	44.0%	V320V (GTT→GTC)	nhaB ←	sodium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	785,972	0	A	G	44.0%	18.2 / 64.6	83	V320V (GTT→GTC)	nhaB	sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base A (23/23);  new base G (22/14);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGA  >  minE/785912‑786038                                                             |                                                                   cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTcc                                                              >  1:3037841/1‑67 (MQ=255) cTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTcc                                                              >  1:3189043/1‑67 (MQ=255)       ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                       >  1:1598776/1‑68 (MQ=255)       ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                       >  1:2193250/1‑68 (MQ=255)       ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                       >  1:2031303/1‑68 (MQ=255)       ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                       >  1:253054/1‑68 (MQ=255)       ggTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGAt                                                       >  1:920084/1‑68 (MQ=255)         tGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:2492870/68‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCTACTTCCGCCAGATGc                                                     <  1:1759469/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:221651/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:549189/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:2238912/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:956952/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:1236422/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:2618972/54‑1 (MQ=255)                       gatgTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGc                                                     <  1:1807485/54‑1 (MQ=255)                          gTTGCCAGAATAATGACTGACAAACCAATCAAGGCAACTTCCGCCAGATGc                                                     <  1:2793286/51‑1 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1019104/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:726204/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:2325073/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:2757496/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:2992367/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:968877/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:3084568/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:878528/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:3409151/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:377684/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:518958/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:703998/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:714785/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1553080/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1167457/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1498562/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1430455/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1955687/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagcca                              >  1:1387481/1‑69 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                                >  1:2070431/1‑67 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCAccagc                                >  1:929546/1‑67 (MQ=255)                               cAGAATAATGACCGACAAACCAATCAAACCGACTTCTGCCAGATGCAGTGACAGCGCAGTCAccagcca                              >  1:1340766/1‑69 (MQ=255)                                aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                              <  1:2766647/68‑1 (MQ=255)                                aGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCAccagcca                              <  1:2974062/68‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:2022944/69‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:1386504/69‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:657516/69‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:1866549/69‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:2034382/69‑1 (MQ=255)                                     aatGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:409368/69‑1 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:793690/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        <  1:1390341/68‑1 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:463654/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:3189941/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:1087916/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:1600856/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:1634500/1‑68 (MQ=255)                                      atGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGcc                        >  1:2885790/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:891224/69‑1 (MQ=255)                                         actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                    <  1:2940257/69‑1 (MQ=255)                                         actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                    <  1:1460705/69‑1 (MQ=255)                                         actgacAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCaata                    <  1:2809468/69‑1 (MQ=255)                                              cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc               >  1:1351573/1‑69 (MQ=255)                                              cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc               >  1:2407936/1‑69 (MQ=255)                                              cAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGc               >  1:2007468/1‑69 (MQ=255)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:1929087/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:556023/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:251461/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:566585/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:3024253/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:1755950/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:539750/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:2794961/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:1541287/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:3133974/54‑1 (MQ=25)                                                 aCCAATCAAACCGACTTCTGCCAGATGCAGTGCCAGCGCAGTCACCAGCCAGAc                           <  1:2989402/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)                                                             |                                                                   CTCATCGGTGACACCGGTCAATGATGTTGCCAGAATAATGACTGACAAACCAATCAAGCCAACTTCCGCCAGATGCAACGCCAGCGCAGTCACCAGCCAGACGCCAATAATCGCCTGGACAATCAGA  >  minE/785912‑786038

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