Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,573,508	T→C	intergenic (‑6/+267)	glgB ← / ← asd	1,4‑alpha‑glucan branching enzyme/aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,573,508	0	T	C	100.0%	112.2 / NA	31	intergenic (‑6/+267)	glgB/asd	1,4‑alpha‑glucan branching enzyme/aspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/17);  total (14/17)

TAGCGTCGGGTAAAAGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGTCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAGGCTGTGAA  >  NC_000913/3573373‑3573638                                                                                                                                        |                                                                                                                                   tAGCGTCGGGTAAAAGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTc                                                                                                                                <  2:353202/140‑1 (MQ=255)         ggTAAAAGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTa                                                                                                                        >  1:129232/1‑140 (MQ=255)         ggTAAAAGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTa                                                                                                                        <  1:214518/140‑1 (MQ=255)              aaGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATtgtg                                                                                                                   <  1:281268/140‑1 (MQ=255)                  gCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCAc                                                                                                               <  2:117791/140‑1 (MQ=255)                   cACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTc                                                                                                                                <  1:199268/122‑1 (MQ=255)                   cACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTc                                                                                                                                >  2:199268/1‑122 (MQ=255)                   cACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCAcc                                                                                                              >  1:163260/1‑140 (MQ=255)                                    gcgGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGAtttt                                                                                             >  2:382718/1‑140 (MQ=255)                                                       ccAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACAt                                                                          >  2:12569/1‑140 (MQ=255)                                                          gTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACAtttt                                                                       >  1:304606/1‑140 (MQ=255)                                                           tACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACAttttt                                                                      <  2:163260/140‑1 (MQ=255)                                                             cGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACAttttttt                                                                    >  2:85387/1‑140 (MQ=255)                                                               gAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGa                                                                  >  2:269431/1‑140 (MQ=255)                                                                     ggATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGaa                                                            <  1:85387/140‑1 (MQ=255)                                                                     ggATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGaa                                                            >  2:5595/1‑140 (MQ=255)                                                                     ggATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGaa                                                            <  1:302799/140‑1 (MQ=255)                                                                         ccGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGgc                                                        >  1:411813/1‑140 (MQ=255)                                                                                aTGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCAc                                                 <  2:304606/140‑1 (MQ=255)                                                                                 tGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCAc                                                 <  1:336203/139‑1 (MQ=255)                                                                                 tGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCAc                                                 >  2:336203/1‑139 (MQ=255)                                                                                  ggCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCActc                                               <  2:129232/140‑1 (MQ=255)                                                                                    ccTGCAATTAGCGCGTTAATCACGTCTCTATCCATACGATTGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGATGTGACATTTTTTTGACCTGAACGGCGCAGCACTCtg                                             <  2:370554/140‑1 (MQ=255)                                                                                      tGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCtgtg                                           <  1:269431/140‑1 (MQ=255)                                                                                        cAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCa                                         <  1:5595/140‑1 (MQ=255)                                                                                                                tATCGATACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTc                 >  1:145327/1‑140 (MQ=255)                                                                                                                     aTACGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTa            <  2:396144/140‑1 (MQ=255)                                                                                                                       aCGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAg          >  2:135455/1‑140 (MQ=255)                                                                                                                       aCGATCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAg          <  1:382718/140‑1 (MQ=255)                                                                                                                          aTCGGACATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAGGCt       >  1:127199/1‑140 (MQ=255)                                                                                                                               aCATGCTTGCCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAGGCTGTGaa  <  2:210644/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   TAGCGTCGGGTAAAAGGGCACGGACTTCCAGTCCCGCGGTGGTTTTATGCATTCCCAGTACGGAAAAAGGATCCGCAAAATGGCCTGCAATTAGCGCGTTAATCACGTCTCTATCGATACGATCGGACATGCTTGTCTTCCTGTTTTATTGTGTCACCCCATCCAATTGGATTTTTGACTTCTGGTTGTGACATTTTTTTGACCTGAACGGCGCAGCACTCTGTGCATCCTCTCTGCGTCGTCCTCACTTCAGGTAAGGCTGTGAA  >  NC_000913/3573373‑3573638

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