Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,664,969	(T)8→7	noncoding (106/106 nt)	gadY →	small regulatory RNA GadY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,664,969	0	T	.	100.0%	72.0 / NA	16	noncoding (106/106 nt)	gadY	small regulatory RNA GadY
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (7/9);  total (7/9)

TTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTTTATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATCCACAGGATACTGCAACTCGCTTAATTGAAAAACCATGTATAACAATAAGTTGC  >  NC_000913/3664743‑3665214                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      ttttGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCt                                                                                                                                                                                                                              >  2:265107/1‑251 (MQ=255)        tAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGAGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTAttcct                                                                                                                                                                                                                       >  1:43158/1‑251 (MQ=255)                       aataaGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCa                                                                                                                                                                                                        <  1:253088/251‑1 (MQ=255)                           aGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGc                                                                                                                                                                                                    >  2:369694/1‑251 (MQ=255)                                        cGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTccc                                                                                                                                                                                       <  1:369694/251‑1 (MQ=255)                                                        caGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACtt                                                                                                                                                                       <  2:40377/251‑1 (MQ=255)                                                                          tataAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTc                                                                                                                                                     >  2:300026/1‑251 (MQ=255)                                                                                          aTCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCg                                                                                                                                     >  2:213074/1‑251 (MQ=255)                                                                                              tATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATc                                                                                                                                  <  1:304649/250‑1 (MQ=255)                                                                                                             ttATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCtgt                                                                                                                  >  2:77918/1‑251 (MQ=255)                                                                                                               atatTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTgg                                                                                                                <  1:45759/251‑1 (MQ=255)                                                                                                                               gagCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCataata                                                                                                <  1:10406/251‑1 (MQ=255)                                                                                                                                                                   aaCGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCtgt                                                             <  1:250746/250‑1 (MQ=255)                                                                                                                                                                        tttATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATc                                                       >  2:125808/1‑251 (MQ=255)                                                                                                                                                                                                               ccccTATGCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATCCACAGGATACTGCAACTCGCTTAATTGAAAAACCATGta                <  1:265107/251‑1 (MQ=255)                                                                                                                                                                                                                      gCCGGGTTTTTTT‑ATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATCCACAGGATACTGCAACTCGCTTAATGGAAAAACCATGTATAACAAt         <  1:77918/251‑1 (MQ=255)                                                                                                                                                                                                                               tttttATGTCTGAGTAAAAATCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCACCAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATCCACAGGATACTGCAACTCGCTTAATTGAAAAACCATGTATAACAATAAGTTGc  >  2:177432/1‑250 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      TTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCTGATCTTATTTCCAGTAAAAGTTATATTTAACTTACTGAGAGCACAAAGTTTCCCGTGCCAACAGGGAGTGTTATAACGGTTTATTAGTCTGGAGACGGCAGACTATCCTCTTCCCGGTCCCCTATGCCGGGTTTTTTTTATGTCTGAGTAAAACTCTATAATCTTATTCCTTCCGCAGAACGGTCAGTGCCGTAAAAATTCCCTTGCGCAACAATACTTGCCGCCGAGTCACGGTTCGACAATCTCTGCGCCGATCGCTCCTGATACTCTGTGGGCGTCATCCCATAATAATTTCGAAAGACGTAAATGAAATACGACACGCTGTGATATCCACAGGATACTGCAACTCGCTTAATTGAAAAACCATGTATAACAATAAGTTGC  >  NC_000913/3664743‑3665214

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