Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,355,219	A→C	6.1%	M327R (ATG→AGG)	yihE ←	predicted kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,355,219	0	A	C	6.1%	218.5 / 9.7	82	M327R (ATG→AGG)	yihE	predicted kinase
Reads supporting (aligned to +/- strand):  ref base A (38/39);  new base C (2/3);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AACGCTAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGTCG  >  minE/2355152‑2355286                                                                    |                                                                    aacgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACAt                                                                    <  1:2159972/69‑1 (MQ=255)  acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAgg                                                                 >  1:3333082/1‑71 (MQ=255)  acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                                 <  1:780513/71‑3 (MQ=255)  acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                                 <  1:3848179/71‑3 (MQ=255)  acgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCAATTAACACCTatg                                                                 <  1:1484006/71‑3 (MQ=255)   cgctaaaACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGgt                                                                <  1:158712/71‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:3674646/69‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:275768/69‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:3124976/69‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                           <  1:1216852/71‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                           <  1:3833444/71‑1 (MQ=255)        aaaCTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTaa                                                           <  1:543450/71‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:3166888/66‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:2092275/66‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:2087683/66‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:2381166/66‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:2773391/66‑1 (MQ=255)           cTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGtt                                                             <  1:3553002/66‑1 (MQ=255)              aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      >  1:1531569/1‑70 (MQ=255)              aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      <  1:1972541/70‑1 (MQ=255)              aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      >  1:3778880/1‑70 (MQ=255)              aaCCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      >  1:767628/1‑70 (MQ=255)                      agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGtt                                             >  1:2889840/1‑71 (MQ=255)                      agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                              >  1:2832942/1‑70 (MQ=255)                      agcgccagcCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGt                                              >  1:1308963/1‑70 (MQ=255)                             gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                     >  1:1592025/1‑40 (MQ=255)                             gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                     >  1:727820/1‑40 (MQ=255)                             gccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACa                                                                     >  1:3384802/1‑40 (MQ=255)                              ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATATGTGTTAATTGCa                                                      <  1:2432538/54‑1 (MQ=255)                              ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      <  1:2424189/54‑1 (MQ=255)                              ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      <  1:2510367/54‑1 (MQ=255)                              ccaAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGa                                     <  1:1509867/71‑1 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:2560214/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:105668/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:3400854/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:2595720/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:2555435/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:2296033/1‑71 (MQ=255)                               caAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGaa                                    >  1:2253/1‑71 (MQ=255)                                 aaTCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCa                                                      <  1:783016/51‑1 (MQ=255)                                   tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                                >  1:1327012/1‑71 (MQ=255)                                   tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                                >  1:473590/1‑71 (MQ=255)                                   tCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtt                                 >  1:2543942/1‑70 (MQ=255)                                    cTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                                  >  1:107521/1‑36 (MQ=255)                                    cTTTTTCATGATCTACTCTCTCCGATTAATACATAg                                                                  >  1:1359769/1‑36 (MQ=255)                                    cTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGa                                     >  1:3080129/1‑65 (MQ=255)                                    cTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                                >  1:3808099/1‑70 (MQ=255)                                    cTTTTTCATGATCTACTCTCTCCGATTAACAGCTATGTATTAAtcg                                                        >  1:109306/1‑44 (MQ=16)                                    cTTTTTCATGATCTACTCTCTCCGATTAACACCTATGTATTAAtcg                                                        >  1:153995/1‑44 (MQ=18)                                       tttCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGcc                            >  1:509583/1‑71 (MQ=255)                                           aTGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGt                         <  1:4061963/70‑1 (MQ=255)                                              aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                         <  1:1086580/51‑1 (MQ=255)                                              aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                         <  1:1648194/51‑1 (MQ=255)                                              aTCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTg                                         <  1:3447282/51‑1 (MQ=255)                                                   ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                      <  1:3367610/49‑1 (MQ=255)                                                   ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                      <  1:3719996/49‑1 (MQ=255)                                                   ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAg                                      <  1:3074027/49‑1 (MQ=255)                                                   ctctctCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATaaaa                <  1:1339758/71‑1 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATaaaa                >  1:2529207/1‑69 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt              >  1:331933/1‑71 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt              >  1:3856149/1‑71 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt              >  1:3035723/1‑71 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGt              >  1:379008/1‑71 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCCATAAAAGt              >  1:3105047/1‑71 (MQ=255)                                                     ctctcCGATTAATACATAGGTGTTAATAGCAAAGGGGGTTCTTGTAGAACTTTTGcc                            >  1:1722492/1‑57 (MQ=255)                                                      tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg               <  1:2705431/69‑1 (MQ=255)                                                      tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg               <  1:65672/69‑1 (MQ=255)                                                      tctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAg               <  1:3209840/69‑1 (MQ=255)                                                       ctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                               <  1:2147930/52‑1 (MQ=255)                                                       ctcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                               <  1:1800615/52‑1 (MQ=255)                                                        tcCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACtttt                               <  1:3861975/51‑1 (MQ=255)                                                          cGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg         >  1:601610/1‑71 (MQ=255)                                                            aTTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACttt                                >  1:270491/1‑46 (MQ=255)                                                              tAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTg     >  1:4099036/1‑71 (MQ=255)                                                              tAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTg     >  1:889286/1‑71 (MQ=255)                                                                aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                  <  1:3935261/40‑1 (MQ=255)                                                                aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                  <  1:3280183/40‑1 (MQ=255)                                                                aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACt                                  <  1:3740095/40‑1 (MQ=255)                                                                aTACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGt    <  1:3990424/70‑1 (MQ=255)                                                                aTACATAGGTGTTAATTGCAAAGGGGGTTCTTCTAGAACt                                  <  1:1645190/40‑1 (MQ=255)                                                                  aCATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGTCg  >  1:988840/1‑70 (MQ=255)                                                                   cATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg         >  1:3858115/1‑62 (MQ=255)                                                                   cATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg         >  1:3227191/1‑62 (MQ=255)                                                                   cATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCg         >  1:2897831/1‑62 (MQ=255)                                                                   cATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCATGTTCTATAAAAGTCGCCg         >  1:1545882/1‑62 (MQ=255)                                                                    |                                                                    AACGCTAAAACTAAACCAGCCAGCGCCAGCCAAATCTTTTTCATGATCTACTCTCTCCGATTAATACATAGGTGTTAATTGCAAAGGGGGTTCTTGTAGAACTTTTGCCTGTTCTATAAAAGTCGCCGTCTGTCG  >  minE/2355152‑2355286

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