Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	822,285	C→T	intergenic (+149/‑18)	unknown → / → unknown	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE/Sodium‑dependent anion transporter family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	822,285	0	C	T	100.0%	93.6 / NA	36	intergenic (+149/‑18)	unknown/unknown	DNA‑binding domain of ModE / Molybdate‑binding domain of ModE/Sodium‑dependent anion transporter family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/16);  total (20/16)

GTGCTACACTCACCAAAACTCAACATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTCTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAACCGCCCTTTTT  >  1/822042‑822506                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              gtgCTACACTCACCAAAACTCAACATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTAttttt                                                                                                                                                                                                                        >  1:520987/1‑251 (MQ=255)                     cAACATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTacac                                                                                                                                                                                                    <  1:827406/251‑1 (MQ=255)                     cAACATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTaca                                                                                                                                                                                                     >  1:712782/1‑250 (MQ=255)                      aaCATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTAcacc                                                                                                                                                                                                   <  1:141322/251‑1 (MQ=255)                          tatTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACtt                                                                                                                                                                                               >  1:179822/1‑251 (MQ=255)                                 cttgagcAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTc                                                                                                                                                                                        <  1:557005/251‑1 (MQ=255)                                            ggTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACt                                                                                                                                                                                         <  1:894061/239‑1 (MQ=255)                                                  gATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCaa                                                                                                                                                                        <  1:571461/250‑1 (MQ=255)                                                                cGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAg                                                                                                                                                         <  1:395454/251‑1 (MQ=255)                                                                               ttGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGAt                                                                                                                                                                                                         <  1:95176/188‑1 (MQ=255)                                                                                    aCCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATg                                                                                                                                     >  1:612656/1‑251 (MQ=255)                                                                                         ttattaGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACATTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTTTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGCGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGAATGATTTTCGTTATTGATGTCAtt                                                                                                                                >  1:123510/1‑251 (MQ=255)                                                                                                                    agcTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGAttt                                                                                                                                                 <  1:641222/207‑1 (MQ=255)                                                                                                                     gcTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCtt                                                                                                     >  1:687518/1‑250 (MQ=255)                                                                                                                              aaTTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATGCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTtgatg                                                                                           <  1:425343/251‑1 (MQ=255)                                                                                                                                 tATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTATTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAAc                                                                                        <  1:867952/251‑1 (MQ=255)                                                                                                                                               cTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGt                                                                          >  1:335631/1‑251 (MQ=255)                                                                                                                                               cTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGt                                                                          >  1:256100/1‑251 (MQ=255)                                                                                                                                               cTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGt                                                                          >  1:224632/1‑251 (MQ=255)                                                                                                                                                    gCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTgg                                                                      >  1:394379/1‑250 (MQ=255)                                                                                                                                                         aCATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCattat                                                                >  1:617742/1‑251 (MQ=255)                                                                                                                                                                          cATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCAt                                                <  1:802841/250‑1 (MQ=255)                                                                                                                                                                                         gcgcATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAg                                >  1:76237/1‑251 (MQ=255)                                                                                                                                                                                                 tctctcAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATaaa                        <  1:539881/251‑1 (MQ=255)                                                                                                                                                                                                 tctctcAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATaaa                        <  1:768454/251‑1 (MQ=255)                                                                                                                                                                                                    ctcAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGt                     >  1:196803/1‑251 (MQ=255)                                                                                                                                                                                                        ggTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTc                                                                                                                       <  1:525877/149‑1 (MQ=255)                                                                                                                                                                                                        ggTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAAtg                                                                                                          <  1:371291/162‑1 (MQ=255)                                                                                                                                                                                                        ggTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAg                  >  1:411296/1‑250 (MQ=255)                                                                                                                                                                                                         gTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCg                >  1:864140/1‑251 (MQ=255)                                                                                                                                                                                                              ttCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAAc            >  1:172253/1‑250 (MQ=255)                                                                                                                                                                                                              ttCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAAc            >  1:120464/1‑250 (MQ=255)                                                                                                                                                                                                                   gttgtGACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAACCGcc        >  1:580873/1‑249 (MQ=255)                                                                                                                                                                                                                       tgACATATTATTTTATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAACCGCCCttttt  <  1:485193/251‑1 (MQ=255)                                                                                                                                                                                                                                 ttttATGATTTTTGTCTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATcg                                                                              >  1:175139/1‑165 (MQ=255)                                                                                                                                                                                                                                                cTTTTTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGaaa                                                                                                                                                           >  1:263656/1‑73 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                              GTGCTACACTCACCAAAACTCAACATATTGAGCTTGAGCAACAGGTTTGGATTTCCGTCGATCCGGAGCAAATCATTCTTGCTACCCTTTATTAGTCTTTTTTTCCATGAAGAGCAGCTACGATGAATTATCATTTGGAGCACTTTTGCCATACATGGCATAAAAGTCGCATTTTGTTGACGTAGCGCATAGTCTCTCAGGTATATTCCTGTTGTGACATATTATTTTATGATTTTTGTCTTTCTATTTTTAGGGTGATTTATGATTACACCACTTTTTACTCCATCACAGAAAAGCAAAATAAACAAGAAAAGTATGATTTTCGTTATTGATGTCATTCTTTTTGTCATTCTTCTTAATGTGCTTCCCTTTGATGAACAAGCCAATCGCGGTTTGGCATTATTGGCGTTTGTCGGCATTTTATGGTTAACCGAGGCCATAAATGTTAGCGTAACCGCCCTTTTT  >  1/822042‑822506

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