Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,208	A→G	39.7%	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	39.7%	43.2 / 72.5	68	Y33H (TAT→CAT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (13/14);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTT  >  minE/980140‑980277                                                                     |                                                                      aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                       <  1:2655852/69‑1 (MQ=255) aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAg                                                                      <  1:1881416/70‑1 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                 <  1:3522901/71‑1 (MQ=255)      cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACa                                                                   <  1:2167497/68‑1 (MQ=255)      cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAac                                                                <  1:1969351/71‑1 (MQ=255)        aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:366689/1‑70 (MQ=255)        aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:1466793/1‑70 (MQ=255)        aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:1819891/1‑70 (MQ=255)        aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAac                                                                >  1:2721994/1‑69 (MQ=255)         aTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAAcaca                                                             >  1:3767170/1‑71 (MQ=255)          ttCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACAt                                                            >  1:3633205/1‑71 (MQ=255)            cAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                           <  1:3826936/70‑1 (MQ=255)             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCtt                                                         >  1:262427/1‑71 (MQ=255)             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCt                                                          <  1:1022604/70‑1 (MQ=255)             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCtt                                                         >  1:3312039/1‑71 (MQ=255)             aaCTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCt                                                          <  1:2038605/70‑1 (MQ=255)                ttCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAac                                                                <  1:1321959/61‑1 (MQ=255)                ttCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                       >  1:392052/1‑70 (MQ=255)                 tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgct                                                     >  1:3336069/1‑71 (MQ=255)                 tCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgc                                                      <  1:3686542/70‑1 (MQ=255)                   cTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                           <  1:1613245/63‑1 (MQ=255)                   cTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                           <  1:369345/63‑1 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                            <  1:2814689/59‑1 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACAt                                                            <  1:397175/59‑1 (MQ=255)                      aTTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCAt                                                <  1:1300839/71‑1 (MQ=255)                          tcAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtgct                                                     >  1:3478178/1‑62 (MQ=255)                            aGTTAATGTTTTGAAATCCTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCtt                                             <  1:2524586/68‑1 (MQ=38)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGcc                                                  >  1:3467360/1‑61 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGc                                        >  1:2636296/1‑71 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        <  1:3820552/71‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        <  1:1732959/71‑1 (MQ=255)                                 aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                     <  1:1122383/71‑1 (MQ=255)                                 aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                     <  1:3160479/71‑1 (MQ=255)                                     tttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAAcaca                                                             >  1:3895556/1‑43 (MQ=255)                                     tttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgct                                                     >  1:82831/1‑51 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCt                                                          >  1:2817759/1‑45 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:2133592/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:266283/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:621809/70‑1 (MQ=255)                                       tGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctg                                                    <  1:2134788/50‑1 (MQ=255)                                         aaaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                             <  1:3570037/71‑1 (MQ=255)                                         aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                             <  1:3829845/71‑1 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        >  1:107367/1‑59 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        >  1:3214165/1‑59 (MQ=255)                                          aaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                             <  1:79450/70‑1 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                           <  1:1164842/71‑1 (MQ=255)                                           aTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                           <  1:74153/71‑1 (MQ=255)                                            tCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                           >  1:2807331/1‑70 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                          <  1:1274742/69‑1 (MQ=255)                                                 gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                             >  1:317576/1‑47 (MQ=255)                                                 gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                     >  1:3253011/1‑71 (MQ=255)                                                    ctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaa                                >  1:2839880/1‑57 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                    >  1:2051136/1‑52 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                             >  1:1881350/1‑59 (MQ=255)                                                       aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg               >  1:2917685/1‑71 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg               >  1:1779617/1‑71 (MQ=255)                                                        aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg              >  1:1237321/1‑71 (MQ=255)                                                        aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg              >  1:1634379/1‑71 (MQ=255)                                                         ttGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                      <  1:1176248/46‑1 (MQ=255)                                                            cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa          <  1:2521743/71‑1 (MQ=255)                                                            cTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTa          <  1:3571031/71‑1 (MQ=255)                                                              gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                           <  1:2008069/52‑1 (MQ=255)                                                               agaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaaa                               >  1:1054258/1‑47 (MQ=255)                                                                gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtg                                     >  1:3158922/1‑40 (MQ=255)                                                                gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaaa                         >  1:2162911/1‑52 (MQ=255)                                                                gaAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc       >  1:3664975/1‑70 (MQ=255)                                                                 aaaTGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCa     >  1:1007396/1‑71 (MQ=255)                                                                  aaTGGACAAACACATCTTTGCTGCCATCTTTCGGAGTGATGAAACCAAAACCTTTTTCAg               <  1:494509/60‑1 (MQ=25)                                                                    tAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt  <  1:2211814/71‑1 (MQ=255)                                                                     |                                                                      AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTT  >  minE/980140‑980277

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