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%	103.1 / NA	33	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 (20/13);  total (20/13)

TGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCTGGAAA  >  NC_000913/1614658‑1614924                                                                                                                                    |                                                                                                                                        tGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                                  <  1:799803/139‑1 (MQ=255)  ggCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCgag                                                                                                                                 <  1:799151/139‑1 (MQ=255)                aaaTTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCt                                                                                                                   <  2:165128/139‑1 (MQ=255)                  aTTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGa                                                                                                                 >  1:403159/1‑139 (MQ=255)                            ttGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGga                                                                                                       <  2:668494/139‑1 (MQ=255)                              gCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGgaga                                                                                                     >  2:287833/1‑139 (MQ=255)                              gCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAACCGAGAATGCTa                                                                                                                          <  2:669160/118‑1 (MQ=255)                              gCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAACCGAGAATGCTa                                                                                                                          >  1:669160/1‑118 (MQ=255)                                            cATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGtt                                                                                       >  1:36000/1‑139 (MQ=255)                                             aTCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTg                                                                                      >  2:792820/1‑139 (MQ=255)                                              tCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGgag                                                                                                      >  1:330354/1‑122 (MQ=255)                                              tCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGgag                                                                                                      <  2:330354/122‑1 (MQ=255)                                               cGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGcc                                                                                    >  2:83313/1‑139 (MQ=255)                                                            ttATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagca                                                                       >  2:206099/1‑139 (MQ=255)                                                            ttATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagca                                                                       <  2:596881/139‑1 (MQ=255)                                                                  gTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAAccc                                                                 >  2:733977/1‑139 (MQ=255)                                                                                aTACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCgaga                                                                                                                                <  1:377947/62‑1 (MQ=255)                                                                                aTACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCgaga                                                                                                                                >  2:377947/1‑62 (MQ=255)                                                                                   cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTg                                                >  2:411390/1‑139 (MQ=255)                                                                                          ttACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATcgcg                                         >  2:429082/1‑139 (MQ=255)                                                                                            aCTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTg                                       <  2:403159/139‑1 (MQ=255)                                                                                             ctctTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGAAGTTCATCGCGGTCGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGc                                      >  1:351751/1‑139 (MQ=255)                                                                                               ctTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCg                                    <  1:411390/139‑1 (MQ=255)                                                                                                  gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTa                                                                                                                          <  1:218462/50‑1 (MQ=255)                                                                                                  gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTa                                                                                                                          >  2:218462/1‑50 (MQ=255)                                                                                                  gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg                                 >  1:395710/1‑139 (MQ=255)                                                                                                  gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg                                 >  2:212026/1‑139 (MQ=255)                                                                                                                       aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTgcgc            >  2:119572/1‑139 (MQ=255)                                                                                                                       aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTgcgc            >  1:266579/1‑139 (MQ=255)                                                                                                                            cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCt                                                                              >  2:233363/1‑68 (MQ=255)                                                                                                                            cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCt                                                                              <  1:233363/68‑1 (MQ=255)                                                                                                                            cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt       >  2:112476/1‑139 (MQ=255)                                                                                                                                 tCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTTCCCGTCTGCGCCAGCTGGaaa  <  2:351751/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        TGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCTGGAAA  >  NC_000913/1614658‑1614924

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