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

CACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAC  >  NC_000913/3818449‑3818698                                                                                                                                        |                                                                                                                   cacacaCGGTGACGTGAGTCGTAAGGAGAGGCAAACCATTCAGTAACGTGGCGGAACCAGGCCTTAGGAGTTGGCGACTGCGCAGGTCAGCATCAAGCAGTTGGAAAATAGTGTAAAAATTACAGAGAGGAAATAACAg                                                                                                                 <  1:1026756/139‑1 (MQ=255)   cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                               >  1:787063/1‑139 (MQ=255)    acacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaaa                                                                                                              >  1:329414/1‑139 (MQ=255)     cacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAAc                                                                                                             >  1:227190/1‑139 (MQ=255)         ttGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaa                                                                                                         >  1:323779/1‑139 (MQ=255)                 ggTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCttt                                                                                                 >  2:383549/1‑139 (MQ=255)                              ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                                    <  2:323779/139‑1 (MQ=255)                                    ccATTCATGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGGAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  1:798156/1‑139 (MQ=255)                                    ccATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  2:826139/139‑1 (MQ=255)                                                  ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  2:890813/1‑125 (MQ=255)                                                  ggCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  1:890813/125‑1 (MQ=255)                                                   gCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAgg                                                               <  2:787063/139‑1 (MQ=255)                                                          cATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                        <  2:329414/139‑1 (MQ=255)                                                            tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  1:735180/115‑1 (MQ=255)                                                            tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  2:735180/1‑115 (MQ=255)                                                            tGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGc                                                      <  1:383549/139‑1 (MQ=255)                                                             gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                     <  2:227190/139‑1 (MQ=255)                                                                 gaggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGa                                                 >  1:1157081/1‑139 (MQ=255)                                                                  aggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              >  2:308628/1‑109 (MQ=255)                                                                  aggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                              <  1:308628/109‑1 (MQ=255)                                                                  aggagTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGaa                                                >  2:1074326/1‑139 (MQ=255)                                                                     agTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAacac                                             <  1:7015/139‑1 (MQ=255)                                                                      gTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAAcaca                                            <  2:169967/139‑1 (MQ=255)                                                                       ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                           >  2:796729/1‑139 (MQ=255)                                                                        tGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAAt                                          >  2:415849/1‑139 (MQ=255)                                                                         gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGc                                                      >  1:589750/1‑126 (MQ=255)                                                                         gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGc                                                      <  2:589750/126‑1 (MQ=255)                                                                             agcccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTa                                     <  1:980983/137‑1 (MQ=255)                                                                                       tCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCAc                           <  1:546129/139‑1 (MQ=255)                                                                                         aGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACct                         >  2:428642/1‑139 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                 >  2:1109801/1‑92 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                 <  1:1109801/92‑1 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                        >  1:288371/1‑101 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                        <  2:288371/101‑1 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAG‑CAAGGTCCAGGTGCAGTGAacac                                             <  2:500921/111‑1 (MQ=255)                                                                                                aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAG‑CAAGGTCCAGGTGCAGTGAacac                                             >  1:500921/1‑111 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAg                                                                                                                 >  2:799143/1‑39 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAg                                                                                                                 <  1:799143/39‑1 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTa                                                                       <  2:378786/81‑1 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTa                                                                       >  1:378786/1‑81 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTccccc             >  1:1174388/1‑139 (MQ=255)                                                                                                     ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTccccc             >  2:657969/1‑139 (MQ=255)                                                                                                      tGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAGACTCTTTTTGTTGAGTGAGACACCCGGTGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            <  2:198074/139‑1 (MQ=255)                                                                                                       ggAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCATTGAACACAATGGCTACATATTGCACCTTTCTTTTCCCCCt            <  2:743630/138‑1 (MQ=255)                                                                                                       ggAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            >  1:743630/1‑138 (MQ=255)                                                                                                              aGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTg                                                       <  1:1149393/88‑1 (MQ=255)                                                                                                              aGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTg                                                       >  2:1149393/1‑88 (MQ=255)                                                                                                               gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAg                                                          >  2:270362/1‑84 (MQ=255)                                                                                                               gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAg                                                          <  1:270362/84‑1 (MQ=255)                                                                                                               gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTaa   >  2:188157/1‑139 (MQ=255)                                                                                                                tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            >  1:877918/1‑129 (MQ=255)                                                                                                                tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            >  1:1151305/1‑129 (MQ=255)                                                                                                                tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            <  2:877918/129‑1 (MQ=255)                                                                                                                tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            <  2:1151305/129‑1 (MQ=255)                                                                                                                tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAc  <  1:415849/139‑1 (MQ=255)                                                                                                                                aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            >  2:105823/1‑113 (MQ=255)                                                                                                                                aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt            <  1:105823/113‑1 (MQ=255)                                                                                                                                   aaaTAACAGAAACTAAAAACTCTTTTTGTTGATTTAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                      >  2:544693/1‑84 (MQ=255)                                                                                                                                   aaaTAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                      <  1:544693/84‑1 (MQ=255)                                                                                                                                        |                                                                                                                   CACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAC  >  NC_000913/3818449‑3818698

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