Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	992,086	A→G	intergenic (‑719/+142)	unknown ← / ← unknown	hypothetical protein/Thiol:disulfide oxidoreductase associated with MetSO reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	992,086	0	A	G	100.0%	98.3 / NA	38	intergenic (‑719/+142)	unknown/unknown	hypothetical protein/Thiol:disulfide oxidoreductase associated with MetSO reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/17);  total (21/17)

ATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCACGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCAAATCCGTATCAAGCAATACGCTGGAAGGGTA  >  1/991850‑992326                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 atatTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAatat                                                                                                                                                                                                                                    <  1:365690/251‑1 (MQ=255)                      atGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATg                                                                                                                                                                                                                  <  1:649125/248‑1 (MQ=255)                                                                                            aTTTGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCaaaaa                                                                                                                                            <  1:102065/248‑1 (MQ=255)                                                                                              ttGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTTCCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGAt                                                                                                                                                                                                                             <  1:113629/165‑1 (MQ=255)                                                                                                           tAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGAc                                                                                                                           <  1:780113/250‑1 (MQ=255)                                                                                                                           gTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGcac                                                                                                                                        >  1:409031/1‑221 (MQ=255)                                                                                                                                                tAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAg                                                                                      >  1:796884/1‑250 (MQ=255)                                                                                                                                                 aGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGcc                                                                                                                                                                                                             >  1:880318/1‑130 (MQ=255)                                                                                                                                                     taACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGc                                                                                     >  1:803869/1‑246 (MQ=255)                                                                                                                                                       aCTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAatat                                                                                                                                                                                                                                    >  1:772214/1‑101 (MQ=255)                                                                                                                                                       aCTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATgg                                                                                                                                   >  1:804784/1‑198 (MQ=255)                                                                                                                                                       aCTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGaattaat                                                                                                >  1:788438/1‑233 (MQ=255)                                                                                                                                                       aCTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTaa                                                                              >  1:574950/1‑251 (MQ=255)                                                                                                                                                        cTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAAttt                                                                                              >  1:867714/1‑234 (MQ=255)                                                                                                                                                                 aTTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATc                                                                                                                                                                                                                     >  1:545965/1‑106 (MQ=255)                                                                                                                                                                         aTGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTgg                                                                                                                                                                                                                               <  1:522494/88‑1 (MQ=255)                                                                                                                                                                         aTGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGAt                                                                                                                                                                                                                             >  1:154849/1‑90 (MQ=255)                                                                                                                                                                              ggaggTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTcc                                                        >  1:198437/1‑250 (MQ=255)                                                                                                                                                                                 ggTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGAt                                                                                                                                                                                                                             <  1:622077/82‑1 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAatat                                                                                                                                                                                                                                    >  1:367241/1‑74 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGAt                                                                                                                                                                                                                             >  1:877435/1‑81 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATg                                                                                                                                                                                                                  <  1:505038/92‑1 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGcc                                                                                                                                                                                                             <  1:348748/97‑1 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATAtt                                                                                                                                                                                               <  1:464222/111‑1 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTa                                                   >  1:200125/1‑251 (MQ=255)                                                                                                                                                                                  gTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTa                                                   >  1:358727/1‑251 (MQ=255)                                                                                                                                                                                    ataGTACCATCTCAATATTTGCATGAATCTAAAGACGCAAATTATGCAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAacac                                                                                                                                                                                          >  1:394334/1‑114 (MQ=255)                                                                                                                                                                                    ataGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGata                                                                                                                                                                                                   <  1:494339/105‑1 (MQ=255)                                                                                                                                                                                       gTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAgg                                                                                                       >  1:599695/1‑194 (MQ=255)                                                                                                                                                                                                    atTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAatat                                                                                                                                                                                                                                    >  1:222376/1‑56 (MQ=255)                                                                                                                                                                                                         cATGAATCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCAAATc                             <  1:740286/250‑1 (MQ=255)                                                                                                                                                                                                             aaTCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGata                                                                                                                                                                                                   <  1:776630/80‑1 (MQ=255)                                                                                                                                                                                                             aaTCTAAAGACGAAAATTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCAAATCCGTAt                        <  1:881610/251‑1 (MQ=255)                                                                                                                                                                                                                         aaaaTTATGTAAATAAAACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCa                                 <  1:20905/230‑1 (MQ=255)                                                                                                                                                                                                                                   aaataaaACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTgcgc                                                                       <  1:748764/182‑1 (MQ=255)                                                                                                                                                                                                                                   aaataaaACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCAAATCCGTATCAAGCAATACGCTGGAAGGGTa  >  1:389630/1‑251 (MQ=255)                                                                                                                                                                                                                                    aataaaACCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAgg                                                                                                       >  1:101986/1‑149 (MQ=255)                                                                                                                                                                                                                                         aaCCGCGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTaaa                                                             <  1:410219/186‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                 ATATTAACTGTTGAAGATAATATGATAATGAGTTATCTAACTATATGATGATATACCAAAATGGTTTATTTTATTTTGATGTAGAAATTGAATTTGATATAGTGCGTAGATAACAAATAAAAGTTCTTTATCAAATAAAAAATTAGTATAACTATACCTTATTTCAAAATGAAGGAGGTATAGTACCATCTCAATATTTGCATGAATCTAAAGACGAAAATTATGTAAATAAAACCACGACTTTACAATATTGTGGATTACTTATCATGTTGCCTGCTTAGATATATTAACACTCTAAAAAATGTTATTTTGGGGTGCAGATCAAAGTGCGGTCAAAAAGCACAATGGATTTTGACCGCACTTTTTCGTATTAAGGAATTAATTTCCTGCCAGCTGTTTAATTTGCGCCGTGCCTAAATGTCCCGGTACGGTTTTTTGGATATTCAAATCCGTATCAAGCAATACGCTGGAAGGGTA  >  1/991850‑992326

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑