Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,419	C→T	Q70Q (CAG→CAA)	iclR ←	DNA‑binding transcriptional repressor IclR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,419	0	C	T	100.0%	105.5 / NA	33	Q70Q (CAG→CAA)	iclR	DNA‑binding transcriptional repressor IclR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/17);  total (16/17)

GCGCGGACATTCGCATCAGATGCGTA‑CACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA  >  NC_000913/4223184‑4223660                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    gcgcGGACATTCGCATCAGATGCGTA‑CACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTg                                                                                                                                                                                                                                      >  2:79515/1‑251 (MQ=255)      gactttCGCATAAGATGC‑TATCACTGTACTGGGTCGAAAAAAGTCTCTGCGTGATCGTTTTGAGCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAg                                                                                                                                                                                                                                 <  2:49602/247‑1 (MQ=255)                    aggCGTA‑CACTGTCCCGGGTCAAAAATAATCGCTTCGTGATCGTTTTGAACAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCg                                                                                                                                                                                                                   <  2:151408/249‑1 (MQ=255)                                                aatCGCTTCGTGATCGCTTTGATTAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCCTTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGTCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGTGTTGGGTAACCCGGCTTGTTgcg                                                                                                                                                                                        >  2:274256/1‑251 (MQ=255)                                                     tgtggtatgcgtttgaaacAAGAACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCa                                                                                                                                                                                     <  2:256390/231‑1 (MQ=255)                                                          tGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCgg                                                                                                                                                                              >  2:91897/1‑250 (MQ=255)                                                           gATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACAATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGAATCGTGGTTAGCAGGCGGTGGGTGGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTg                                                                                                                                                                             >  2:59643/1‑251 (MQ=255)                                                            atcgttttgacaaaGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGa                                                                                                                                                                            <  1:16294/239‑1 (MQ=255)                                                               gCTTTGATCACGCACCGCCATATTGACCGTTTCGCCCGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCCGTTCCGTGAGTg                                                                                                                                                                         >  2:364320/1‑251 (MQ=255)                                                                    aagcaacccccGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCcaca                                                                                                                                                                    <  2:63969/242‑1 (MQ=255)                                                                      tCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACt                                                                                                                                                                  >  2:246917/1‑251 (MQ=255)                                                                          ccACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCa                                                                                                                                                              <  2:105080/250‑1 (MQ=255)                                                                            aCCGCCATAATAAACGTCTCGACAGACTCTTCAATTTCACTGCGCAGGATAGGGGCAACAATCGCTAACAACTTACGGATCTGGAGAACGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGCAACCTTGCTGTTGCATCCTGGTTAGCAGGCGGAGGGTCGTGGAACTGGGTAAACCGGCTTGTTGCGCCAGATCCGTGAGTGCCACACTGCCaat                                                                                                                                                            >  2:371289/1‑249 (MQ=255)                                                                             ccGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTg                                                                                                                                                           <  1:348983/251‑1 (MQ=255)                                                                             ccGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATGGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAATCCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGTTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGGTCCGGGAGTGCCTCACTGCCATTg                                                                                                                                                           <  2:238602/251‑1 (MQ=255)                                                                              cGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTTCCACACTGCCATTgg                                                                                                                                                          >  2:289403/1‑251 (MQ=255)                                                                               gCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGa                                                                                                                                                         <  1:218564/251‑1 (MQ=255)                                                                                       gACCGTTTAGCCATACTCTTCCCTTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGACAGTTCCGTGAGTTCCACACTGCCATTGGATTCGGCa                                                                                                                                                  >  2:287060/1‑250 (MQ=255)                                                                                             ttCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAAGTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCAc                                                                                                                                            <  2:136652/250‑1 (MQ=255)                                                                                                         ttCCATTAGATTGCGCAGGATCGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTc                                                                                                                                >  2:371566/1‑250 (MQ=255)                                                                                                                aGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCAc                                                                                                                         <  2:363595/250‑1 (MQ=255)                                                                                                                               gggTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAAc                                                                                                         <  1:183421/251‑1 (MQ=255)                                                                                                                                   gAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTg                                                                                                      <  1:247667/250‑1 (MQ=255)                                                                                                                                          cGCTAACAAATTACGGCTCGGGATAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTc                                                                                               <  2:159758/250‑1 (MQ=255)                                                                                                                                              aaCAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTgg                                                                                          >  1:310283/1‑251 (MQ=255)                                                                                                                                                 aaaTTTCGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGGCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGc                                                                                       <  2:222187/251‑1 (MQ=255)                                                                                                                                                        ggCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGc                                                                                 >  1:113153/1‑250 (MQ=255)                                                                                                                                                                                                        cccAATGTCCCAGTTCGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTGC‑GCGTTTTGCGGGGATggggggggccctgac                                >  1:553/1‑238 (MQ=255)                                                                                                                                                                                                                       cGCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCTCTTTTTTCTg                <  1:89084/250‑1 (MQ=255)                                                                                                                                                                                                                        gCCAACCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCTCTTTTTTCTGTa              <  2:196425/251‑1 (MQ=255)                                                                                                                                                                                                                            aCCTGACGCACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTCAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTTTTTCGGGTTTCGAGGGAAAGTGGGCG‑GCCCAGACAGTCTCTTTTTTTTGtgtcct          >  1:278021/1‑246 (MQ=255)                                                                                                                                                                                                                                    cACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTTCCGCGTTTTTCCGGAATGGGGGCG‑ACCATGACAGTCCCtttttttttgtttgtgggaatat  >  2:194240/1‑234 (MQ=255)                                                                                                                                                                                                                                    cACGAAACCTTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGGTTTGCGGGGATGGGGGCGGACC‑TGACAGGCTCTTTTTTTTTTATCGTGGGAAtt   >  1:115362/1‑249 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    GCGCGGACATTCGCATCAGATGCGTA‑CACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCCATTAGATTGCGCAGGATAGGGTGAACAATCGCTAACAAATTACGGCTCTGGAGAAAGCTGCTGCCGACCATAAAGGCATGTGCGCCGATTGCCCAATGTCCCAGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA  >  NC_000913/4223184‑4223660

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