Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	988,133	C→G	11.1%	L792V (CTG→GTG)	ynfE →	oxidoreductase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	988,133	0	C	G	11.1%	181.6 / 11.2	63	L792V (CTG→GTG)	ynfE	oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base G (4/3);  total (35/28)
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.84e-01

AATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATA  >  Exported/988064‑988201                                                                      |                                                                     aaTATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACt                                                                     <  1:1216463/71‑1 (MQ=255)     tGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGc                                                                  <  1:2258357/70‑1 (MQ=255)     tGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCg                                                                 <  1:3391165/71‑1 (MQ=255)       tCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCg                                                                 >  1:1325714/1‑69 (MQ=255)       tCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGa                                                                <  1:858155/70‑1 (MQ=255)            cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                          <  1:2818790/71‑1 (MQ=255)            cGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                          <  1:2699108/71‑1 (MQ=255)                 aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                     <  1:1239386/71‑1 (MQ=255)                 aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                     >  1:1221467/1‑71 (MQ=255)                 aaaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                     >  1:2635735/1‑71 (MQ=255)                  aaTCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCg                                                     <  1:151833/70‑1 (MQ=255)                             gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTaa                                          >  1:402700/1‑70 (MQ=255)                             gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:1474647/1‑71 (MQ=255)                             gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:1501389/1‑71 (MQ=255)                             gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:1840596/1‑71 (MQ=255)                             gcggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:183870/1‑71 (MQ=255)                              cggcTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         <  1:407268/70‑1 (MQ=255)                                 cTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTgg                                     >  1:2683691/1‑71 (MQ=255)                                  tgtgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTa                                           <  1:1443779/64‑1 (MQ=255)                                  tgtgtgAATACGCTAACCACTCTGCGACCTTCACCACt                                                                     <  1:320147/38‑1 (MQ=255)                                    tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:1096262/1‑64 (MQ=255)                                    tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:1398391/1‑64 (MQ=255)                                    tgtgAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                   <  1:187340/70‑1 (MQ=255)                                       gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                          <  1:3130649/44‑1 (MQ=255)                                       gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                          <  1:1682497/44‑1 (MQ=255)                                       gAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGaa                                                          <  1:2379727/44‑1 (MQ=255)                                             gCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaa                          >  1:752801/1‑70 (MQ=255)                                              cTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         >  1:34722/1‑70 (MQ=255)                                              cTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         >  1:2829146/1‑70 (MQ=255)                                              cTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         >  1:2632320/1‑70 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgaga                               >  1:1753739/1‑61 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgaga                               >  1:1679553/1‑61 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCgag                                >  1:2270126/1‑60 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                           <  1:984386/65‑1 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGa                           <  1:2491769/65‑1 (MQ=255)                                                 aCCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAat                     <  1:1409491/71‑1 (MQ=255)                                                   cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAatat                   <  1:3014067/71‑1 (MQ=255)                                                   cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAatat                   <  1:3082596/71‑1 (MQ=255)                                                   cACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAata                    <  1:938746/70‑1 (MQ=255)                                                    aCTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAatat                   >  1:3415244/1‑70 (MQ=255)                                                       ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                   >  1:1805583/1‑51 (MQ=255)                                                       ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                   >  1:601571/1‑51 (MQ=255)                                                       ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                   >  1:3194947/1‑51 (MQ=255)                                                       ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTc                                   >  1:3132404/1‑51 (MQ=255)                                                       ctGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc               >  1:2093287/1‑71 (MQ=255)                                                          cgcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAAcc               >  1:1174709/1‑68 (MQ=255)                                                          cgcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCAc            <  1:3730349/71‑1 (MQ=255)                                                           tcccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         <  1:747772/56‑1 (MQ=25)                                                           gcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAat                     >  1:2092111/1‑61 (MQ=255)                                                           gcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAat                     >  1:729636/1‑61 (MQ=255)                                                           gcCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAat                     >  1:2631930/1‑61 (MQ=255)                                                            cccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         <  1:1291168/56‑1 (MQ=37)                                                            cccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaa                         <  1:1215909/56‑1 (MQ=37)                                                             ccTTCGCCGGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       >  1:1040233/1‑57 (MQ=21)                                                             ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       >  1:2372890/1‑57 (MQ=37)                                                             ccTTCGCCAGTGGTGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       >  1:2219659/1‑57 (MQ=37)                                                             ccTTCGCCAGTGATGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       >  1:1638431/1‑57 (MQ=25)                                                             ccTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGa          >  1:3396006/1‑70 (MQ=255)                                                              cTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGAc         >  1:3181267/1‑70 (MQ=255)                                                               ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       <  1:95858/55‑1 (MQ=255)                                                               ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaaa                       <  1:146156/55‑1 (MQ=255)                                                               ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaa                          <  1:3175784/52‑1 (MQ=255)                                                               ttCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACaa       >  1:3740456/1‑71 (MQ=255)                                                                 caccacTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:2989082/1‑35 (MQ=255)                                                                 caccacTGGCGAAGGGAAACCCGCAGCACACTAAt                                         >  1:613182/1‑35 (MQ=255)                                                                 caccacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAAc      >  1:3820309/1‑70 (MQ=255)                                                                  accacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCa    >  1:2750904/1‑71 (MQ=255)                                                                    cacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGaaaa                        <  1:3236255/49‑1 (MQ=255)                                                                    cacTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATa  <  1:611858/71‑1 (MQ=255)                                                                     acTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATa  >  1:202251/1‑70 (MQ=255)                                                                      |                                                                     AATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCACTGGCGAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATA  >  Exported/988064‑988201

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