Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,817,594	A→G	intergenic (+55/‑166)	yicC → / → dinD	UPF0701 family 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,594	0	A	G	100.0%	98.7 / NA	28	intergenic (+55/‑166)	yicC/dinD	UPF0701 family protein YicC/DNA damage‑inducible protein D
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/13);  total (15/13)

CCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTA  >  NC_000913/3817463‑3817703                                                                                                                                    |                                                                                                              ccGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCAtt                                                                                                       >  2:226706/1‑140 (MQ=255)     aGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTttc                                                                                                   >  2:372712/1‑140 (MQ=255)              cTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCa                                                                                          <  2:207734/140‑1 (MQ=255)                 cGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACa                                                                                       >  1:17510/1‑140 (MQ=255)                                    ttgattgaGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTg                                                                    >  1:347046/1‑140 (MQ=255)                                    ttgattgaGCAGATGCGCGAGCAGATTCAGAACATAGAATAACTCGTAACGCCACTTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTg                                                                    >  2:15808/1‑140 (MQ=255)                                     tgattgaGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGc                                                                   <  1:226706/140‑1 (MQ=255)                                      gattgaGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGcc                                                                  <  2:376359/140‑1 (MQ=255)                                                                    cATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGccc                                    <  2:17510/140‑1 (MQ=255)                                                                      tCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAg                                  >  1:395757/1‑140 (MQ=255)                                                                      tCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAg                                  >  1:184953/1‑140 (MQ=255)                                                                       cGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGc                                 >  2:192133/1‑140 (MQ=255)                                                                                cGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTgaga                        <  2:380249/140‑1 (MQ=255)                                                                                   aaCGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCATACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGcaca                              <  2:344952/131‑1 (MQ=255)                                                                                   aaCGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGcaca                              >  1:344952/1‑131 (MQ=255)                                                                                   aaCGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACg                     >  2:91769/1‑140 (MQ=255)                                                                                   aaCGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACg                     >  1:373929/1‑140 (MQ=255)                                                                                    aCGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGt                    >  2:311594/1‑140 (MQ=255)                                                                                      gCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGca                                >  1:106989/1‑126 (MQ=255)                                                                                      gCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGca                                <  2:106989/126‑1 (MQ=255)                                                                                       ccAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTgag                 <  1:192133/140‑1 (MQ=255)                                                                                       ccAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTgag                 <  2:294829/140‑1 (MQ=255)                                                                                       ccAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTgag                 <  1:311594/140‑1 (MQ=255)                                                                                                 cTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACaa       <  2:184953/140‑1 (MQ=255)                                                                                                 cTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACaa       <  1:56411/140‑1 (MQ=255)                                                                                                   tttCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACt     >  1:215592/1‑140 (MQ=255)                                                                                                    ttCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTg    <  2:215592/140‑1 (MQ=255)                                                                                                    ttCCGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTg    >  1:2357/1‑140 (MQ=255)                                                                                                      ccGCCTTCACAAATGCCGCCACTCAAACAGGGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGta  >  2:304914/1‑140 (MQ=255)                                                                                                                                    |                                                                                                              CCGAAGTGACAAACTCCGCCATCGAGCTGAAAGTGTTGATTGAGCAGATGCGCGAGCAGATTCAGAACATCGAATAACTCGTAACGCCAATTCTTACTTTTCCGCCTTCACAAATGCCGCCACTCAAACAGAGCGGCATTTTTCTTCCCCGCAACATTCAATTCTGTTTTGCGTGCCTGCTCCAGATTTTGCGATGTTTTTTTGCCCAGCACACTGAGAACGTGAGATACTCACAACTGTA  >  NC_000913/3817463‑3817703

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