Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	479,151	T→C	I16T (ATA→ACA)	unknown →	Sulfur carrier protein FdhD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	479,151	0	T	C	100.0%	99.3 / NA	36	I16T (ATA→ACA)	unknown	Sulfur carrier protein FdhD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

TTTTGATAGTTATAAATTCGTATCTGAGTAGATAATTTACGATCGGTTTTGTTAAAATCACATGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAATAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAGTATTGATATTGAAGAAGTTTGTAAT  >  1/478906‑479371                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             ttttGATAGTTATAAATTCGTATCTGAGTAGATAATTTACGATCGGTTTTGTTAAAATCACATGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATg                                                                                                                                                                                                                         <  1:811957/251‑1 (MQ=255)                    gTATCTGAGTAGATAATTTACGATCGGTTTTGTTAAAATCACATGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAaca                                                                                                                                                                                                       >  1:534560/1‑250 (MQ=255)                          gAGTAGATAATTTACGATCGGTTTTGTTAAAATCACATGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAatcat                                                                                                                                                                                                >  1:757498/1‑251 (MQ=255)                                                             caTGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAgg                                                                                                                                                             <  1:365517/251‑1 (MQ=255)                                                                                     cTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTCTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAAt                                                                                                                                      <  1:899099/250‑1 (MQ=255)                                                                                          taaatGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATaa                                                                                                                                       >  1:539641/1‑244 (MQ=255)                                                                                                         gTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCAcc                                                                                                                                                                                            <  1:558116/176‑1 (MQ=255)                                                                                                         gTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAatg                                                                                                                 >  1:288747/1‑251 (MQ=255)                                                                                                                    gTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGata                                                                                                                                  >  1:832594/1‑223 (MQ=255)                                                                                                                    gTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTc                                                                                                                              >  1:784261/1‑227 (MQ=255)                                                                                                                                        gggTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATaa                                                                                                                                       >  1:263931/1‑198 (MQ=255)                                                                                                                                        gggTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGa                                                                                                                                    >  1:20889/1‑201 (MQ=255)                                                                                                                                                  ccTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAAc                                                                                                                                                                                  <  1:864383/145‑1 (MQ=255)                                                                                                                                                   cTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAgg                                                                                        >  1:302751/1‑234 (MQ=255)                                                                                                                                                            cTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTcc                                                                                                                             >  1:313320/1‑188 (MQ=255)                                                                                                                                                                    aTCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCat                                                      >  1:122838/1‑251 (MQ=255)                                                                                                                                                                        ataAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGTTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCattatt                                                  >  1:400578/1‑251 (MQ=255)                                                                                                                                                                           aGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGaa                                               >  1:503298/1‑251 (MQ=255)                                                                                                                                                                            gTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAag                                                                                                                                                                <  1:750282/137‑1 (MQ=255)                                                                                                                                                                            gTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGaaa                                              >  1:581258/1‑251 (MQ=255)                                                                                                                                                                            gTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGaaa                                              >  1:224932/1‑251 (MQ=255)                                                                                                                                                                             taataaTAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACgg                                                                                                                      <  1:704125/178‑1 (MQ=255)                                                                                                                                                                             taataaTAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGAt                                                                                      <  1:642162/210‑1 (MQ=255)                                                                                                                                                                             taataaTAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGaaaa                                             >  1:304888/1‑251 (MQ=255)                                                                                                                                                                                   taGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGACAAACCGAg                                       >  1:59380/1‑251 (MQ=255)                                                                                                                                                                                     gTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAatgatg                                                                                                              >  1:105614/1‑178 (MQ=255)                                                                                                                                                                                     gTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGc                                                                   <  1:680656/221‑1 (MQ=255)                                                                                                                                                                                     gTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTg                                     >  1:74806/1‑251 (MQ=255)                                                                                                                                                                                             gggATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTtat                             >  1:594036/1‑251 (MQ=255)                                                                                                                                                                                              ggATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTtata                            >  1:141912/1‑251 (MQ=255)                                                                                                                                                                                              ggATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTtat                             <  1:225779/250‑1 (MQ=255)                                                                                                                                                                                               gATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGG‑TTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAg                           <  1:697352/250‑1 (MQ=255)                                                                                                                                                                                                                 aTCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAGTATTGATATTGaagaag          <  1:70583/250‑1 (MQ=255)                                                                                                                                                                                                                     caAAAAAAGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAatg                                                                                                                 <  1:641136/143‑1 (MQ=255)                                                                                                                                                                                                                       aaaaaaGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAGTATTGATATTGAAGAAGTTTGTaa   <  1:441903/251‑1 (MQ=255)                                                                                                                                                                                                                        aaaaaGTAATTTTTTTTATCAAAAAATTAACAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAGTATTGATATTGAAGAAGTTTGTAAt  <  1:657398/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             TTTTGATAGTTATAAATTCGTATCTGAGTAGATAATTTACGATCGGTTTTGTTAAAATCACATGATATTTAAGCAAAATCGAAGCTAAATAAATGATGTGATATGTTAAATATCTGTTATACATTTGTTAGCGGTGGGTAGAATACCTATTTTTTCTACGGAAATCTATAAGTAATAATAGTTATAATGGGATAAAATTTTGAACTGGATCACAAAAAAAGTAATTTTTTTTATCAAAAAATTAATAAATGAGAAACAAAATCAAAACACAATCATCACCTCCGAGAAACAAGAAACCCTCGCAAAAGAGGTGCCTGTGGCGCTGGTTTATAATGATATTTCCCATACGGTAATGATGTGTTCGCCCCGTGATTTAGAGGATTTCGCTATGGGATTTTCGCTATCGGAAGGCATTATTGAAAAACCGAGTGATATTTATAGTATTGATATTGAAGAAGTTTGTAAT  >  1/478906‑479371

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