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%	133.4 / NA	42	*275* (TAG→TAA)	dinD	DNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (18/24);  total (18/24)

ACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCA  >  NC_000913/3818450‑3818715                                                                                                                                       |                                                                                                                                    acacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGa                                                                                                                                 >  2:781466/1‑139 (MQ=255)  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                >  2:147477/1‑139 (MQ=255)  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                <  2:365075/139‑1 (MQ=255)  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                >  2:482108/1‑139 (MQ=255)               gggTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtt                                                                                                                   >  2:1237422/1‑139 (MQ=255)                 gTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaa                                                                                                                         >  2:52656/1‑131 (MQ=255)                 gTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaa                                                                                                                         <  1:52656/131‑1 (MQ=255)                 gTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtttt                                                                                                                 <  2:1274650/139‑1 (MQ=255)                         gTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttgatt                                                                                                         <  1:482108/139‑1 (MQ=255)                         gTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttgatt                                                                                                         <  1:421284/139‑1 (MQ=255)                           gAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttgattga                                                                                                       >  2:312808/1‑139 (MQ=255)                                  aCCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACAccc                                                                                                >  1:1411795/1‑139 (MQ=255)                                        cAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCg                                                                                          <  2:1411795/139‑1 (MQ=255)                                             aCTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGc                                                                                     >  1:1136672/1‑139 (MQ=255)                                               ttGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                                                   <  1:147477/139‑1 (MQ=255)                                                tGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                                  <  2:1136672/139‑1 (MQ=255)                                                 ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAg                                                                                 >  2:144701/1‑139 (MQ=255)                                                     gAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTcc                                                                             <  2:648524/139‑1 (MQ=255)                                                         cATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                                         <  1:439376/139‑1 (MQ=255)                                                            gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                                      <  1:781466/139‑1 (MQ=255)                                                            gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                                      <  2:1346336/139‑1 (MQ=255)                                                                    agTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAacac                                                              <  2:676052/139‑1 (MQ=255)                                                                     gTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAAcaca                                                             <  1:1237422/139‑1 (MQ=255)                                                                     gTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAAcaca                                                             <  1:315598/139‑1 (MQ=255)                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTc                                                                              >  2:1353361/1‑121 (MQ=255)                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTc                                                                              <  1:1353361/121‑1 (MQ=255)                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttga                                                                                                           <  1:803083/90‑1 (MQ=255)                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttga                                                                                                           >  2:803083/1‑90 (MQ=255)                                                                                cGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                            >  2:103144/1‑129 (MQ=255)                                                                                cGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                            <  1:103144/129‑1 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaa                                                                                                                          <  2:402546/52‑1 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaa                                                                                                                          >  1:402546/1‑52 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                                                 >  1:1025728/1‑109 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                                                 <  2:1025728/109‑1 (MQ=255)                                                                                                     tGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt                             <  1:144701/139‑1 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAgg                                                                          >  2:1154576/1‑85 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAgg                                                                          <  1:1154576/85‑1 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                            >  2:129137/1‑99 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                            <  1:129137/99‑1 (MQ=255)                                                                                                              gtgtAAAAATTAC‑‑AGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAcc                  >  1:680053/1‑139 (MQ=255)                                                                                                                               aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGc   >  2:1372393/1‑139 (MQ=255)                                                                                                                                agaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCa  <  1:1372393/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    ACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCA  >  NC_000913/3818450‑3818715

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: ‑