Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,818,584	G→A	*275* (TAG→TAA)	dinD →	DNA damage‑inducible protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,818,584	0	G	A	100.0%	129.4 / NA	41	*275* (TAG→TAA)	dinD	DNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)

ACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAC  >  NC_000913/3818452‑3818698                                                                                                                                     |                                                                                                                   acacTTTGACGTGGGTCACAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaaa                                                                                                              >  2:22237/1‑139 (MQ=255)                 cGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttg                                                                                              <  2:310903/139‑1 (MQ=255)                   cAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttgtt                                                                                            >  2:122746/1‑139 (MQ=255)                          aGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgag                                                                                     >  1:199277/1‑139 (MQ=255)                           ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                                    <  2:118048/139‑1 (MQ=255)                           ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                                    <  1:22237/139‑1 (MQ=255)                           ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                                    <  2:761313/139‑1 (MQ=255)                                  cATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGa                                                                             >  1:219720/1‑139 (MQ=255)                                    ttCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATg                                                                           >  1:154824/1‑139 (MQ=255)                                           aCTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGc                                                                    <  2:219720/139‑1 (MQ=255)                                            cTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGcc                                                                   <  2:703213/139‑1 (MQ=255)                                            cTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGcc                                                                   <  1:1119887/139‑1 (MQ=255)                                              tGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                 >  2:919588/1‑139 (MQ=255)                                               ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  1:1410848/125‑1 (MQ=255)                                               ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  2:1410848/1‑125 (MQ=255)                                               ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAg                                                                >  2:197467/1‑139 (MQ=255)                                               ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAg                                                                >  2:1379040/1‑139 (MQ=255)                                                   gAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTcc                                                            >  2:324642/1‑139 (MQ=255)                                                        aTGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTg                                                       <  2:394598/139‑1 (MQ=255)                                                         tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGc                                                      >  2:937609/1‑139 (MQ=255)                                                          gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                     >  1:1091505/1‑139 (MQ=255)                                                           ccTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAg                                                    <  1:1379040/139‑1 (MQ=255)                                                                ggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACt                                                                                                            <  2:1119789/78‑1 (MQ=255)                                                                ggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACt                                                                                                            >  1:1119789/1‑78 (MQ=255)                                                                    ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTc                                                             >  1:43703/1‑121 (MQ=255)                                                                    ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTc                                                             <  2:43703/121‑1 (MQ=255)                                                                         gACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                      <  1:919588/139‑1 (MQ=255)                                                                         gACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                      <  1:122746/139‑1 (MQ=255)                                                                           ccccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtttt                                                                                                <  2:1441095/79‑1 (MQ=255)                                                                           ccccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtttt                                                                                                >  1:1441095/1‑79 (MQ=255)                                                                                                  ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  2:1013854/74‑1 (MQ=255)                                                                                                  ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  1:1013854/1‑74 (MQ=255)                                                                                                  ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTccccc             <  2:666641/139‑1 (MQ=255)                                                                                                  ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTccccc             >  1:391746/1‑139 (MQ=255)                                                                                                            gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAgg                                                               >  1:257951/1‑79 (MQ=255)                                                                                                            gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAgg                                                               <  2:257951/79‑1 (MQ=255)                                                                                                            gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTaa   >  2:89877/1‑139 (MQ=255)                                                                                                             tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTTCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAc  <  1:430366/139‑1 (MQ=255)                                                                                                             tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAc  <  1:937609/139‑1 (MQ=255)                                                                                                                             aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                                <  1:183012/77‑1 (MQ=255)                                                                                                                             aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                                >  2:183012/1‑77 (MQ=255)                                                                                                                                     |                                                                                                                   ACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAC  >  NC_000913/3818452‑3818698

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