Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,208	A→G	11.1%	Y33H (TAT→CAT)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	980,208	0	A	G	11.1%	252.4 / 19.8	90	Y33H (TAT→CAT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (35/45);  new base G (2/8);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATT  >  minE/980140‑980276                                                                     |                                                                     aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                      <  1:2447076/69‑1 (MQ=255) aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                      <  1:2280126/69‑1 (MQ=255) aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                      <  1:1852619/69‑1 (MQ=255) aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTATGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                      <  1:1227352/69‑1 (MQ=255)  aTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   <  1:736782/71‑1 (MQ=255)   tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   >  1:3430386/1‑70 (MQ=255)   tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   >  1:1232744/1‑70 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                >  1:2865965/1‑71 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaa                                                                 <  1:1752815/70‑1 (MQ=255)      cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAac                                                               <  1:979700/71‑1 (MQ=255)       aaaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                              <  1:3421969/71‑1 (MQ=255)       aaaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAaca                                                              <  1:947633/71‑1 (MQ=255)         aTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGa                                                                    >  1:399770/1‑63 (MQ=255)         aTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGa                                                                    >  1:2097277/1‑63 (MQ=255)             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCtt                                                        <  1:934290/71‑1 (MQ=255)               cTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   >  1:2990362/1‑58 (MQ=255)               cTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   >  1:2496475/1‑58 (MQ=255)               cTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                      >  1:2058155/1‑71 (MQ=255)                  ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                      <  1:2683498/68‑1 (MQ=255)                  ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                      <  1:3294111/68‑1 (MQ=255)                   cTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctg                                                   <  1:752571/70‑1 (MQ=255)                    tGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   <  1:1915313/53‑1 (MQ=255)                    tGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   <  1:2688280/53‑1 (MQ=255)                    tGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                   <  1:1645844/53‑1 (MQ=255)                    tGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaa                                                                 <  1:2117209/55‑1 (MQ=255)                    tGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaa                                                                 <  1:2989196/55‑1 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaa                                                                 <  1:214750/53‑1 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCAt                                               >  1:3424598/1‑71 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCAt                                               >  1:2937598/1‑71 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCg                                         <  1:1700214/71‑1 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCg                                         >  1:3693572/1‑71 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCg                                         >  1:2761843/1‑71 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCg                                         >  1:1162814/1‑71 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                           <  1:2069027/51‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                           <  1:2064953/51‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                           <  1:316227/51‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTTCTGCCATCTTTCGGc                                       <  1:2546637/71‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCgg                                        <  1:3475706/70‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                       <  1:3151706/71‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                       <  1:1619732/71‑1 (MQ=255)                                aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCg                                      <  1:1948884/70‑1 (MQ=255)                                 aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                     <  1:3129800/70‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                  >  1:3315966/1‑70 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                          <  1:578706/46‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                           <  1:2314251/61‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                           <  1:717348/61‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                  >  1:2316499/1‑70 (MQ=255)                                     tttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                 >  1:2412359/1‑70 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:1110907/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:1671398/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCtt                                                        >  1:2826060/1‑46 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:539108/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:1787966/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:1574133/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                   >  1:3299124/1‑67 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCCTTCGGCGtgatga                                <  1:2636730/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCCTTCGGCGtgatga                                <  1:2632424/70‑1 (MQ=255)                                        gAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctg                                                   >  1:2521112/1‑49 (MQ=255)                                        gAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctg                                                   >  1:3089467/1‑49 (MQ=255)                                        gAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgct                                                    >  1:558555/1‑48 (MQ=255)                                         aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                             >  1:2938397/1‑70 (MQ=255)                                         aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                             >  1:82546/1‑70 (MQ=255)                                         aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                             >  1:3358330/1‑70 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                            >  1:588817/1‑70 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                           <  1:3600398/71‑1 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                           <  1:39410/71‑1 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                           <  1:3296875/71‑1 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCa                           <  1:3354565/71‑1 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                             >  1:35912/1‑68 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAc                             >  1:2756414/1‑68 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                          <  1:3532434/71‑1 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGAGTGATGAAAc                             >  1:2949511/1‑68 (MQ=255)                                               ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                      <  1:2986436/71‑1 (MQ=255)                                               ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                      <  1:2912701/71‑1 (MQ=255)                                               ttGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAcc                      <  1:1195259/71‑1 (MQ=255)                                                 gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                    >  1:2256026/1‑71 (MQ=255)                                                   tctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                   <  1:1807807/70‑1 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:2438203/55‑1 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                <  1:324109/55‑1 (MQ=255)                                                      gAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATGTTTCGGCGTGATGAAACCaaa                         >  1:1341911/1‑61 (MQ=38)                                                      gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                   >  1:1440872/1‑67 (MQ=255)                                                      gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATATTTCGGCGTGATGAAACCAAAACCttt                   >  1:1281004/1‑67 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCGAAACCTTTTTc                <  1:2127203/69‑1 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                <  1:1026343/69‑1 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa               <  1:642689/70‑1 (MQ=255)                                                            cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa         <  1:547626/71‑1 (MQ=255)                                                              gagaAATAGACAAACACATCTTTGCTGCCATCTTTCg                                         <  1:100266/37‑1 (MQ=255)                                                              gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa       >  1:256963/1‑71 (MQ=255)                                                                gaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc      >  1:1501188/1‑70 (MQ=255)                                                                 aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg             >  1:1246287/1‑62 (MQ=255)                                                                   aTAGACAAACACATCTTTTCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt  <  1:2231721/71‑1 (MQ=255)                                                                   aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg             <  1:601740/60‑1 (MQ=255)                                                                   aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt  >  1:2115194/1‑71 (MQ=255)                                                                   aTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt  >  1:363315/1‑71 (MQ=255)                                                                     |                                                                     AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATT  >  minE/980140‑980276

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