Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,817,618	Δ1 bp	intergenic (+79/‑142)	yicC → / → dinD	putative RNase adaptor protein YicC/DNA damage‑inducible protein D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,817,618	0	A	.	100.0%	90.9 / NA	20	intergenic (+79/‑142)	yicC/dinD	putative RNase adaptor protein YicC/DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (8/12);  total (8/12)

AAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCC‑ACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGC  >  NC_000913/3817375‑3817815                                                                                                                                                                                                                                                     |                                                                                                                                                                                                      aaaaaaGCAGCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAAtt                                                                                                                                                                                                >  2:213254/1‑251 (MQ=255)        aagCCGTTGGTCGTCGTCTGGACTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGtttt                                                                                                                                                                                         >  2:288378/1‑251 (MQ=255)                                 tATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCaa                                                                                                                                                                                                  >  1:427766/1‑217 (MQ=255)                                        cAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTGGCGATGtttttt                                                                                                                                                         <  1:368457/251‑1 (MQ=255)                                                  aCCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGcac                                                                                                                                               <  2:89850/251‑1 (MQ=255)                                                                     tGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACt                                                                                                                            >  1:374473/1‑251 (MQ=255)                                                                        gTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGGAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTGCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTcac                                                                                                                         <  1:159501/251‑1 (MQ=255)                                                                             aGTCTATCAATGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTg                                                                                                                    <  2:263765/251‑1 (MQ=255)                                                                                     aaTGCCGAAGTGACAAACTCCGCGATTGAGCTGAAAGTGTTGATTGAGCAGATTCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATGTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATaaa                                                                                                            <  1:56867/251‑1 (MQ=255)                                                                                                           cgattgAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAATATTGAGTAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTActttct                                                                                      <  1:58230/246‑1 (MQ=255)                                                                                                                               attgaGCAGATGCGCGAGCAGATTCAGAATATTGAATAACGCGTAACACTAAGTCTTACTTGTCAGCCTTCACAAGTGCCGCCCACTCAAACAGAGCGGCATTTTGCTTCCCCGCAAC‑TTCAATGCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTTAGATACTCACAACTGTATATAAATACAGGTACAGAATTACTTCCTGTGCAATCAATATCCCATgg                                                                  <  2:389470/251‑1 (MQ=255)                                                                                                                                    gCAGATGCTCGAGCAGATTCAGAATATTGAGGAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTGTTTTTGCCCAGCCCACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATGTACTTTCTTTGCAATCGATATCACATGGAGTgg                                                             <  1:141969/251‑1 (MQ=255)                                                                                                                                                          aatattgagtAACTCGTAACACTAATTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAAc                                       >  2:193880/10‑251 (MQ=255)                                                                                                                                                                           aacactaaTTCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAAATacgg                      <  1:107278/245‑4 (MQ=255)                                                                                                                                                                                    tCTTACTTTTCAGCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAAATacggcattatgg              <  1:193880/250‑12 (MQ=255)                                                                                                                                                                                             tcagCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTCTGAAGAAATacggcattatggcactgaagg     <  1:165763/247‑21 (MQ=255)                                                                                                                                                                                                gCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAAATacggcattatggcactgaagggca  >  2:328706/1‑227 (MQ=255)                                                                                                                                                                                                gCCTTCACAAATGCCGCCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAAATacggcattatggcactgaagggca  >  2:354197/1‑227 (MQ=255)                                                                                                                                                                                                                                          ttcCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATgaac                                                    <  2:135247/158‑1 (MQ=255)                                                                                                                                                                                                                                          ttcCCCGCAAC‑TTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATCGATATCACATGGAGTGGGCAATgaac                                                    >  1:135247/1‑158 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                      AAAAAAGAAGCGGTTGGTCGTCGTCTGGATTTTATGATGCAGGAGTTCAACCGCGAGTCGAACACTCTTGCGTCGAAGTCTATCAATGCCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCC‑ACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTATATAAATACAGTTACAGATTTACTTTCTTTGCAATTGATATCACATGGAGTGGGCAATGAACGAACATCATCAACCTTTTGAAGAGATAAAACTGATTAATGCAAACGGAGC  >  NC_000913/3817375‑3817815

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