Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,490,141	C→T	*333* (TAG→TAA)	idnR ←	transcriptional repressor, 5‑gluconate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,490,141	0	C	T	100.0%	130.3 / NA	41	*333* (TAG→TAA)	idnR	transcriptional repressor, 5‑gluconate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/18);  total (23/18)

GATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGC  >  NC_000913/4490066‑4490275                                                                            |                                                                                                                                       ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGAt                                                                                                                     <  1:32910‑M2/34‑1 (MQ=255) ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGAt                                                                                                                     >  2:32910‑M2/62‑95 (MQ=255) ttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTa                                                                         >  1:167407‑M2/62‑139 (MQ=255)   cactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggttg                                                                                  >  2:671469‑M2/60‑128 (MQ=255)   cactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTgttgtggttg                                                                                  <  1:671469‑M2/69‑1 (MQ=255)     ctcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGGGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtttt                                                                     <  2:880923‑M2/82‑1 (MQ=255)      tcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAccc                                                                                                        <  2:519223‑M2/47‑1 (MQ=255)      tcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATAccc                                                                                                        >  1:519223‑M2/57‑103 (MQ=255)                tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTaaa                                                                                                                                    >  2:451639‑M2/47‑65 (MQ=255)                tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTaaa                                                                                                                                    <  1:451639‑M2/19‑1 (MQ=255)                       gtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgt                                                                                            >  2:943717‑M2/40‑98 (MQ=255)                       gtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAgt                                                                                            <  1:943717‑M2/59‑1 (MQ=255)                       gtttcctacctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATa                                                                                                                    <  1:396396‑M2/35‑1 (MQ=255)                       gtttcctacctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATa                                                                                                                    >  2:396396‑M2/40‑74 (MQ=255)                        tttcctgcctggagtatagtcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                             <  1:721377‑M2/58‑1 (MQ=255)                        tttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                             >  2:721377‑M2/39‑96 (MQ=255)                          tcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACAttt                                                >  2:556163‑M2/37‑139 (MQ=255)                           cctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTg                                               >  1:698441‑M2/36‑139 (MQ=255)                               cctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                                           <  1:556163‑M2/108‑1 (MQ=255)                                 tggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                        <  1:392791‑M2/95‑1 (MQ=255)                                 tggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                        >  2:392791‑M2/30‑124 (MQ=255)                                    agtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                                      >  2:717571‑M2/27‑139 (MQ=255)                                        taggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGa                                  >  2:101769‑M2/23‑139 (MQ=255)                                             aacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAat                             >  2:266243‑M2/18‑139 (MQ=255)                                              acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                             >  2:748648‑M2/17‑74 (MQ=255)                                              acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                             <  1:748648‑M2/58‑1 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                                          <  2:343329‑M2/77‑1 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGtt                                                                          >  1:343329‑M2/16‑92 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                                <  2:171859‑M2/119‑1 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                                >  1:171859‑M2/16‑134 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat                           >  1:52628‑M2/16‑139 (MQ=255)                                               cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat                           >  2:304928‑M2/16‑139 (MQ=255)                                                gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                          >  2:720192‑M2/15‑139 (MQ=255)                                                gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                          >  1:110611‑M2/15‑139 (MQ=255)                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGtgttg                                                                                        >  2:383673‑M2/10‑72 (MQ=255)                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGtgttg                                                                                        <  1:383673‑M2/63‑1 (MQ=255)                                                     gcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAAc                     <  2:110611‑M2/130‑1 (MQ=255)                                                         agagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACgggg                 <  1:717571‑M2/134‑1 (MQ=255)                                                                     aaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCt     >  2:420540/1‑139 (MQ=255)                                                                        tAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGc  <  1:304928/139‑1 (MQ=255)                                                                        tAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGc  <  1:420540/139‑1 (MQ=255)                                                                            |                                                                                                                                       GATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTACGCTGGC  >  NC_000913/4490066‑4490275

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