Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,248	A→G	12.8%	G19G (GGT→GGC)	cspI ←	cold shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	980,248	0	A	G	12.8%	181.7 / 16.4	70	G19G (GGT→GGC)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (28/33);  new base G (3/6);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

AAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTT  >  minE/980180‑980317                                                                     |                                                                      aaaTCGTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                                     <  1:3570037/71‑1 (MQ=255) aaaTCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                                     <  1:3829845/71‑1 (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)         gCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                             >  1:3253011/1‑71 (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)                    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)                        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)                           atggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAtt                                           <  1:987032/68‑1 (MQ=255)                            tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                                     >  1:1273661/3‑60 (MQ=255)                            tggACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                                     >  1:232045/3‑60 (MQ=255)                            tAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                          <  1:2211814/71‑1 (MQ=255)                             gcaCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAttt                                          >  1:2548954/3‑70 (MQ=255)                              gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                                   >  1:1701962/1‑44 (MQ=255)                              gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                                   >  1:2141863/1‑44 (MQ=255)                              caCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttta                                        <  1:658298/70‑3 (MQ=255)                                cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAAcc                                                                     <  1:2518770/40‑1 (MQ=255)                                cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                       >  1:1534400/1‑70 (MQ=255)                                cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                                      <  1:2728011/71‑1 (MQ=255)                                cAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACt                                      <  1:1450160/71‑1 (MQ=255)                                 aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACTa                                     <  1:3525091/71‑1 (MQ=255)                                  aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                                    <  1:3226106/71‑1 (MQ=255)                                  aaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaa                                    <  1:1430417/71‑1 (MQ=255)                                   acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                                     >  1:2235155/1‑53 (MQ=255)                                   acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                                   >  1:3760251/1‑71 (MQ=255)                                     acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAcc                                 <  1:2092435/68‑1 (MQ=9)                                     acgtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                                  >  1:3505348/4‑70 (MQ=9)                                     acaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAcc                                 <  1:3379269/71‑1 (MQ=255)                                             gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAttt                         >  1:1066236/1‑71 (MQ=255)                                             gctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                          <  1:2909648/70‑1 (MQ=255)                                              gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt                        <  1:2536224/69‑1 (MQ=255)                                              gagcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt                        <  1:3804793/69‑1 (MQ=255)                                               tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAttt                                          <  1:2182225/52‑1 (MQ=255)                                               tgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCAt                                            <  1:531124/50‑1 (MQ=255)                                               agcCATCGTCAGGAGTGATGAAGCCGAAGCCTTTGTCAGCGTTGAACCATTTTACGATACCAGTCAtttt                        <  1:2864534/69‑1 (MQ=255)                                                 ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                       <  1:833754/53‑1 (MQ=255)                                                 ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      <  1:2761567/70‑1 (MQ=255)                                                  cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                        <  1:3158695/35‑1 (MQ=255)                                                  cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaa                                                 >  1:2866869/1‑42 (MQ=255)                                                  cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                                        >  1:400185/1‑51 (MQ=255)                                                  cATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                    <  1:1984849/71‑1 (MQ=255)                                                   aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                       <  1:2663073/35‑1 (MQ=255)                                                   aTCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                  <  1:2536476/56‑1 (MQ=255)                                                    tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                                      <  1:1487346/35‑1 (MQ=255)                                                    tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAgg                                                      <  1:1740353/35‑1 (MQ=255)                                                     cTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTg                       >  1:1881533/1‑65 (MQ=255)                                                      tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                          >  1:3310296/1‑61 (MQ=255)                                                      tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACa                >  1:1521667/1‑71 (MQ=255)                                                         ctgCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATa              <  1:2200645/68‑1 (MQ=255)                                                         cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa             >  1:1203435/1‑71 (MQ=255)                                                         cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATaa             >  1:3118615/1‑71 (MQ=255)                                                           gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAtt                          <  1:3418694/56‑1 (MQ=255)                                                           gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt           <  1:3671569/71‑1 (MQ=255)                                                            cGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat            <  1:3395416/69‑1 (MQ=255)                                                            agTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAAc                                  >  1:825724/2‑47 (MQ=16)                                                                  gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                      >  1:3460537/1‑53 (MQ=255)                                                                  gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCt    >  1:3817724/1‑71 (MQ=255)                                                                  gaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCt    >  1:2399905/1‑71 (MQ=255)                                                                    aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt  >  1:3280759/1‑71 (MQ=255)                                                                     |                                                                      AAATCGTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTT  >  minE/980180‑980317

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