Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	980,248	A→G	6.9%	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	6.9%	199.8 / 6.1	72	G19G (GGT→GGC)	cspI	cold shock protein
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base G (3/2);  total (38/34)
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.03e-01

GTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTT  >  minE/980185‑980318                                                                |                                                                       gTTGCTCTGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                 >  1:2114059/1‑71 (MQ=255) gTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAAc                                                                 >  1:472130/1‑71 (MQ=255)    gatctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                              >  1:951510/3‑71 (MQ=255)     ctctGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                             <  1:3571260/71‑1 (MQ=255)     ctctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                             <  1:2223545/71‑1 (MQ=255)     ctctGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                             <  1:3768585/71‑1 (MQ=255)        tGAATTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                           >  1:831058/1‑70 (MQ=255)        tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtttt                                                            <  1:1909561/69‑1 (MQ=255)        tGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                             <  1:2728527/68‑1 (MQ=255)         gAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttttt                                                           >  1:2104409/1‑69 (MQ=255)          aaTTGCTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAg                                                        >  1:3376469/1‑71 (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)               cTGAGAAATGGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGtt                                                    <  1:3483792/70‑1 (MQ=255)                tGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                                   <  1:2606954/54‑1 (MQ=255)                tGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTc                                                          <  1:1210051/63‑1 (MQ=255)                 gagaAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTaaa                                                 >  1:2814006/1‑71 (MQ=255)                    aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         <  1:3005820/60‑1 (MQ=255)                    aaaTAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         <  1:3001028/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)                         gACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                                         >  1:3146072/1‑71 (MQ=255)                          aCAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGt                                                     >  1:232536/1‑58 (MQ=255)                          aCAAACACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCAtttt                                          <  1:2770728/69‑2 (MQ=255)                            aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         >  1:681447/1‑52 (MQ=255)                            aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         >  1:1616963/1‑52 (MQ=255)                            aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         >  1:2651381/1‑52 (MQ=255)                            aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCa                                                         >  1:1036952/1‑52 (MQ=255)                            aaaCACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAAc                                                <  1:2570456/61‑1 (MQ=255)                            aaaCACATCTTTACTACCATCAACAGGAGAAATAAAGCCGAAACCTTTATCAGCGTTAAACCATTTTACt                                       >  1:304054/1‑70 (MQ=255)                              acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCaa                                                                    <  1:892674/39‑1 (MQ=255)                              acacATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                                    <  1:3288749/71‑1 (MQ=255)                               cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCaaa                                                                   <  1:832539/39‑1 (MQ=255)                               cacaTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                   >  1:3043582/1‑71 (MQ=255)                                 caTCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                   >  1:3386561/1‑69 (MQ=255)                                  gtCTTTGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAg                                >  1:2390159/2‑71 (MQ=9)                                   tCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCt                                                               <  1:179116/39‑1 (MQ=255)                                    cTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtt                                                              <  1:2327577/39‑1 (MQ=255)                                     tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCtttt                                                            >  1:3404931/1‑40 (MQ=255)                                     tttGCTGCCATCTTTCGGCGTGATGAAACCAAAACCttt                                                             >  1:658817/1‑39 (MQ=255)                                     tttGCTGCCATCATCAGGAGTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCa                             >  1:1965061/1‑71 (MQ=12)                                       tgctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAggg                                                      >  1:472922/1‑44 (MQ=255)                                       tgctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCa                                         >  1:43453/1‑57 (MQ=255)                                       tgctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGt                               <  1:3571953/67‑1 (MQ=255)                                       tgctgcCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                            >  1:566309/1‑70 (MQ=255)                                            ccATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                       <  1:1421616/70‑1 (MQ=255)                                               tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTc                                          <  1:2476967/48‑1 (MQ=255)                                               tCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAAc                                   >  1:3936028/1‑55 (MQ=255)                                                 tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCAc                                        >  1:3815115/1‑48 (MQ=255)                                                 tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCAt                            <  1:3776268/60‑1 (MQ=255)                                                 tttCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTa                     >  1:2373808/1‑67 (MQ=255)                                                   tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa                   <  1:282312/67‑1 (MQ=255)                                                   tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                <  1:1702742/70‑1 (MQ=255)                                                   tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                <  1:1131223/70‑1 (MQ=255)                                                   tCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACAt                <  1:1899011/70‑1 (MQ=255)                                                    cGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGa                   <  1:2707690/66‑1 (MQ=255)                                                     ggCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat             <  1:3668892/71‑1 (MQ=255)                                                      gagTGATAAAGCCAAAACCTTTATCTGCGTTAAACCATTTTACTAAACCAGTCATTTTATTAGACATAttt            <  1:3732265/69‑1 (MQ=255)                                                      gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat             <  1:2375025/70‑1 (MQ=255)                                                      gCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAat             >  1:1403659/1‑70 (MQ=255)                                                         tgatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                       <  1:1443865/57‑1 (MQ=255)                                                          gatgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAg                    <  1:141036/59‑1 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAACCCAGTCATTTTGTTAGACATAATTATTAc       >  1:1118278/1‑71 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTc                              >  1:2430700/1‑48 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc       >  1:1791914/1‑71 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc       >  1:3689580/1‑71 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc       >  1:681991/1‑71 (MQ=255)                                                           atgaAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTAc       >  1:2435792/1‑71 (MQ=255)                                                              aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTaaa                                    <  1:2146255/39‑1 (MQ=255)                                                              aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTa                                      <  1:1102343/37‑1 (MQ=255)                                                              aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGt                       >  1:3279076/1‑52 (MQ=255)                                                              aaaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtt    <  1:3966583/71‑1 (MQ=255)                                                               aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                  >  1:3125964/1‑56 (MQ=255)                                                               aaCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCttt   <  1:2234344/71‑1 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                  >  1:1521040/1‑55 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAc                  >  1:3450140/1‑55 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGAGATAATTATTACCtttt  >  1:3562023/1‑71 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAatt            <  1:275604/61‑1 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt  <  1:2667077/71‑1 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt  <  1:2858166/71‑1 (MQ=255)                                                                aCCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCtttt  <  1:1097479/71‑1 (MQ=255)                                                                |                                                                       GTTGCTCTGAATTGCTGAGAAATAGACAAACACATCTTTGCTGCCATCTTTCGGCGTGATGAAACCAAAACCTTTTTCAGGGTTAAACCATTTCACTAAACCAGTCATTTTGTTAGACATAATTATTACCTTTT  >  minE/980185‑980318

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