Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	1,074,318	G→T	8.7%	A39E (GCG→GAG)	ycdM ←	predicted monooxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	1,074,318	0	G	T	8.7%	186.8 / 2.9	57	A39E (GCG→GAG)	ycdM	predicted monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (21/31);  new base T (2/3);  total (23/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.97e-01

TCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCT  >  AP009048/1074249‑1074385                                                                      |                                                                    tcgaaatGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGgcgc                                                                    <  1:2392253/71‑1 (MQ=255)      atGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTggg                                                               >  1:4363561/1‑71 (MQ=255)       tGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTggg                                                               <  1:1437670/70‑1 (MQ=255)          tAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGa                                                           <  1:1425164/71‑1 (MQ=255)          tAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGa                                                           <  1:4209646/71‑1 (MQ=255)            gTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGaa                                                          <  1:2998646/70‑1 (MQ=255)             tGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGACGCGTGGGTCGAAAt                                                        <  1:3609670/71‑1 (MQ=255)               cTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATg                                                       <  1:4572624/70‑1 (MQ=255)                tCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTgg                                                                <  1:264928/60‑1 (MQ=255)                tCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGt                                                              >  1:2691814/1‑62 (MQ=255)                 ccGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATg                                                       >  1:3935963/1‑68 (MQ=255)                  cGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATg                                                       <  1:1752000/67‑1 (MQ=255)                  cGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagc                                                    >  1:1903320/1‑70 (MQ=255)                   gCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGa                                                      >  1:1515496/1‑67 (MQ=255)                        ttGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgt                                             >  1:669410/1‑71 (MQ=255)                          gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                            <  1:990390/70‑1 (MQ=255)                          gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                            <  1:4566625/70‑1 (MQ=255)                          gCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                            <  1:1530703/70‑1 (MQ=255)                           cacaATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttgt                                          <  1:2575368/71‑1 (MQ=255)                           cacaATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                           <  1:3612660/70‑1 (MQ=255)                           cacaATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                           <  1:119367/70‑1 (MQ=255)                            acaATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTgg                                                                >  1:1226184/1‑48 (MQ=255)                              aaTGGCTTTATTCAGTTCAACGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGc                                       <  1:3991488/71‑1 (MQ=255)                              aaTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGa                                                           <  1:3362233/51‑1 (MQ=255)                               aTGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                           <  1:694686/66‑1 (MQ=255)                                  gCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCaa                                    <  1:1482451/70‑1 (MQ=255)                                  gCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAAt                                   <  1:2272631/71‑1 (MQ=255)                                  gCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAAt                                   <  1:3902942/71‑1 (MQ=255)                                  gCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAAt                                   <  1:3911682/71‑1 (MQ=255)                                    tttATTCAGTTCAAAGGTCGGCATGTACTGCGGcgcg                                                                   <  1:4419610/37‑1 (MQ=255)                                     ttATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAAt                                                        <  1:4478678/47‑1 (MQ=255)                                      tATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAgg                                >  1:3321662/1‑70 (MQ=255)                                      tATTCAGTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAg                                 <  1:2358821/69‑1 (MQ=255)                                      tATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAgg                                >  1:1601217/1‑70 (MQ=255)                                       aTTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGt                               <  1:3034979/70‑1 (MQ=255)                                       aTTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTa                              >  1:3825357/1‑71 (MQ=255)                                       aTTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTa                              >  1:115163/1‑71 (MQ=255)                                        ttCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTAc                             <  1:1513692/71‑1 (MQ=255)                                         tCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgtt                                            <  1:625920/55‑1 (MQ=255)                                           aGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGt                                                              <  1:908679/35‑1 (MQ=255)                                           aGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGaa                          <  1:2032728/71‑1 (MQ=255)                                           aGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGGCGAAATGAGCCAGCCGTTGTTGCCAATAgg                                <  1:3228937/65‑1 (MQ=255)                                            gTTCAAAGGTCGGCATGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAAt                         <  1:2042128/71‑1 (MQ=255)                                             ttCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATa                        >  1:4607881/1‑71 (MQ=255)                                               cAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACg                      >  1:850587/1‑71 (MQ=255)                                               cAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACg                      >  1:1204988/1‑71 (MQ=255)                                                aaaGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagc                                                    >  1:3375628/1‑40 (MQ=255)                                                aaaGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGc                     <  1:2327938/71‑1 (MQ=255)                                                 aaGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGagccagc                                                >  1:2826281/1‑43 (MQ=255)                                                  aGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGcc                    >  1:3905249/1‑70 (MQ=255)                                                  aGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGcc                    >  1:1764512/1‑70 (MQ=255)                                                      cGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCgttg                                           <  1:3230192/43‑1 (MQ=255)                                                       ggCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAAtttt              >  1:1989939/1‑71 (MQ=255)                                                        gCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtc             <  1:3581539/71‑1 (MQ=255)                                                           tGTACTGCGGCTCGTGGGTCGAAATGAGCCAGCCg                                              >  1:16790/1‑35 (MQ=255)                                                           tGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcat           >  1:1175043/1‑70 (MQ=255)                                                            gTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatca         >  1:3228625/1‑71 (MQ=255)                                                             tACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATAc                       >  1:2545783/1‑56 (MQ=255)                                                               cTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTtcatcatca      <  1:273388/71‑1 (MQ=255)                                                                 gcggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGaa                          >  1:4320853/1‑49 (MQ=255)                                                                 gcggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGa                           >  1:4349650/1‑48 (MQ=255)                                                                 gcggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGa                           >  1:1612329/1‑48 (MQ=255)                                                                  cggcgCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAAcc   >  1:822921/1‑71 (MQ=255)                                                                    gcgcGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACct  <  1:584863/70‑1 (MQ=255)                                                                      |                                                                    TCGAAATGGTAGTGCTCCGCTTTTTGCACAATGGCTTTATTCAGTTCAAAGGTCGGCATGTACTGCGGCGCGTGGGTCGAAATGAGCCAGCCGTTGTTGCCAATAGGTACGAATACGCCAATTTTCATCATCAACCT  >  AP009048/1074249‑1074385

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