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%	63.5 / NA	18	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/11);  total (7/11)

TGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG  >  NC_000913/1978013‑1978458                                                                                                                                                                                                                                       |                                                                                                                                                                                                                        tGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCgg                                                                                                                                                                                                      <  2:160729/250‑1 (MQ=255) tGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACGTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCgg                                                                                                                                                                                                      <  1:283767/250‑1 (MQ=255)       gTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGccc                                                                                                                                                                                               <  2:15592/251‑1 (MQ=255)             gCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTa                                                                                                                                                                                          >  1:58942/1‑250 (MQ=255)                aGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTaaaaa                                                                                                                                                                                      <  1:50736/251‑1 (MQ=255)                                  agagTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATGTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATa                                                                                                                                                                    <  1:58291/251‑1 (MQ=255)                                                   cGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCaa                                                                                                                                                   <  1:41651/251‑1 (MQ=255)                                                                    cATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCttt                                                                                                                                  <  1:233884/251‑1 (MQ=255)                                                                                                  cTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGaact                                                                                                    >  2:359592/1‑249 (MQ=255)                                                                                                     tGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATgg                                                                                                  <  1:165211/250‑1 (MQ=255)                                                                                                                      cGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGCATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTa                                                                                 >  2:213729/1‑250 (MQ=255)                                                                                                                                               aGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATGGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTg                                                       >  2:318514/1‑251 (MQ=255)                                                                                                                                                             aTGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCTATAAGCTGCAATAAg                                                                                                                                               <  1:107250/149‑1 (MQ=255)                                                                                                                                                             aTGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAg                                                                                                                                               >  2:107250/1‑149 (MQ=255)                                                                                                                                                                          aCTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCCCCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACg                            >  2:60234/1‑251 (MQ=255)                                                                                                                                                                            tCGGAGGTATGCATTATTCCCACACAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGCTAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACgg                           >  2:406679/1‑250 (MQ=255)                                                                                                                                                                               gAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACGCGCACATACAACGGAggg                       <  2:134478/251‑1 (MQ=255)                                                                                                                                                                                                    caGAATAACCAACTTTATTTTTATGCGGTTTCACCTCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATg  <  1:196981/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                        TGATTGGTTTCTGCCAGCTTAACCATTTGCGGAAGAGTCAGTGCCGCTAACGTTGTCGCCATTTCTTCATTTATGCCGAGACGAAACATAGCGGACGCTTTGTCCTGAACAATCAAACGCTGTGCAAGTAGTAAATATGACAAGTTGATGTCATAAATGTGTTTCAGCAACTCGGAGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATG  >  NC_000913/1978013‑1978458

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