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

ACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCATT  >  NC_000913/3818450‑3818717                                                                                                                                       |                                                                                                                                      acacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGa                                                                                                                                   <  2:1356060/139‑1 (MQ=255)  cacacTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                  <  2:240312/139‑1 (MQ=255)       tttGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaa                                                                                                                             >  2:117050/1‑139 (MQ=255)          gACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGCACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaaa                                                                                                                          >  2:178580/1‑139 (MQ=255)          gACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaaa                                                                                                                          >  2:1496838/1‑139 (MQ=255)          gACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTaaaaa                                                                                                                          >  1:1278569/1‑139 (MQ=255)                   cGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttg                                                                                                                 <  2:1162052/139‑1 (MQ=255)                     cAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTttgtt                                                                                                               >  2:547737/1‑139 (MQ=255)                         gTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGttgatt                                                                                                           <  2:1278569/139‑1 (MQ=255)                            aGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgag                                                                                                        >  1:1442312/1‑139 (MQ=255)                            aGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgag                                                                                                        >  2:444207/1‑139 (MQ=255)                             ggCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTgaga                                                                                                       <  1:547737/139‑1 (MQ=255)                                              cTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGcc                                                                                      <  1:444207/139‑1 (MQ=255)                                               ttGGCGGAACCATGCCTGAGGAGTTGCAGACCCCGCAGGTCAGCAGCAAGCAGTTGCAAATCAGTGTAAAAATCACAGAGAAGACATAACAGAAACGAAAAACTCTTTTTGTTGATTGAGACACCCGCTGCGTAAGCCa                                                                                     <  1:178580/139‑1 (MQ=255)                                                      aaCCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                  <  2:345766/87‑1 (MQ=255)                                                      aaCCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGaa                                                                                                                                  >  1:345766/1‑87 (MQ=255)                                                      aaCCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCa                                                                              <  1:1496838/139‑1 (MQ=255)                                                            gCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCa                                                                        >  1:369853/1‑139 (MQ=255)                                                             ccTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                                    >  1:207129/1‑126 (MQ=255)                                                             ccTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCaa                                                                                    <  2:207129/126‑1 (MQ=255)                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                              >  1:365680/1‑139 (MQ=255)                                                                      ttGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACaa                                                              >  1:788715/1‑139 (MQ=255)                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAgg                                                                                  >  2:1203578/1‑117 (MQ=255)                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAgg                                                                                  <  1:1203578/117‑1 (MQ=255)                                                                        gCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATg                                                            >  2:1308467/1‑139 (MQ=255)                                                                              cccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                                           >  1:1069371/1‑118 (MQ=255)                                                                              cccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGt                                                                           <  2:1069371/118‑1 (MQ=255)                                                                               ccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtt                                                                                                                     >  2:1199173/1‑75 (MQ=255)                                                                               ccGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCtt                                                                                                                     <  1:1199173/75‑1 (MQ=255)                                                                                  cAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATAtt                                                  >  2:344565/1‑139 (MQ=255)                                                                                   aGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTg                                                 <  2:1442312/139‑1 (MQ=255)                                                                                     gTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCa                                               >  2:1158151/1‑139 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGTAGAAATAACATAAACTAAAAACACTTTTTGTTGATTGAGACAACCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATAGCTAAATATTGCACCTTTCtttt                                     >  1:806746/1‑139 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCtttt                                     >  1:442715/1‑139 (MQ=255)                                                                                               aGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCtttt                                     >  1:247041/1‑139 (MQ=255)                                                                                                    ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATTCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTccccc                                <  2:957386/139‑1 (MQ=255)                                                                                                    ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                                                     <  1:1272569/86‑1 (MQ=255)                                                                                                    ttGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                                                     >  2:1272569/1‑86 (MQ=255)                                                                                                     tGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt                               <  2:365680/139‑1 (MQ=255)                                                                                                     tGGAAAACAGTGTAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt                               <  2:788715/139‑1 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt                               <  2:407143/130‑1 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCt                               >  1:407143/1‑130 (MQ=255)                                                                                                              gtgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTaa                      >  1:1282633/1‑139 (MQ=255)                                                                                                               tgtAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCAATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAAc                     <  1:1158151/139‑1 (MQ=255)                                                                                                                 tAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                                                 >  1:1456719/1‑61 (MQ=255)                                                                                                                 tAAAAATTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCg                                                                                                 <  2:1456719/61‑1 (MQ=255)                                                                                                                  aaaaaTTACAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTa                  >  1:1029111/1‑139 (MQ=255)                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                                         >  1:1332895/1‑93 (MQ=255)                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCt                                                         <  2:1332895/93‑1 (MQ=255)                                                                                                                         aCAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTc           >  2:1428160/1‑139 (MQ=255)                                                                                                                          cAGAGAAGAAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCt          <  2:1029111/139‑1 (MQ=255)                                                                                                                             agaagaAATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTAt       >  2:863039/1‑139 (MQ=255)                                                                                                                               aagaaATAACAGAAACTGAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGccc                                                                                     <  2:990563/59‑2 (MQ=255)                                                                                                                               aagaaATAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCa                                                                                     >  1:990563/1‑59 (MQ=255)                                                                                                                                  aaaTAACAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCAtt  >  1:843897/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                      ACACACTTTGACGTGGGTCGCAAGGTGAGGCAAACCATTCAGGAACTTGGCGGAACCATGCCTGAGGAGTTGCCGACCCCGCAGGTCAGCATCAAGCAGTTGGAAAACAGTGTAAAAATTACAGAGAAGAAATAGCAGAAACTAAAAACTCTTTTTGTTGATTGAGACACCCGATGCGTAAGCCAAGGTCCAGGTGCAGTGAACACAATGGCTAAATATTGCACCTTTCTTTTCCCCCTCAGTTTTAACCTATTTTTTCTTATGCATT  >  NC_000913/3818450‑3818717

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