Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,647,725	A→G	R42R (CGT→CGC)	dicC ←	Qin prophage; DNA‑binding transcriptional regulator for DicB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,647,725	0	A	G	100.0%	184.3 / NA	51	R42R (CGT→CGC)	dicC	Qin prophage; DNA‑binding transcriptional regulator for DicB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/22);  total (29/22)

GGCATCCAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAAC  >  NC_000913/1647595‑1647857                                                                                                                                   |                                                                                                                                     ggCATCCAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATcgcg                                                                                                                             >  2:148704/1‑140 (MQ=255)     cccAGACCCCCACGCGCCTGTTCAGGAGCGTTTTTCATTTTCCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGAc                                                                                                                         <  2:109080/139‑1 (MQ=255)      ccAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGAcc                                                                                                                        >  2:338234/1‑140 (MQ=255)      ccAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGAcc                                                                                                                        >  1:72792/1‑140 (MQ=255)       cAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCt                                                                                                                       >  1:8114/1‑140 (MQ=255)                aCGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCggg                                                                                                                  <  1:898357/136‑1 (MQ=255)                aCGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCggg                                                                                                                  >  2:898357/1‑136 (MQ=255)                    gCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAAt                                                                                                          >  1:1012881/1‑140 (MQ=255)                      ttGTTAAGGAGTTATTTCCTTTGTACAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATcc                                                                                                        <  2:586307/139‑1 (MQ=255)                       tGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCTGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCGTCCTGTACGCGCATCGCGCGACCTTCGGGAACTAAATccc                                                                                                       >  2:922122/1‑140 (MQ=255)                        tttCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATcccc                                                                                                      <  1:338234/139‑1 (MQ=255)                              gAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCa                                                                                                >  2:844176/1‑140 (MQ=255)                                   aTTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTc                                                                                                  >  1:309540/1‑133 (MQ=255)                                   aTTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTc                                                                                                  <  2:309540/133‑1 (MQ=255)                                   aTTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCtat                                                                                           >  2:29029/1‑140 (MQ=255)                                                     ccccgcccGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGaa                                                                         >  2:804234/1‑140 (MQ=255)                                                       ccgcccgcTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGaa                                                                         <  1:657034/138‑1 (MQ=255)                                                       ccgcccgcTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGaa                                                                         >  2:657034/1‑138 (MQ=255)                                                              cTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACctgctg                                                                <  1:652557/140‑1 (MQ=255)                                                                                atcatAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATG‑‑‑CCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTg                                           >  1:246995/1‑140 (MQ=255)                                                                                  catAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGttttt                                            >  1:850122/1‑140 (MQ=255)                                                                                   atAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATG‑‑‑CCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAAc                                        <  2:246995/140‑1 (MQ=255)                                                                                     aaaCTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGaa                                         >  2:472795/1‑140 (MQ=255)                                                                                     aaaCTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGaa                                         >  1:929200/1‑140 (MQ=255)                                                                                      aaCTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAAc                                        <  1:148704/140‑1 (MQ=255)                                                                                      aaCTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAAc                                        <  1:988048/140‑1 (MQ=255)                                                                                      aaCTTTGGGATCATACTCAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAAc                                        <  2:929200/140‑1 (MQ=255)                                                                                        cTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCg                                      <  1:844176/140‑1 (MQ=255)                                                                                          ttGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGaa                                    <  2:72792/140‑1 (MQ=255)                                                                                           tGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCg                                      >  2:881304/1‑137 (MQ=255)                                                                                           tGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCg                                      <  1:881304/137‑1 (MQ=255)                                                                                           tGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGaaa                                   >  1:155500/1‑140 (MQ=255)                                                                                                  aTACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACaaaa                            >  1:915895/1‑140 (MQ=255)                                                                                                   tACTGAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAag                           >  1:873224/1‑140 (MQ=255)                                                                                                       gAAGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCg                       >  2:75222/1‑140 (MQ=255)                                                                                                         aGCTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTc                     >  1:774360/1‑140 (MQ=255)                                                                                                           cTCCCCGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAg                   >  1:664323/1‑140 (MQ=255)                                                                                                             ccccGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGtt                 >  1:108008/1‑140 (MQ=255)                                                                                                             ccccGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGtt                 >  1:989649/1‑140 (MQ=255)                                                                                                             ccccGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGtt                 >  2:412597/1‑140 (MQ=255)                                                                                                              cccGCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGttt                >  2:752514/1‑140 (MQ=255)                                                                                                                 gCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTaa             <  1:752514/140‑1 (MQ=255)                                                                                                                 gCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTaa             <  1:855166/140‑1 (MQ=255)                                                                                                                 gCCAGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTaa             <  2:155500/140‑1 (MQ=255)                                                                                                                    aGATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCa          <  2:669200/140‑1 (MQ=255)                                                                                                                     gATGCCTCCTGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCAt         >  2:541667/1‑140 (MQ=255)                                                                                                                           tcctGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTaaaa   <  1:72600/140‑1 (MQ=255)                                                                                                                           tcctGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTaaaa   <  1:412597/140‑1 (MQ=255)                                                                                                                           tcctGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTaaaa   <  1:472795/140‑1 (MQ=255)                                                                                                                           tcctGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTaaaa   <  2:108008/140‑1 (MQ=255)                                                                                                                            cctGTAGGCGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAac  >  1:212382/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                     GGCATCCAGAGCCTCACGAGCCTGTTCAGGAGTGATTTTCATTGTTCAACCGCCCCGCCCGCTTCGTCTTACGATATTCATCATAAACTTTGGGATCATACTGAAGCTCCCCGCCAGATGCCTCCTGTAGACGCATCGCGCGACCTTCGGGAACTAAATCCCCTTTCCAGCTATAAAGCGAAGCCAAACGAATACCTGCTGCTTGTGCAAGTTTTGTTTTTGAACCGAAATACAAAAGAGCGTCAGTTTTAAGCATTTAAAAC  >  NC_000913/1647595‑1647857

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