Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,978,243	G→T	intergenic (‑46/+275)	flhD ← / ← insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,978,243	0	G	T	100.0%	47.6 / NA	19	intergenic (‑46/+275)	flhD/insB5	DNA‑binding transcriptional dual regulator FlhD/IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (7/12);  total (7/12)

CTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTC  >  NC_000913/1978030‑1978451                                                                                                                                                                                                                      |                                                                                                                                                                                                                 cTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAgn                                                                                                                                                                                          <  1:78553/238‑2 (MQ=255)             cGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCAtt                                                                                                                                                                   >  1:73811/1‑249 (MQ=255)                    gTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGACCAATCAAACGCTGTGTCAGTAGTCAATATGACACGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGCATAACCAACTTTATTTTTATGCGGTTTCAACTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGTATAGtt                                                                                                                                                          >  1:299722/1‑251 (MQ=255)                             ncTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGc                                                                                                                                                 >  1:82465/2‑251 (MQ=255)                                               tcttcATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTGTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCTGTAAAAAATAATGAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAacca                                                                                                                                <  1:133582/250‑1 (MQ=255)                                                                        tANCGGNCNNTNNGTCNTGAACAATCAAACGCTGTGCANNTAGTAAAGANGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCNACCCAGAATNNCCNACTTGNTNNNTNNNCGNTTTCACCTCACCCCGTGATNTCNNCGGGAANCCCCGGTAAAAAATAAGTAGCATTAg                                                                                                                                                                 <  2:73811/192‑1 (MQ=255)                                                                                             cAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGTGTCTTTacac                                                                                 <  1:74612/251‑1 (MQ=255)                                                                                                         gtgCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAGCCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATggg                                                                                          <  1:325934/230‑1 (MQ=255)                                                                                                          caccgcatggcaaTATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTGCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGAttt                                                                    <  1:139952/240‑1 (MQ=255)                                                                                                             aaGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGGAATGTGTCGTTACACTTATCTAAGATTTTTc                                                                 <  1:112457/251‑1 (MQ=255)                                                                                                              agtagtAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTcn                                                                <  1:261666/251‑2 (MQ=255)                                                                                                                           aCAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTGTATGCGGTTTCACCTCACCCCGTGATGTCg                                                                                                                                                                                                   <  1:257262/107‑1 (MQ=255)                                                                                                                                       cATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCg                                                                                                                                                                                                           <  1:132732/87‑1 (MQ=255)                                                                                                                                                      gCAACTCGGAGGTATGCATTATTCCCNCCCAGANNAGNCAACTTTATTTTTATGCGGNTTNNNCNNNCNACGTGNTGNNGCCGGGAAGNCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGNACCACCTTTNNGGTTTAATATGTCCTTACAAANAGANNANNGNCttntcn                                                                                  >  2:77821/1‑189 (MQ=255)                                                                                                                                                          cTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgg                    >  1:23264/1‑251 (MQ=255)                                                                                                                                                          cTCGGAGGTATGCATTATTCCCACCCAGAATAACCAAATTTATTTTTATACGGTTTCACCTAAACCCGTGATCTCGACGGGAGGACCCGGTAACAAATAATAAGCATTAGAATAGTTCCCATACGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATg                                                                                            >  1:218126/1‑179 (MQ=255)                                                                                                                                                                     gCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTa                                                                                                        >  1:135194/1‑156 (MQ=255)                                                                                                                                                                             cccacccaGAATAACCAACTTTATTTTTATGCGGTTTCACCTCCCCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAAGAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGTGGGCTGCGATTTTc  <  1:77821/250‑1 (MQ=255)                                                                                                                                                                                                           gCGGTTTCACCTCACCCCGTGATGTCGCCGGGACGACCCGGTAAAAAATAATTAGCATTAGAATAGTTGGGATAAGCTGCAATAAGCAGAACCCCCTTTTTGGTTTAAACTGTACTTACAAATAGAACAGGGTCTTTACACTTATCt                                                                           >  1:269738/1‑147 (MQ=255)                                                                                                                                                                                                              gTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAac                                                                                                                                  <  1:190200/89‑1 (MQ=255)                                                                                                                                                                                                                      |                                                                                                                                                                                                                 CTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTC  >  NC_000913/1978030‑1978451

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