Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,632,711	G→A	R454C (CGT→TGT)	guaB ←	inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,632,711	0	G	A	100.0%	78.0 / NA	26	R454C (CGT→TGT)	guaB	inosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

GTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACC  >  NC_000913/2632472‑2632957                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       gctgattccactgatacacgacgTCCAGCATGAGGGTCCGATGCTTATGAATGTTTTCCGTCATTGAAGCTTATTCCGAGGCACGTGAAACAGATAATATAAATAGCCCGAACTAAAGTCGGGCGGAGAGAATCAGGAGTCCAGACGGTCGATCGGGGACTCTTTAGTACTGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGATGATACGTACAAACTCCGCTTTAGTACACAGTtcgtcg                                                                                                                                                                                                                                              >  2:21841/23‑250 (MQ=255)        gtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTAc                                                                                                                                                                                                                                       <  1:230648/250‑1 (MQ=255)              agaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTAATCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGcc                                                                                                                                                                                                                                <  2:31825/251‑1 (MQ=255)                    gAAGTCCAGAATGAGGCAGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATAGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTc                                                                                                                                                                                                                            >  2:245218/1‑249 (MQ=255)                            gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATAc                                                                                                                                                                                                                  >  2:174967/1‑251 (MQ=255)                             aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACa                                                                                                                                                                                                                 <  1:60853/251‑1 (MQ=255)                                                         tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAAGCGCCCGACAGGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTg                                                                                                                                                                                      <  2:87544/250‑1 (MQ=255)                                                              cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGAt                                                                                                                                                                                <  1:174967/251‑1 (MQ=255)                                                                    cGCTTATTCCGAGGCAAGTGAAACAGATAATACAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCttt                                                                                                                                                                          >  1:212164/1‑251 (MQ=255)                                                                                  cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGTTAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTtat                                                                                                                                                            >  2:3083/1‑251 (MQ=255)                                                                                                  tataAATCGCCCGACATGAAATCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACTCGACCTTc                                                                                                                                            >  1:89370/1‑251 (MQ=255)                                                                                                      actCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGCAGTTCGGGGCCTCTTTAGTCATGTTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCCCCGCTGTTACGTACAAACTCCGCTTTAGTACACAGTCCGTCGATAGTACCACAGCCGGTCCGACCCATACAGGAGCGCAGGCCGCCCATCTGATGGAGAATGATCTCTTTCAGCCGACCTTTATGGGCACCGCGACCTTCGAt                                                                                                                                         >  2:245736/3‑250 (MQ=255)                                                                                                       aTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGGCCTTTATAGGCTACGCGACCTTCGATAc                                                                                                                                       >  1:186896/1‑251 (MQ=255)                                                                                                                                       aGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGTAAGTGGCTTTCCTGAATGCCCGCAACGCTGATACGTACAAACTCCGCTTTAGTACAAAGTTAGTCGCTCGTACCACAGCCGGTCCGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAAGCTACGCGACATTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTAt                                                                                                         >  1:166744/1‑249 (MQ=255)                                                                                                                                       aGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATTCCCGCACCGCTGATACGTACAAACTCCGCGTTAGTACACAGTTCGTCGATAGTACCCCAGCCGGTCAGACCCAGACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCGCTTTCAGGCGACCGTTATAGGCTACGCGACCTTCTATACCTTCGGGCTCCAGTTTGTCGTCAGGGGTATCg                                                                                                       <  2:265257/251‑1 (MQ=255)                                                                                                                                                      aGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGGCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACg                                                                                         <  1:19296/250‑1 (MQ=255)                                                                                                                                                                        tAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCa                                                                                                                >  2:145040/1‑209 (MQ=255)                                                                                                                                                                        tAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCa                                                                                                                <  1:145040/209‑1 (MQ=255)                                                                                                                                                                                  gtcgtGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGCTAGTACCACAGCCGGTCAGACCCCTACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTATGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGcccc                                                            >  1:129282/1‑250 (MQ=255)                                                                                                                                                                                    cgtGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAg                                                           <  1:232784/250‑1 (MQ=255)                                                                                                                                                                                                         ggCCCGCTCCGCTGATACCTACAAACTCCCCGTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCATGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTTTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGTAACCTTTGGACATATCGCCCAGGGAACCCAGACCTCGGTAAGAt                                     <  2:92742/250‑1 (MQ=255)                                                                                                                                                                                                                     tGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACgg                         >  1:27766/1‑251 (MQ=255)                                                                                                                                                                                                                          cGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTg                    <  2:212164/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGtt            <  1:245218/249‑1 (MQ=255)                                                                                                                                                                                                                                           aGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCTATTTCAc   <  1:74169/251‑1 (MQ=255)                                                                                                                                                                                                                                            gTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc  >  1:246190/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       GTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACC  >  NC_000913/2632472‑2632957

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