Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,981	G→C	T23S (ACC→AGC)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,981	0	G	C	100.0%	87.8 / NA	29	T23S (ACC→AGC)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (9/20);  total (9/20)

CAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACA  >  NC_000913/69743‑70212                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                        atgcctgttcccggcGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAATCCa                                                                                                                                                                                                                             <  2:160251/237‑1 (MQ=255)                  cATCAATCGGGGAGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTCTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCTCAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGcaca                                                                                                                                                                                                            <  2:192789/251‑1 (MQ=255)                                   gTCGAGCTGGTACTGTCAACGCCAATCCCGACCGCAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACAGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATc                                                                                                                                                                                           <  2:259799/251‑1 (MQ=255)                                   gTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATc                                                                                                                                                                                           >  2:143514/1‑251 (MQ=255)                                        gCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGc                                                                                                                                                                                      <  1:375419/251‑1 (MQ=255)                                          cGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCa                                                                                                                                                                                    <  1:89625/251‑1 (MQ=255)                                                    aaCGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCg                                                                                                                                                                          >  1:241955/1‑251 (MQ=255)                                                      cGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGcc                                                                                                                                                                        <  1:315005/251‑1 (MQ=255)                                                             cccGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGttt                                                                                                                                                                 <  2:210124/251‑1 (MQ=255)                                                                                                  gCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACGGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAg                                                                                                                               <  2:94501/248‑1 (MQ=255)                                                                                                     aGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTag                                                                                                                            >  1:267690/1‑248 (MQ=255)                                                                                                        cgggTTTTCAGTGCCGCTTCCATTGACTCAATGTAGGCACGCTGATGATGACGGAGCTTGTTATGCGGAGCATCACAAACTTGCCTTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGTCGCGGGCCACCGCCACAGCTCGCGCAGAATCACTGCCAAAATCGAGTGCAAGTACAATCGACATCGTTTCACTCCATTCAAAAAAACGGTTATGGAGAATCAGTatatatta                                                                                                                      <  2:121768/249‑8 (MQ=255)                                                                                                                    gCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTGTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCTTTTCACTCCATCCAAAAAAACGGTTATGGAGAAACAGGAGATAGTTGCGATAAAAAGc                                                                                                          <  2:45942/251‑1 (MQ=255)                                                                                                                               ttaCTCAATGTATTCACGCGGATGATGACGGAACTGTTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGTATACCACTCTCCGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACGCCATCCAAAAAAACGGGTATGGAGAAACAGTCGAGAGTTTCGATAAAAAGCGTCAGGTAgg                                                                                                <  2:192881/248‑1 (MQ=255)                                                                                                                                                  ggATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTc                                                                                                                                                                <  2:142207/167‑1 (MQ=255)                                                                                                                                                  ggATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTc                                                                                                                                                                >  1:142207/1‑167 (MQ=255)                                                                                                                                                      gatgaCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATaaaaa                                                                        <  1:238148/251‑1 (MQ=255)                                                                                                                                                         gaCGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGACACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCACTTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGCAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGc                                                                     >  1:183673/1‑251 (MQ=255)                                                                                                                                                                    ttATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGc                                                          >  2:313489/1‑251 (MQ=255)                                                                                                                                                                          ggggCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAgtg                                                    <  2:164719/251‑1 (MQ=255)                                                                                                                                                                                      aaTTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCa                                                                                                                                                                                                                             <  1:313826/70‑1 (MQ=255)                                                                                                                                                                                      aaTTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCa                                                                                                                                                                                                                             >  2:313826/1‑70 (MQ=255)                                                                                                                                                                                                 ttGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAAt                              <  1:49456/249‑1 (MQ=255)                                                                                                                                                                                                   gCCAACGGGGATACCACTCGACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAAtgt                            <  2:112739/250‑1 (MQ=255)                                                                                                                                                                                                    ccAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTgg                          >  1:294770/1‑251 (MQ=255)                                                                                                                                                                                                                 cACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGc                                                          <  1:81295/206‑1 (MQ=255)                                                                                                                                                                                                                 cACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGc                                                          >  2:81295/1‑206 (MQ=255)                                                                                                                                                                                                                       aCGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATta        <  1:313489/250‑1 (MQ=255)                                                                                                                                                                                                                            ggtggCGATCTCTTCACCGCTGGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGaca  <  1:193577/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                        CAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCGACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCCACCGCCAAAGCTCGCACAGAATCACTGCCAAAATCGAGGCCAATTGCAATCGCCATCGTTTCACTCCATCCAAAAAAACGGGTATGGAGAAACAGTAGAGAGTTGCGATAAAAAGCGTCAGGTAGGATCCGCTAATCTTATGGATAAAAATGCTATGGCATAGCAAAGTGTGACGCCGTGCAAATAATCAATGTGGACTTTTCTGCCGTGATTATAGACA  >  NC_000913/69743‑70212

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