Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,620,766	C→T	G388G (GGC→GGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,766	0	C	T	100.0%	85.3 / NA	28	G388G (GGC→GGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/16);  total (12/16)

CTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAA  >  NC_000913/4620526‑4621007                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gacggatggctggccgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTgg                                                                                                                                                                                                                                         >  2:233258‑M2/224‑251 (MQ=255)   cggatggctggccgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCa                                                                                                                                                                                                                                       <  1:249767‑M2/30‑1 (MQ=255)                gctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTa                                                                                                                                                                                                                          <  2:163401‑M2/43‑1 (MQ=255)                   gtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTAc                                                                                                                                                                                                                         <  1:143203‑M2/44‑1 (MQ=255)                   gtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATAGTTTGGTAc                                                                                                                                                                                                                         <  2:313684‑M2/44‑1 (MQ=255)                         atcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGat                                                                                                                                                                                                                  >  2:415458‑M2/200‑250 (MQ=255)                            ccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATgg                                                                                                                                                                                                              <  1:284963‑M2/55‑1 (MQ=255)                            ccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaaaactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATgg                                                                                                                                                                                                              >  2:204477‑M2/197‑251 (MQ=255)                                  gaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGCCTCTGGCAAACTATTTTGGTACTTCTGATATGGAATATg                                                                                                                                                                                                        >  1:464069‑M2/191‑251 (MQ=255)                                  gaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATg                                                                                                                                                                                                        <  2:511372‑M2/61‑1 (MQ=255)                                         aagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGc                                                                                                                                                                                                 >  1:412230‑M2/184‑251 (MQ=255)                                               cacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGtt                                                                                                                                                                                           <  2:483774‑M2/74‑1 (MQ=255)                                                   taatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacacgatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTg                                                                                                                                                                                        >  1:591080‑M2/174‑250 (MQ=255)                                                     atcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgatcataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAATCCATGTTCTGata                                                                                                                                                                                     <  1:193515‑M2/80‑2 (MQ=255)                                                     atcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGAt                                                                                                                                                                                      >  2:377701‑M2/172‑250 (MQ=255)                                                       cagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCtgataat                                                                                                                                                                                   <  1:62792‑M2/82‑4 (MQ=255)                                                                  atacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAaagccatgttctgataatgggggcggctt                                                                                                                                                                        >  1:105492‑M2/159‑237 (MQ=255)                                                                  atacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAagccatgttctgataatgggtgcggct                                                                                                                                                                         <  2:167743‑M2/92‑14 (MQ=255)                                                                           gaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGaatatggcaaagccatgttctgataatgggtgcggcttgatgccctc                                                                                                                                                               >  2:57547‑M2/150‑228 (MQ=255)                                                                                                                                                                                                                agttggagtcattaccaccttcgcggatattggTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcg                          <  2:158780‑M1/235‑157 (MQ=255)                                                                                                                                                                                                                    ggagtcattaccaccttcgcggataTTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggggggggtgaaaacctgttgataaggataaaacaatggctacaccacacattaatgcaggaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagt                      >  2:601573‑M1/13‑91 (MQ=255)                                                                                                                                                                                                                      agtcattaccaccttcgcggaTATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccccctaccccagggggagggggaaaacatgttgataaggataaaacaaaggctaccccacacattaatgcagaaatgggcggtttcgctgacgtaggttttattccaagcgacccgctgcgtgcgaagtat                    >  1:567934‑M1/11‑89 (MQ=255)                                                                                                                                                                                                                         cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccccagggggagggggaaaacaagttgataaggataaaacaatggctaccccacccattaatgcagaaaagggcgatttcgctgacgtagttttgctgccaggcgaccccctgcgggcgaaggatatt                 >  1:183732‑M1/8‑86 (MQ=255)                                                                                                                                                                                                                          attaccaccttcgCGGATAGTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctgccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattg                <  2:314312‑M1/245‑167 (MQ=255)                                                                                                                                                                                                                              ccaccTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctga            <  1:462901‑M1/249‑171 (MQ=255)                                                                                                                                                                                                                              ccaccTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctga            <  2:591464‑M1/249‑171 (MQ=255)                                                                                                                                                                                                                                     cgcgGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATAGGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacgatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtaggtttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaacttt      <  2:432537‑M1/250‑177 (MQ=255)                                                                                                                                                                                                                                        ggATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttcctt  <  2:251231‑M1/251‑181 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  CTGGGGCCACCGTCGCGACGTCGCCGGTTATGCCGCGGGTCTGGAACTGTTCGACCGCCGTCTGCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAA  >  NC_000913/4620526‑4621007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑