Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,208	A→G	18.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	18.1%	214.0 / 41.4	94	Y33H (TAT→CAT)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (44/33);  new base G (8/9);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AATACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTT  >  minE/980140‑980277                                                                     |                                                                      aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATa                                                                       <  1:3573158/69‑1 (MQ=255) aaTACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGa                                                                     >  1:2284359/1‑71 (MQ=255)   tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACa                                                                   <  1:1384619/71‑1 (MQ=255)   tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAc                                                                    >  1:2593799/1‑70 (MQ=255)   tACCAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACa                                                                   <  1:3495720/71‑1 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGa                                                                     >  1:1233850/1‑67 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                 >  1:1160507/1‑71 (MQ=255)     ccAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                 >  1:1206121/1‑71 (MQ=255)      cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACa                                                                   >  1:4088097/1‑68 (MQ=255)      cAAATTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACaaa                                                                 >  1:721144/1‑70 (MQ=255)        aaTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:3851220/1‑70 (MQ=255)         aTTCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAAcaca                                                             >  1:3740763/1‑71 (MQ=255)           tCAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                           >  1:3862870/1‑71 (MQ=255)            cAACTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATc                                                           <  1:2528144/70‑1 (MQ=255)               cTTCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                       >  1:4070436/1‑71 (MQ=255)                ttCCTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAAcaca                                                             <  1:3094662/64‑1 (MQ=255)                  ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctg                                                    >  1:1804681/1‑71 (MQ=255)                  ccTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctc                                                    >  1:2579687/1‑70 (MQ=255)                   cTGATTCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:2318074/1‑59 (MQ=255)                       ttCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATc                                               >  1:875312/1‑71 (MQ=255)                       ttCTCAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATc                                               >  1:121/1‑71 (MQ=255)                         ctcAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                              <  1:2435060/70‑1 (MQ=255)                         ctcAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCt                                              <  1:3402535/70‑1 (MQ=255)                          tcAGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCttt                                                        <  1:2182833/59‑1 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgc                                                      >  1:3890769/1‑59 (MQ=255)                            aGTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTc                                           <  1:2729400/70‑1 (MQ=255)                             gTTAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGAcc                                                                   <  1:3767553/45‑2 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCgg                                         <  1:789007/70‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCgg                                         <  1:1988618/70‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACa                                                                   <  1:31514/44‑1 (MQ=255)                              ttAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        <  1:1807274/71‑1 (MQ=255)                               tAATGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgc                                                      >  1:2197148/1‑56 (MQ=255)                                aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAaca                                                               >  1:3863843/1‑46 (MQ=255)                                aaTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAac                                                                >  1:2013423/1‑45 (MQ=255)                                 aTGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATc                                               <  1:419249/61‑1 (MQ=255)                                  tGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtgctgc                                                   <  1:1784762/56‑1 (MQ=255)                                  tGTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                    >  1:2174084/1‑71 (MQ=255)                                   gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtg                                                       >  1:94600/1‑51 (MQ=255)                                   gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACa                                                                   <  1:3589931/39‑1 (MQ=255)                                   gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        <  1:631902/66‑1 (MQ=255)                                   gTTTTGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGc                                        <  1:2133910/66‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTtg                                                       >  1:2145898/1‑50 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCa                                                 <  1:286058/56‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                  <  1:1950379/71‑1 (MQ=255)                                    ttttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                   >  1:440008/1‑70 (MQ=255)                                     tttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatg                                  >  1:2649215/1‑70 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCt                                              <  1:1476242/57‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:3677156/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 >  1:3872196/1‑70 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:1131090/70‑1 (MQ=255)                                      ttGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgatga                                 <  1:1027584/70‑1 (MQ=255)                                       tGAAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaa                                <  1:1728456/70‑1 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                        >  1:2114059/1‑71 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                       >  1:3257370/1‑40 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTtg                                                       >  1:1544146/1‑40 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCtt                                             <  1:3888601/50‑1 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtgat                                   <  1:4011060/60‑1 (MQ=255)                                              gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                        >  1:472130/1‑71 (MQ=255)                                               gtGCTCTGAATTGCTGAGAAATAGACAAACACATCt                                                          >  1:3403036/2‑36 (MQ=255)                                                tGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                    >  1:1596230/1‑57 (MQ=255)                                                 gatctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                     >  1:951510/3‑71 (MQ=255)                                                 gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTc                                           >  1:7342/1‑49 (MQ=255)                                                  ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                    <  1:3571260/71‑1 (MQ=255)                                                  ctctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                    <  1:2223545/71‑1 (MQ=255)                                                  ctctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                    <  1:3768585/71‑1 (MQ=255)                                                   tctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                            >  1:647152/1‑46 (MQ=255)                                                    ctGAATTGCTGAGAAATAGACAAACACATCTTtgc                                                      >  1:2153759/1‑35 (MQ=255)                                                    ctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCg                                       >  1:3476832/1‑50 (MQ=255)                                                     tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCttt                                            >  1:1919532/1‑44 (MQ=255)                                                     tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                  >  1:831058/1‑70 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATc                                               <  1:2764071/41‑1 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCttt                                            >  1:4116653/1‑44 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtttt                   <  1:1909561/69‑1 (MQ=255)                                                     tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                    <  1:2728527/68‑1 (MQ=255)                                                      gAATTGCTGAGAAATAGACAAACACATCTTtgctgc                                                   >  1:2336661/1‑36 (MQ=255)                                                      gAATTGCTGAGAAATAGACAAACACATCTTtgctgc                                                   >  1:742342/1‑36 (MQ=255)                                                      gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                  >  1:2104409/1‑69 (MQ=255)                                                       aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg               >  1:3376469/1‑71 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGcc                                                  >  1:2659747/1‑36 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                 <  1:166648/69‑1 (MQ=255)                                                       aaTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                 <  1:2260311/69‑1 (MQ=255)                                                       aaTTGCTGAGAAATAGAAAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg               >  1:1303568/1‑71 (MQ=255)                                                        aTTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg              >  1:3111000/1‑71 (MQ=255)                                                           gCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGt                                      >  1:68485/1‑44 (MQ=255)                                                            cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGtt           <  1:3483792/70‑1 (MQ=255)                                                             tGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                          <  1:2606954/54‑1 (MQ=255)                                                             tGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                 <  1:1210051/63‑1 (MQ=255)                                                              gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGaa                                <  1:1860063/47‑1 (MQ=255)                                                              gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa        >  1:2814006/1‑71 (MQ=255)                                                               agaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGtga                                    >  1:3962688/1‑42 (MQ=255)                                                                 aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                <  1:3001028/60‑1 (MQ=255)                                                                 aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                <  1:3005820/60‑1 (MQ=255)                                                                  aaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                >  1:3921384/1‑59 (MQ=255)                                                                  aaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAt    >  1:1910277/1‑71 (MQ=255)                                                                    tAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt  <  1:656965/71‑1 (MQ=255)                                                                     aGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                    >  1:875802/1‑52 (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: ‑