Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	354,768	A→G	intergenic (+176/‑154)	prpE → / → codB	propionyl‑CoA synthetase/cytosine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,768	0	A	G	100.0%	89.2 / NA	26	intergenic (+176/‑154)	prpE/codB	propionyl‑CoA synthetase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/12);  total (14/12)

TGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTT  >  NC_000913/354631‑354897                                                                                                                                          |                                                                                                                                  tgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGcc                                                                                                                                 >  1:23290/1‑140 (MQ=255)    gCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCAc                                                                                                                              <  2:828557/140‑1 (MQ=255)    gCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCAc                                                                                                                              <  1:827523/140‑1 (MQ=255)           ccTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGTTGCCTGATACGAAGCTAACGCGTCTAATCAGGCCTGCCCACTGttttt                                                                                                                       >  2:346490/1‑140 (MQ=38)           ccTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGttttt                                                                                                                       >  2:529639/1‑140 (MQ=255)            cTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTa                                                                                                                      <  2:557557/140‑1 (MQ=255)              gATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTaca                                                                                                                    <  2:23290/140‑1 (MQ=255)                      gCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAAt                                                                                                            <  1:78212/140‑1 (MQ=255)                                          cTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCa                                                                                        <  1:591679/140‑1 (MQ=255)                                          cTACGGACCTGAACCATAGGTCGGAGAAGACGCTCGCCTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCa                                                                                        <  1:339065/140‑1 (MQ=39)                                            aCGGACCTGAACCGTAGGACGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGACTGCCCACCGTTTTTACACCGATAACTTTTTCCCCACCTTTTTGCACt                                                                                      >  2:339065/1‑140 (MQ=37)                                              ggACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACtca                                                                                    <  1:20846/140‑1 (MQ=255)                                                            ggTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAat                                                                      >  2:19708/1‑140 (MQ=255)                                                               cGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAatata                                                                   >  1:721410/1‑140 (MQ=255)                                                                            cgcgTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGaaa                                                      >  1:733870/1‑140 (MQ=255)                                                                                tCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCACTAATATAAAAAATATTGTTCCCCAAGAAAACGa                                                  >  2:421709/1‑140 (MQ=39)                                                                                 cGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGAt                                                 <  2:733870/140‑1 (MQ=255)                                                                                 cGCATCCGACACCATGCTCAGACGCATGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGAt                                                 <  1:365553/140‑1 (MQ=255)                                                                                  gCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGAtt                                                <  1:19708/140‑1 (MQ=255)                                                                                              aTGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCtt                                    >  1:431607/1‑140 (MQ=255)                                                                                                gCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCa                                  >  1:763868/1‑140 (MQ=255)                                                                                                     gATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGa                             >  2:621286/1‑140 (MQ=255)                                                                                                         ccTGATGCGACGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAg                         >  1:28010/1‑140 (MQ=255)                                                                                                                   cGCTGACGCGTCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGcggcgg               >  2:167338/1‑140 (MQ=255)                                                                                                                             tCTTATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTggg     <  2:763868/140‑1 (MQ=255)                                                                                                                                tATCAGGCCTGCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGttt  >  2:820040/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                  TGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTACGGACCTGAACCGTAGGTCGGATAAGGCGCTCGCGTCGCATCCGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTT  >  NC_000913/354631‑354897

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