Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,658	T→C	D46D (GAT→GAC)	creA →	PF05981 family protein CreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,658	0	T	C	100.0%	99.6 / NA	32	D46D (GAT→GAC)	creA	PF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/14);  total (18/14)

GGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGA  >  NC_000913/4635420‑4635888                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       gggCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGAtgtg                                                                                                                                                                                                                            <  1:69142/251‑1 (MQ=255)      tGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCGGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGaaa                                                                                                                                                                                                                         <  1:470648/249‑1 (MQ=255)            aaaTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCAc                                                                                                                                                                                                                 >  1:94436/1‑251 (MQ=255)              aTAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGCAGCCTTTGACGATCCCGATGGGAAAAATGTCCCCt                                                                                                                                                                                                               >  2:407159/1‑251 (MQ=255)              aTAAAAGAATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTCATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATACTGCTGGGCCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCACGGATCACAACATTGTTGTGGAAGCCTTTGCCGTTCTCGATGTGAAAAATGTCACCt                                                                                                                                                                                                               >  2:407153/1‑251 (MQ=255)                                  ttCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGAAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTTACCTGTTATGTGAGCCGGGCGaa                                                                                                                                                                                            >  1:448404/1‑250 (MQ=255)                                       aaTCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAAcagg                                                                                                                                                                                      <  1:86379/251‑4 (MQ=255)                                          cAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTTATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGTCCCGTATCTCAAAATTGTTGTTGAAGCCTTTGACTATCCCTATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAAcatgtg                                                                                                                                                                                    <  1:407153/250‑6 (MQ=255)                                          cAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAggtg                                                                                                                                                                                    <  1:407159/250‑1 (MQ=255)                                                  tATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTaa                                                                                                                                                                             >  1:260620/1‑249 (MQ=255)                                                   aTTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAgg                                                                                                                                                                          >  1:480406/1‑251 (MQ=255)                                                      tACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGTCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAgggg                                                                                                                                                                        <  2:502590/250‑1 (MQ=255)                                                            gACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAt                                                                                                                                                                                >  1:492728/1‑236 (MQ=255)                                                            gACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTAt                                                                                                                                                                                <  2:492728/236‑1 (MQ=255)                                                                           cGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCATGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGGGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATAc                                                                                                                                                  >  1:155404/1‑251 (MQ=255)                                                                             tgtAATAGATAACAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAACCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACGGTATTTAAAATGATCGGCCCGGATCACAAGATTGTTGTGGAAGCCTTTGACGTTCCCGTTGGGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAACCAGGTGGGATAAAAGGGGGATTGGGTCTGGCGGAAGATAcac                                                                                                                                                >  2:406292/1‑249 (MQ=255)                                                                                          aaTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGctat                                                                                                                                   >  1:151905/1‑248 (MQ=255)                                                                                               tAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGACATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCATTTGACGATCCCGATGTGCAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGAGGTATTAAGGGGGGATTGGGTCTGGCGGAAGATACTTCCGATGAGGATATTTCtt                                                                                                                              >  1:347903/1‑251 (MQ=255)                                                                                                 aCAATATGAAGTACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGAAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGt                                                                                                                            >  2:317702/1‑251 (MQ=255)                                                                                                      atGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCa                                                                                                                       <  2:480406/251‑1 (MQ=255)                                                                                                               aaGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTTGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGcc                                                                                                              <  1:105852/251‑1 (MQ=255)                                                                                                                   aTTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGAtt                                                                                                          <  1:64243/251‑1 (MQ=255)                                                                                                                        aTCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGaaa                                                                                                                                                                                                                         <  1:449458/135‑1 (MQ=255)                                                                                                                        aTCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGaaa                                                                                                                                                                                                                         >  2:449458/1‑135 (MQ=255)                                                                                                                                tttAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATc                                                                                             >  1:461606/1‑251 (MQ=255)                                                                                                                                      ccTGATAATGCTGGGGCCATTGGCTCAAGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTAtt                                                                                        >  1:440897/1‑250 (MQ=255)                                                                                                                                                      ccATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTc                                                                        >  1:479289/1‑250 (MQ=255)                                                                                                                                                         ttGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAggg                                                                    >  1:36057/1‑251 (MQ=255)                                                                                                                                                          tGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAggg                                                                    <  1:402690/250‑1 (MQ=255)                                                                                                                                                                                             gTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCtt                                <  1:504660/251‑1 (MQ=255)                                                                                                                                                                                                       tGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGttt                      >  1:290688/1‑250 (MQ=255)                                                                                                                                                                                                                           aTTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTatga  <  1:113258/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                       GGGCTTGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGA  >  NC_000913/4635420‑4635888

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