Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,424,386	G→A	*229* (TAG→TAA)	ulaF →	L‑ribulose 5‑phosphate 4‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,424,386	0	G	A	100.0%	142.9 / NA	45	*229* (TAG→TAA)	ulaF	L‑ribulose 5‑phosphate 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (26/19);  total (26/19)

GATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCT  >  NC_000913/4424254‑4424517                                                                                                                                     |                                                                                                                                    gATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                               <  2:876744/139‑1 (MQ=255)   tGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcgcg                                                                                                                             <  1:325153/139‑1 (MQ=255)         aagTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCgg                                                                                                                       >  2:699327/1‑139 (MQ=255)            tGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                               <  2:150930/128‑1 (MQ=255)            tGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACAcg                                                                                                                               >  1:150930/1‑128 (MQ=255)                    aTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATTGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGaaa                                                                                                                    >  2:786372/1‑131 (MQ=255)                    aTGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGaaa                                                                                                                    <  1:786372/131‑1 (MQ=255)                      ggCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTtccttcct                                                                                                          <  1:519088/139‑1 (MQ=255)                            gATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGgag                                                                                                    <  1:179587/139‑1 (MQ=255)                                 cccGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAAc                                                                                               >  1:295513/1‑139 (MQ=255)                                                     cTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGt                                                                           >  2:153399/1‑139 (MQ=255)                                                        aaTCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                     <  2:819978/110‑1 (MQ=255)                                                        aaTCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                     >  1:819978/1‑110 (MQ=255)                                                           cacaTCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCa                                                                     <  1:153399/139‑1 (MQ=255)                                                                  aCAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaa                                                              <  1:976837/139‑1 (MQ=255)                                                                        tCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGccccc                                                        >  1:475877/1‑139 (MQ=255)                                                                                       acTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCACCTGCGAAGGGGgcaaa                                         >  1:360477/1‑139 (MQ=255)                                                                                        cTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGGAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCTTGCGAAGGGGgcaaag                                        <  2:681116/139‑1 (MQ=255)                                                                                         ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagc                                       <  2:293572/139‑1 (MQ=255)                                                                                         ttCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagc                                       >  2:1106345/1‑139 (MQ=255)                                                                                          tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                      >  2:843333/1‑139 (MQ=255)                                                                                          tCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaagca                                      >  1:31077/1‑139 (MQ=255)                                                                                           cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAAttt                                                                                                                 >  1:782266/1‑63 (MQ=255)                                                                                           cATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAAttt                                                                                                                 <  2:782266/63‑1 (MQ=255)                                                                                             tGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                   >  2:1110723/1‑139 (MQ=255)                                                                                                   aaCACGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACT‑GTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCa                            >  2:233008/1‑139 (MQ=255)                                                                                                    acacGGTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCa                            <  1:866548/139‑1 (MQ=255)                                                                                                        ggTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGt                        >  2:524689/1‑139 (MQ=255)                                                                                                         gTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                       >  2:889420/1‑139 (MQ=255)                                                                                                         gTCCTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                       >  2:1084120/1‑139 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                    <  2:295513/139‑1 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                    >  1:665498/1‑139 (MQ=255)                                                                                                            cTAACGCTTATTACGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                    <  1:1110723/139‑1 (MQ=255)                                                                                                            cTAACGCTTATTACCGGCAGTAGTAAAACACGCGCTGCGGAATTTTCCTTCCTCGGGAGATAACTAGTCTAATTCCGCATCCGTTTTTCAAAAAAAAGCCCCCTGCTAAGCGTTCAAAGCAAACTATTGCAATGTTTCg                    <  2:360477/139‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGc                                                                                 <  1:774320/70‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGc                                                                                 >  2:774320/1‑70 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                        <  2:389321/111‑1 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                        >  1:389321/1‑111 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat            >  2:12063/1‑139 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat            >  1:162832/1‑139 (MQ=255)                                                                                                                    tattaCGGGCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat            >  1:527347/1‑139 (MQ=255)                                                                                                                           ggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt     >  1:1053064/1‑139 (MQ=255)                                                                                                                           ggCAGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt     >  1:536510/1‑139 (MQ=255)                                                                                                                              aGAAGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTgtt  <  1:1084120/139‑3 (MQ=255)                                                                                                                                aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                     >  2:734163/1‑38 (MQ=255)                                                                                                                                aaGTAAAACACGCGCTGCGGAAATTTCCTTCCTCGGga                                                                                                     <  1:734163/38‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    GATGGAAGAAGTGGCGAAAATGGCGTGGATTGCCCGCGGCATTAACCCACAACTCAATCACATCGACAGCTTCCTGATGAATAAACACTTCATGCGTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCT  >  NC_000913/4424254‑4424517

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