Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	538,165	C→T	15.4%	S280L (TCG→TTG)	rhlE →	RNA helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	538,165	0	C	T	15.4%	200.8 / 27.2	78	S280L (TCG→TTG)	rhlE	RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base T (9/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCGCTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACT  >  minE/538096‑538232                                                                      |                                                                    gcgcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCg                                                                    <  1:1324558/71‑1 (MQ=255)  cgcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGc                                                                   <  1:2578264/71‑1 (MQ=255)  cgcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGc                                                                   <  1:3255997/71‑1 (MQ=255)  cgcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGc                                                                   <  1:3475077/71‑1 (MQ=255)   gcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGc                                                                   >  1:2243410/1‑70 (MQ=255)   gcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCa                                                                  <  1:1174989/71‑1 (MQ=255)   gcTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCa                                                                  <  1:1608624/71‑1 (MQ=255)      aaCCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAgg                                                               >  1:3236068/1‑71 (MQ=255)      aaCCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAgg                                                               >  1:3487546/1‑71 (MQ=255)       aCCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGt                                                              <  1:2311813/71‑1 (MQ=255)       aCCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAg                                                                <  1:1273465/69‑1 (MQ=255)        ccATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTg                                                             <  1:2122477/71‑1 (MQ=255)         cATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTgc                                                            <  1:3534596/71‑1 (MQ=255)         cATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTgc                                                            <  1:3573632/71‑1 (MQ=255)         cATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTgc                                                            <  1:3578192/71‑1 (MQ=255)         cATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTg                                                             >  1:2878993/1‑70 (MQ=255)           tCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAg                                                                <  1:704797/65‑1 (MQ=255)           tCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTgcgc                                                          >  1:1771227/1‑71 (MQ=255)             tGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGTGcgcg                                                         >  1:2664533/1‑70 (MQ=255)             tGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGcgcg                                                         >  1:3508611/1‑70 (MQ=255)               gCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGTGCGCGTAc                                                      >  1:1842568/1‑71 (MQ=255)                 tGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTc                                                    >  1:3542610/1‑71 (MQ=255)                   aaCAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGTgcgc                                                          >  1:3364447/1‑63 (MQ=255)                   aaCAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGTGCGCGTACTCGt                                                  >  1:279809/1‑71 (MQ=255)                   aaCAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATTGCAAGGTGCGCGTACTCGt                                                  >  1:2135679/1‑71 (MQ=255)                    aCAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTg                                                 >  1:2500122/1‑71 (MQ=255)                    aCAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTg                                                 >  1:1008625/1‑71 (MQ=255)                      aGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTgcg                                               <  1:3801996/71‑1 (MQ=255)                        cTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCa                                                                  <  1:3387406/50‑1 (MQ=255)                        cTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCa                                                                  <  1:3630538/50‑1 (MQ=255)                        atCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGAGCGTACTCGTGCgct                                             >  1:1712199/2‑71 (MQ=255)                          cAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCa                                                                  >  1:711118/1‑48 (MQ=255)                           aaTAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTgc                                                            <  1:3508869/53‑1 (MQ=255)                            aTAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCgctggct                                         <  1:164878/71‑1 (MQ=255)                              aaaGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGa                                       >  1:3351502/1‑71 (MQ=255)                              aaaGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGa                                       >  1:183288/1‑71 (MQ=255)                               aaGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGAt                                      >  1:900126/1‑71 (MQ=255)                                     gCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTaaa                                <  1:3582435/71‑1 (MQ=255)                                        tCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCg                             <  1:1140324/71‑1 (MQ=255)                                        tCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCg                             <  1:1523292/71‑1 (MQ=255)                                         ccGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTgcgc                                              >  1:3415966/1‑53 (MQ=255)                                          cGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTaa                                 >  1:2842645/1‑65 (MQ=255)                                            tAGTGCGGCGATCCACGGCAATAAATGGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGc                          >  1:1452004/1‑70 (MQ=255)                                            tAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGc                          >  1:1690004/1‑70 (MQ=255)                                              gTGCGGCGATCCACGGCAATAAATTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGa                        <  1:3104667/70‑1 (MQ=255)                                               tGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGata                      >  1:2169333/1‑71 (MQ=255)                                               tGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTATAAATCGGGCGata                      <  1:2766400/71‑1 (MQ=255)                                                 cggcgATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGata                      >  1:1928869/1‑69 (MQ=255)                                                 aggcgATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATAtt                    >  1:4106265/2‑71 (MQ=255)                                                  ggcgATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGAt                                      <  1:3803125/52‑1 (MQ=255)                                                     gATCCACGGCAATAAATCGCAAGCTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCgtg                >  1:3074136/1‑71 (MQ=255)                                                      aTCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCgtg                >  1:2026926/1‑70 (MQ=255)                                                         cACGGCAATAAATTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACt            >  1:716885/1‑71 (MQ=255)                                                           cGGCAATAAATAGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtgg          >  1:918918/1‑71 (MQ=255)                                                            ggCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtggt         <  1:3593635/71‑1 (MQ=255)                                                            ggCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtggt         <  1:31352/71‑1 (MQ=255)                                                            ggCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtggt         <  1:862987/71‑1 (MQ=255)                                                             gCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtggtg        <  1:1893703/71‑1 (MQ=255)                                                              cAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACtggtg        <  1:993976/70‑1 (MQ=255)                                                               aaTAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCgg                            <  1:822109/49‑1 (MQ=255)                                                               aaTAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGc      <  1:1601819/71‑1 (MQ=255)                                                                 tAAATTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCaa    <  1:3530274/71‑1 (MQ=255)                                                                 tAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCgtg                <  1:1193268/59‑1 (MQ=255)                                                                 tAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCaa    >  1:2275795/1‑71 (MQ=255)                                                                  aaaTTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGata                      >  1:2255236/1‑52 (MQ=255)                                                                  aaaTTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGata                      >  1:2226398/1‑52 (MQ=255)                                                                  aaaTTGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCgt                 >  1:138728/1‑57 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCg                         >  1:3544853/1‑48 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  <  1:2946124/71‑1 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  >  1:765898/1‑71 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  >  1:627213/1‑71 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  <  1:3267342/71‑1 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  <  1:1798041/71‑1 (MQ=255)                                                                   aaTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  <  1:1011999/71‑1 (MQ=255)                                                                    aTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  >  1:3029261/1‑70 (MQ=255)                                                                    aTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  >  1:2889163/1‑70 (MQ=255)                                                                    aTCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACt  >  1:1665193/1‑70 (MQ=255)                                                                     tCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGatat                     >  1:3341948/1‑50 (MQ=255)                                                                      |                                                                    GCGCTAACCATCTGGCTGAACAGCTCAATAAAGATGGCATCCGTAGTGCGGCGATCCACGGCAATAAATCGCAAGGTGCGCGTACTCGTGCGCTGGCTGATTTTAAATCGGGCGATATTCGTGTACTGGTGGCAACT  >  minE/538096‑538232

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