Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,614,789	C→T	intergenic (‑86/‑15)	sad ← / → yneJ	succinate semialdehyde dehydrogenase, NAD(P)+‑dependent/putative DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,614,789	0	C	T	100.0%	146.4 / NA	46	intergenic (‑86/‑15)	sad/yneJ	succinate semialdehyde dehydrogenase, NAD(P)+‑dependent/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/22);  total (24/22)

GGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCT  >  NC_000913/1614659‑1614919                                                                                                                                   |                                                                                                                                   ggCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCGACTTATTCATGAATCgag                                                                                                                            <  1:784770/139‑1 (MQ=255)           aTCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATgg                                                                                                                  >  1:820281/1‑139 (MQ=255)                aaTTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          >  1:511392/1‑126 (MQ=255)                aaTTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          <  2:511392/126‑1 (MQ=255)                         aGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGTGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           <  1:427885/132‑1 (MQ=255)                         aGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGTGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           >  2:427885/1‑132 (MQ=255)                               cGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATg                                                                                              >  1:903/1‑139 (MQ=255)                                        ggTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAcgcg                                                                                     <  2:820281/139‑1 (MQ=255)                                            aTCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          <  2:312806/98‑1 (MQ=255)                                            aTCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          >  1:312806/1‑98 (MQ=255)                                            aTCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTg                                                                                 >  2:181952/1‑139 (MQ=255)                                                                  tCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATTAATAATATTGAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCa                                                           >  2:662172/1‑139 (MQ=255)                                                                   cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagca                                                                  >  2:1036961/1‑131 (MQ=255)                                                                   cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagca                                                                  <  1:1036961/131‑1 (MQ=255)                                                                     tGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          >  1:769562/1‑73 (MQ=255)                                                                     tGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          <  2:769562/73‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAAt                                                                                                                                <  1:914093/54‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAAt                                                                                                                                >  2:914093/1‑54 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTa                                                                                                                     >  1:578596/1‑65 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTa                                                                                                                     <  2:578596/65‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTAt                                                                                                                    <  1:1006835/66‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTAt                                                                                                                    >  2:1006835/1‑66 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGa                                                                                                            <  1:174530/74‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGa                                                                                                            >  2:174530/1‑74 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGaa                                                                            >  1:78214/1‑106 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGaa                                                                            <  2:78214/106‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTg                                                                        <  1:1015131/110‑1 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTg                                                                        >  2:1015131/1‑110 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAgag                                           >  2:974453/1‑137 (MQ=255)                                                                                  cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTg                                           >  1:204717/1‑139 (MQ=255)                                                                                        tttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATcgc                                     <  1:181952/139‑1 (MQ=255)                                                                                              ctTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCg                               <  1:442181/139‑1 (MQ=255)                                                                                                 gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg                            <  2:903/139‑1 (MQ=255)                                                                                                  cTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCGACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGa                           <  1:662172/139‑1 (MQ=255)                                                                                                    ttAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAt                         >  1:798760/1‑139 (MQ=255)                                                                                                            aataatATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTAc                 >  2:485987/1‑139 (MQ=255)                                                                                                                      aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           >  1:816103/1‑39 (MQ=255)                                                                                                                      aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           <  2:816103/39‑1 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAg    <  2:557630/137‑1 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAg    >  1:557630/1‑137 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt  >  2:994222/1‑139 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt  >  2:473225/1‑139 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt  <  1:994222/139‑1 (MQ=255)                                                                                                                           cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt  <  1:473225/139‑1 (MQ=255)                                                                                                                            ttATTCATGAATCGAGAAGGCTATGGATCTGACCCAACTGGAGATGTTCAACGCgg                                                                                    >  1:707347/1‑56 (MQ=255)                                                                                                                            tcaTTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCgg                                                                                    <  2:707347/54‑1 (MQ=255)                                                                                                                                   |                                                                                                                                   GGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCT  >  NC_000913/1614659‑1614919

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