Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,792,331	(CG)4→3	coding (1059‑1060/1449 nt)	gabD →	NADP(+)‑dependent succinate‑semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,792,324	0	G	.	100.0%	141.8 / NA	34	coding (1052/1449 nt)	gabD	NADP(+)‑dependent succinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (15/19);  total (15/19)
*	NC_000913	2,792,325	0	C	.	100.0%	141.8 / NA	34	coding (1053/1449 nt)	gabD	NADP(+)‑dependent succinate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (15/19);  total (15/19)

TGCAGCAGGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTC  >  NC_000913/2792198‑2792457                                                                                                                               ||                                                                                                                                     tGCAGCAGGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGgcgcgcgtggtttg                                                                                                                          <  1:351447/140‑8 (MQ=255)        ggcagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGt                                                                                                                 >  1:59648/1‑140 (MQ=255)         gcagTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTa                                                                                                                <  2:953708/140‑1 (MQ=255)             tGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAgc                                                                                                            >  1:803617/1‑140 (MQ=255)                  aaaCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTa                                                                                                                <  2:346910/131‑1 (MQ=255)                  aaaCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTa                                                                                                                >  1:346910/1‑131 (MQ=255)                   aaCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGa                                                                                                      >  1:878216/1‑140 (MQ=255)                    aCTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGaa                                                                                                     >  1:128718/1‑140 (MQ=255)                        cacaTCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAAcgcg                                                                                                 <  2:290745/140‑1 (MQ=255)                            tCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGcggcgg                                                                                             <  2:128718/140‑1 (MQ=255)                                cGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCaa                                                                                          <  2:337708/139‑1 (MQ=255)                                cGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCaa                                                                                          >  1:337708/1‑139 (MQ=255)                                cGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAAc                                                                                         >  2:46538/1‑140 (MQ=255)                                        tGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCa                                                                                 <  2:525539/140‑1 (MQ=255)                                          gATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGC‑AAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGcc                                                                              <  2:383273/140‑1 (MQ=255)                                             aaCGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGa                                                                            <  2:878216/140‑1 (MQ=255)                                                ggCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCa                                                                         <  2:979200/140‑1 (MQ=255)                                                  cGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCAtt                                                                       >  1:430200/1‑140 (MQ=255)                                                  cGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCAtt                                                                       >  1:350809/1‑140 (MQ=255)                                                     cACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGcc                                                                              >  1:236788/1‑130 (MQ=255)                                                     cACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGcc                                                                              <  2:236788/130‑1 (MQ=255)                                                                  cTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCgg                                                       >  1:471642/1‑140 (MQ=255)                                                                   tGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGc                                                      <  1:46538/140‑1 (MQ=255)                                                                   tGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGc                                                      <  1:276022/140‑1 (MQ=255)                                                                           gAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCa                                              <  2:59648/140‑1 (MQ=255)                                                                                 gCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAgtgt                                        >  2:106007/1‑140 (MQ=255)                                                                                    gTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGa                                     <  2:350809/140‑1 (MQ=255)                                                                                    gTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGa                                     <  2:430200/140‑1 (MQ=255)                                                                                      aGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGaaa                                   <  2:865851/140‑1 (MQ=255)                                                                                          aaaGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAaagaag                               <  2:803617/140‑1 (MQ=255)                                                                                                    agCATATTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCgg                     >  2:374151/1‑140 (MQ=255)                                                                                                       atatTGCCGATGCGCTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGccc                  <  2:688839/140‑1 (MQ=255)                                                                                                                     cTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGt    >  2:241636/1‑140 (MQ=255)                                                                                                                     cTGGAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGt    >  2:417682/1‑140 (MQ=255)                                                                                                                       ggAGAAAG‑‑GCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTc  <  1:374151/140‑1 (MQ=255)                                                                                                                               ||                                                                                                                                     TGCAGCAGGCAGTGAGCAAACTGCACATCGGCGACGGGCTGGATAACGGCGTCACCATCGGGCCGCTGATCGATGAAAAAGCGGTAGCAAAAGTGGAAGAGCATATTGCCGATGCGCTGGAGAAAGGCGCGCGCGTGGTTTGCGGCGGTAAAGCGCACGAACGCGGCGGCAACTTCTTCCAGCCGACCATTCTGGTGGACGTTCCGGCCAACGCCAAAGTGTCGAAAGAAGAGACGTTCGGCCCCCTCGCCCCGCTGTTC  >  NC_000913/2792198‑2792457

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