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%	164.6 / NA	51	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 (26/25);  total (26/25)

ACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCT  >  NC_000913/1614653‑1614919                                                                                                                                         |                                                                                                                                   aCCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGa                                                                                                                                  <  2:925449/139‑1 (MQ=255) aCCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGa                                                                                                                                  >  1:471821/1‑139 (MQ=255)             aTTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             <  2:126809/132‑1 (MQ=255)             aTTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             >  1:126809/1‑132 (MQ=255)                                   gCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGgaga                                                                                                >  1:1121695/1‑139 (MQ=255)                                    ccGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          <  1:1153054/112‑1 (MQ=255)                                    ccGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGaa                                                                                                                          >  2:1153054/1‑112 (MQ=255)                                            aTGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAcg                                                                                       >  1:544918/1‑139 (MQ=255)                                              ggTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAcgcg                                                                                     <  2:544918/139‑1 (MQ=255)                                                          aTCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCt                                                                         >  1:278513/1‑139 (MQ=255)                                                                tttATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagc                                                                   <  1:289580/139‑1 (MQ=255)                                                                      aGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCTACGCGGTTGCCGAAGCTGGCAGCATAAcc                                                             >  1:644559/1‑139 (MQ=255)                                                                         cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATc                                                                                                                               >  2:373919/1‑70 (MQ=255)                                                                         cATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATc                                                                                                                               <  1:373919/70‑1 (MQ=255)                                                                           tGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGc                                                        <  2:1121695/139‑1 (MQ=255)                                                                                 gAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagc                                                                   >  2:1398682/1‑122 (MQ=255)                                                                                 gAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagc                                                                   <  1:1398682/122‑1 (MQ=255)                                                                                      tACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgcag                                                   <  1:1244337/133‑1 (MQ=255)                                                                                      tACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTgcag                                                   >  2:1244337/1‑133 (MQ=255)                                                                                        cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcag                                                                    <  2:92491/114‑1 (MQ=255)                                                                                        cGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcag                                                                    >  1:92491/1‑114 (MQ=255)                                                                                            gATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATc                                       >  1:701489/1‑139 (MQ=255)                                                                                                 aCTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTg                                  <  2:278513/139‑1 (MQ=255)                                                                                                   tctTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTg                                                                                 <  2:1156047/90‑1 (MQ=255)                                                                                                   tctTGCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTg                                                                                 >  1:1156047/1‑90 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             <  1:576782/42‑1 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             <  1:1010865/42‑1 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             >  2:576782/1‑42 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCga                                                                                                                             >  2:1010865/1‑42 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           >  2:1162536/1‑60 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGAc                                                                                                           <  1:1162536/60‑1 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTTCCGAAGCTGgcagc                                                                   >  1:717275/1‑100 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGgcagc                                                                   <  2:717275/100‑1 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg                            >  1:1232464/1‑139 (MQ=255)                                                                                                       gCTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCg                            >  1:210894/1‑139 (MQ=255)                                                                                                        cTTTAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGa                           <  2:471821/139‑1 (MQ=255)                                                                                                         tttAAAATGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGaa                          <  2:210894/139‑1 (MQ=255)                                                                                                            aaaaTGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTg                                  >  2:185609/1‑128 (MQ=255)                                                                                                            aaaaTGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTg                                  <  1:185609/128‑1 (MQ=255)                                                                                                            aaaaTGAATAATATTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAAtct                       <  2:701489/139‑1 (MQ=255)                                                                                                                      atatTAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGt             <  1:1349256/139‑1 (MQ=255)                                                                                                                          tAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAcgcg                                                                                     >  1:283113/1‑63 (MQ=255)                                                                                                                          tAAGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAAcgcg                                                                                     <  2:283113/63‑1 (MQ=255)                                                                                                                            aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTgg                                                                       <  1:355112/75‑1 (MQ=255)                                                                                                                            aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTgg                                                                       >  2:355112/1‑75 (MQ=255)                                                                                                                            aGCCACTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTgcgc       >  2:974842/1‑139 (MQ=255)                                                                                                                                 cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCaaa                                               <  2:320007/94‑1 (MQ=255)                                                                                                                                 cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCaaa                                               >  1:320007/1‑94 (MQ=255)                                                                                                                                 cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCa                     <  1:194344/120‑1 (MQ=255)                                                                                                                                 cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCa                     >  2:194344/1‑120 (MQ=255)                                                                                                                                 cTTATTCATGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCt  >  1:1183552/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   ACCCGTGGCAGGATTTATCGAAATTGCATGAGTTGCCGGAGTAATGGTCATCGGGGTATCTCCTTTATGAGTCATGGTATGAAGATACGCAGATTTACTCTTGCTTTAAAATGAATAATATTAAGCCACTTATTCACGAATCGAGAATGCTATGGATCTGACCCAACTGGAGATGTTCAACGCGGTTGCCGAAGCTGGCAGCATAACCCAGGCTGCAGCAAAAGTGCATCGCGTGCCGTCGAATCTCACTACCCGTCTGCGCCAGCT  >  NC_000913/1614653‑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: ‑