Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,451,427	A→G	intergenic (‑78/+170)	tynA ← / ← paaZ	copper‑containing amine oxidase/oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,451,427	0	A	G	100.0%	103.4 / NA	30	intergenic (‑78/+170)	tynA/paaZ	copper‑containing amine oxidase/oxepin‑CoA hydrolase/3‑oxo‑5,6‑dehydrosuberyl‑CoA semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/17);  total (13/17)

CTGCCAACGCCAGGGTTGTTTTACGGGCAGAATACAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTG  >  NC_000913/1451298‑1451554                                                                                                                                  |                                                                                                                                cTGCCAACGCCAGGGTTGTTTTACGGGCAGAATACAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTg                                                                                                                              <  1:415874/133‑1 (MQ=255)    ccAACGCCAGGGTTGTTTTACGGGCAGAATACAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATg                                                                                                                           >  1:352717/1‑133 (MQ=255)       aCGCCAGGGTTGTTTTACGGGCAGAATACAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCaa                                                                                                                        <  1:39056/133‑1 (MQ=255)       aCGCCAGGATTGTTTCACTGGCCGAATACAAAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCaa                                                                                                                        <  1:239564/133‑1 (MQ=255)            aGGGTTGTATTACGGGCATAATCCAGAGAGCGGTTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAgggg                                                                                                                   <  1:220021/133‑1 (MQ=255)             gggTTGTTTTACGGGTAGCATATAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGa                                                                                                                  <  1:163690/133‑1 (MQ=255)                                     gagagGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATCGCCAGCGCCTGCCAGGTTGttt                                                                                          >  1:208974/1‑133 (MQ=255)                                         ggggCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGCGATTGCCGGCAGCTGCCAGGTGGTTTGGCa                                                                                      >  1:137439/1‑133 (MQ=255)                                           ggCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGgttgt                                                                                            <  1:82541/125‑1 (MQ=255)                                           ggCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGgttgt                                                                                            >  2:82541/1‑125 (MQ=255)                                                  cATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCa                                                                             <  1:290673/133‑1 (MQ=255)                                                  cATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCa                                                                             <  2:352717/133‑1 (MQ=255)                                                    ttattaACCTCGTCATATGTTGAGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGc                                                                           <  1:213995/133‑1 (MQ=255)                                                    ttattaACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCa                                                                             >  1:31696/1‑131 (MQ=255)                                                    ttattaACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCa                                                                             <  2:31696/131‑1 (MQ=255)                                                             tCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATAc                                                                  <  2:35067/133‑1 (MQ=255)                                                              cGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGGTGTTTGGGAGGGGTACCAGCTTTTCATACa                                                                 >  1:147910/1‑133 (MQ=255)                                                               gTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAg                                                                >  1:397340/1‑133 (MQ=255)                                                               gTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAg                                                                >  2:360865/1‑133 (MQ=255)                                                                       ttgtGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGcc                                                        <  2:682/133‑1 (MQ=255)                                                                       ttgtGTTCTTCTTAGGAACCGCGCCCTGTGGGTGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGcc                                                        <  2:137439/133‑1 (MQ=255)                                                                            ttCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGaa                                                   >  1:34834/1‑133 (MQ=255)                                                                                       aaCCGCGCTCCGTGGCCGGTTTAGCCAGCTTCTCATTATCATTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTac                                        <  2:147910/133‑1 (MQ=255)                                                                                            cgcTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACGGTGGATGGCCTGAACATAGATATACAcatc                                   >  1:143961/1‑133 (MQ=255)                                                                                              ctctGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGt                                          >  1:91880/1‑124 (MQ=255)                                                                                              ctctGTGGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGt                                          <  2:91880/124‑1 (MQ=255)                                                                                                 tgtgGGCGGTTTAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGTACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCa                              <  2:301464/133‑1 (MQ=255)                                                                                                 tgtgGGCGGTATAGTTAGGTTCACATTAACAGTGCTGATGCCAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCa                              <  2:143961/133‑1 (MQ=255)                                                                                                           ttAGTCAGGTTCACATTATCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGTATGCCCTGCAAAAAGAAGTACACATCATGCATAATGTGACa                    >  2:185843/1‑133 (MQ=255)                                                                                                                             tCAGTGCTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTg  >  2:79295/1‑133 (MQ=255)                                                                                                                                  |                                                                                                                                CTGCCAACGCCAGGGTTGTTTTACGGGCAGAATACAGAGAGGGGCTTCCCATTATTAACCTCGTCAGATGTTGTGTTCTTGTTAGCAACCGCGCTCTGTGGGCGGTTTAGTCAGGTTCACATTATCAGTACTGATGCAAAGGGGATTGCCTGCACCTGCCAGGTTGTTTGGCAGGTGTGCCAGCTTTTCATACAGTGGATGCCCTGAAAATAGATGTACACATCATGCATAATGTGACAACGTCACAAAACTTAGTG  >  NC_000913/1451298‑1451554

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