Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	571,521	A→G	intergenic (‑17/‑283)	unknown ← / → unknown	5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine methyltransferase (EC 2.1.1.14)/Transcriptional activator MetR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	571,521	0	A	G	100.0%	101.3 / NA	36	intergenic (‑17/‑283)	unknown/unknown	5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine methyltransferase (EC 2.1.1.14)/Transcriptional activator MetR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (11/25);  total (11/25)

TGTTGCCACTTGCAAGTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTATTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACGACCGCACTTTTTCATCGATTTCTTTATGGATTTTATGGATTTCTTCTATGTAACGTTAAAGTGCGAGCCGTTTT  >  1/571277‑571754                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          tGTTGCCACTTGCAAGTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAAt                                                                                                                                                                                                                                     <  1:690767/251‑1 (MQ=255)   ttGCCACTTGCAAGTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATg                                                                                                                                                                                                                                    <  1:391720/250‑1 (MQ=255)          ttGCAAGTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGAc                                                                                                                                                                                                                            <  1:437116/251‑1 (MQ=255)                gTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTaa                                                                                                                                                                                                                      >  1:324646/1‑251 (MQ=255)                 tCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTaaa                                                                                                                                                                                                                     >  1:292705/1‑251 (MQ=255)                   tAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACg                                                                                                                                                                                                                   <  1:261488/251‑1 (MQ=255)                          tgATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATAc                                                                                                                                                                                                            <  1:837159/251‑1 (MQ=255)                            aTCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCGCCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATGCTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACg                                                                                                                                                                                                           <  1:317004/250‑1 (MQ=255)                            aTCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACg                                                                                                                                                                                                           <  1:528669/250‑1 (MQ=255)                                tAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTggg                                                                                                                                                                                                       <  1:684404/250‑1 (MQ=255)                                tAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGc                                                                                                                                                                                                      <  1:801085/251‑1 (MQ=255)                                                           tAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGtattat                                                                                                                                                                            <  1:146030/250‑1 (MQ=255)                                                            aaTCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGc                                                                                                                                                                          <  1:213974/251‑1 (MQ=255)                                                                      cGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATtaataa                                                                                                                                                                <  1:601666/251‑1 (MQ=255)                                                                      cGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATtaataa                                                                                                                                                                <  1:603875/251‑1 (MQ=255)                                                                      cGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATtaataa                                                                                                                                                                <  1:67893/251‑1 (MQ=255)                                                                            cGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGtat                                                                                                                                                                               >  1:38645/1‑230 (MQ=255)                                                                                            ttGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATgaa                                                                                                                                             <  1:584881/248‑1 (MQ=255)                                                                                                    gcgTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAg                                                                                                                                                     <  1:396057/232‑1 (MQ=255)                                                                                                     cgTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTaaa                                                                                                                                                                                                                                      <  1:390415/150‑1 (MQ=255)                                                                                                          ggCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGa                                                                                                                                                                                                                                   <  1:440192/148‑1 (MQ=255)                                                                                                          ggCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGt                                                                                                                                                                                                                  <  1:470259/165‑1 (MQ=255)                                                                                                          ggCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACg                                                                                                                                                                                                           <  1:781403/172‑1 (MQ=255)                                                                                                           gCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAAt                                                                                                                                                                                                                                     >  1:76963/1‑145 (MQ=255)                                                                                                                    aTGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATtaataa                                                                                                                                                                >  1:839882/1‑205 (MQ=255)                                                                                                                      gTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGAc                                                                                                                                                                                                                            >  1:475937/1‑143 (MQ=255)                                                                                                                      gTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGtatta                                                                                                                                                                             >  1:813943/1‑190 (MQ=255)                                                                                                                                      tCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGc                                                                                                                                                                          >  1:763996/1‑177 (MQ=255)                                                                                                                                          cGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATaaaa                                                                                            <  1:369229/251‑1 (MQ=255)                                                                                                                                            aTTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATaaaaa                                                                                           <  1:350790/250‑1 (MQ=255)                                                                                                                                                          aTAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACg                                                                            >  1:539566/1‑251 (MQ=255)                                                                                                                                                                ttCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGtatta                                                                                                                                                                             >  1:83446/1‑148 (MQ=255)                                                                                                                                                                                                         gAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACGACCGCACTTTTTCATCGATTTCTTTATGGATTTTATGGATttcttc                              <  1:142558/250‑1 (MQ=255)                                                                                                                                                                                                                             gTTGTCATACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAAt                                                                                                                                                                                                <  1:97629/68‑1 (MQ=255)                                                                                                                                                                                                                                    tACTAATTCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACGACCGCACTTTTTCATCGATTTCTTTATGGATTTTATGGATTTCTTCTATGTAACGTTAAAGTGCGGGTCGtttt  >  1:218830/1‑251 (MQ=255)                                                                                                                                                                                                                                           tCCTCAATGTGTTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACGACCGCACTTTTTCATCGATTTCTTTATGGATTTTATGGATTTCTTCTATGTAACGTTAAAGTACGGGTCGtttt  <  1:401275/244‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          TGTTGCCACTTGCAAGTCTAAAATGTGATCGTAGAAGGTGAAATCCGCCACGGCGACGTAATCCGCATTCGCTGCCGCTTGGTGTTTCCAGTTGATTTCGCGTAAGGCTTTGGCGATGTCTAATAAATCGGCTTCTGCGATTTCACCGCGCCAATAACGTTCTTGTGCAAATTTAAGTTCACGTTTTGCGCCCACGCGCGGAAAGCCGGATAAATGAAAAGTTGTCATACTAATTCCTCAATGTATTAAATGATTTAGACGGCTAAACGTCTATACGTGGGCTAGAATGACGAAAACGGGGGTATTATGCAAATTAATAATTTTCTAAAAGTTTATGAAAGAAATTAATGATCCTGTATTTTATTAGTTTGCGGTGGTAATCATAAAAAACGCTTTTAAAAACGACCGCACTTTTTCATCGATTTCTTTATGGATTTTATGGATTTCTTCTATGTAACGTTAAAGTGCGAGCCGTTTT  >  1/571277‑571754

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