Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,667,460	G→A	Q369* (CAG→TAG)	mlc ←	DNA‑binding transcriptional repressor Mlc

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,667,460	0	G	A	100.0%	83.5 / NA	24	Q369* (CAG→TAG)	mlc	DNA‑binding transcriptional repressor Mlc
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

ACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTT  >  NC_000913/1667230‑1667687                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    aCCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCct                                                                                                                                                                                                                 <  2:75829/251‑1 (MQ=255)        tGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACg                                                                                                                                                                                                          >  2:25314/1‑251 (MQ=255)                       aaCCAGGAAAATAATTAACCTTGAAATGCGAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAAGTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGTGCCTGCTGACGGATGCTGTCTGAGAt                                                                                                                                                                                           <  2:205675/251‑1 (MQ=255)                        aCCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATg                                                                                                                                                                                          <  1:58731/251‑1 (MQ=255)                           aGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGAcc                                                                                                                                                                                       >  2:49932/1‑251 (MQ=255)                               aaaTAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCggg                                                                                                                                                                                    >  1:97092/1‑250 (MQ=255)                                aataatTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGa                                                                                                                                                                                   >  1:18713/1‑250 (MQ=255)                                           ttGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTg                                                                                                                                                                       >  1:268202/1‑251 (MQ=255)                                             gAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACTAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGc                                                                                                                                                                      <  1:140096/250‑1 (MQ=255)                                                                    ggCCCAAAGTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAAt                                                                                                                                              <  2:312172/251‑1 (MQ=255)                                                                                         aaTTTTGGTAGCACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaa                                                                                                                          >  1:251433/1‑250 (MQ=255)                                                                                         aaTTTTGGTAGAACAGTTAAAAAATGTTAACCCTGAAACAGACGAATCCACACAGAACCGTTATACATCACGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGACAACTGAGTACTCTCAACGCTGATGTGATAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTtaa                                                                                                                          >  2:170039/1‑250 (MQ=255)                                                                                                           aaaaaaTGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATgg                                                                                                        <  1:293394/250‑1 (MQ=255)                                                                                                                                                aaCCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGACGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATATGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAatg                                                                  >  1:189753/1‑251 (MQ=255)                                                                                                                                                       ataCATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCtt                                                           >  2:306752/1‑251 (MQ=255)                                                                                                                                                               cgTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGt                                                   >  1:242787/1‑251 (MQ=255)                                                                                                                                                                    tttACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGa                                                <  1:295819/249‑1 (MQ=255)                                                                                                                                                                     ttACCAGTGCCGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTAATCTAACCGCGGATGTGCTAACTATACGCAGCAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATATGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCAACATGATGGCAAGAATGCGACAGACATGCGCGCCAACCCCGGTAATGATGTCTTTTGCCCGTAGAt                                               >  2:266683/1‑249 (MQ=255)                                                                                                                                                                                  gcCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccg                                <  1:120714/251‑1 (MQ=255)                                                                                                                                                                                  gcCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATgccg                                >  1:239942/1‑251 (MQ=255)                                                                                                                                                                                                      tAGAAAACTGAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCAc             <  1:46666/250‑1 (MQ=255)                                                                                                                                                                                                               gAGTACTCTCAACGCTGATGTTCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGGCTGAGATGACCGGGAAGAGGATAGCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAAGGAGTCCACGGTTAACg     <  2:200129/249‑1 (MQ=255)                                                                                                                                                                                                               gAGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGt   >  1:136163/1‑251 (MQ=255)                                                                                                                                                                                                                aGTACTCTCAACGCTGATGTGCTAACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGtt  <  2:97092/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    ACCGTTATGCTTCACAAATATAAACCAGGAAAATAATTAACCTTGAAAGTCTAAGTTATGCTTTCCTGGCCCAAATTGAGATAGCGCAAATTTTGGTAGAACAGTTAAAAAATGTTAACCCTGCAACAGACGAATCAACAAAGAACCGTTATACATCGCGTCTTTTACCAGTGCAGCGCCTGCCATCGTGCCCTGGTTAGAAAACTGAGTACTCTCAACGCTGATGTGCTGACTATACGCAGGAAGGGCCTGCTGACGGATGCTGTCTGAGATGACCGGGAAGAGGATATCTGCCGCTTTACTTAACGGTGAGCCAATCAGTATTTTTTGTGGGTTAAATAAATTCACCATGATGGCAAGAATGCGCCCGACATGCGCGCCCACCCCGGTAATGATGTCTTTTGCCAGTAGATCGCCGCGCAATGCCGCCTGACACAATGAGTCCACGGTTAACGGTT  >  NC_000913/1667230‑1667687

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