Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	1,201,459	T→A	intergenic (+7/+34)	dgcT → / ← ECW_RS05820	GGDEF domain‑containing protein/IS3‑like element IS3 family transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	1,201,459	0	T	A	100.0%	145.4 / NA	40	intergenic (+7/+34)	dgcT/ECW_RS05820	GGDEF domain‑containing protein/IS3‑like element IS3 family transposase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (20/20);  total (20/20)

TGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCTTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTT  >  NC_017635/1201412‑1201504                                                |                                              tggtTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCatt                                             <  2:250372/50‑1 (MQ=255)   gtTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCatttt                                           >  2:1814167/1‑50 (MQ=255)     tgtgAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGa                                         >  2:626096/1‑50 (MQ=255)        gAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATcc                                      >  1:1369806/1‑50 (MQ=255)         aGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATcc                                      <  2:391797/49‑1 (MQ=255)         aGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCt                                     <  2:893373/50‑1 (MQ=255)         aGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCt                                     <  2:1710679/50‑1 (MQ=255)         aGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCt                                     <  1:79852/50‑1 (MQ=255)          gggATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTa                                    >  2:1478508/1‑50 (MQ=255)            gATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTAcc                                  >  2:654484/1‑50 (MQ=255)            gATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTAcc                                  >  1:2063392/1‑50 (MQ=255)            gATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTAcc                                  <  2:1343984/50‑1 (MQ=255)            gATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTAcc                                  >  1:1030108/1‑50 (MQ=255)              tgtgGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCa                                >  2:378823/1‑50 (MQ=255)                  gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGt                             >  2:1644168/1‑49 (MQ=255)                  gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTa                            <  2:1689593/50‑1 (MQ=255)                     aaTTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAata                         >  1:723473/1‑50 (MQ=255)                     aaTTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAata                         <  2:1509105/50‑1 (MQ=255)                      aTTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAatat                        >  1:1759203/1‑50 (MQ=255)                      aTTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAatat                        <  2:1868823/50‑1 (MQ=255)                       tttttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATg                       >  2:88246/1‑50 (MQ=255)                       tttttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATg                       >  1:55017/1‑50 (MQ=255)                         tttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATgg                      <  1:1811405/49‑1 (MQ=255)                          ttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGa                     >  1:39351/1‑49 (MQ=255)                           tgtgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGac                    <  2:418500/49‑1 (MQ=255)                           tgtgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGac                    <  2:422339/49‑1 (MQ=255)                             tgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGacac                  <  1:2004041/49‑1 (MQ=255)                             tgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGAcaca                 <  2:1552956/50‑1 (MQ=255)                             tgAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGAcaca                 <  2:852371/50‑1 (MQ=255)                               aGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAg                <  1:706887/49‑1 (MQ=255)                               aGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAg                <  1:1748562/49‑1 (MQ=255)                                  cACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGcc             >  1:65693/1‑49 (MQ=255)                                    ccATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCt            >  1:283423/1‑48 (MQ=255)                                    ccATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTaa          <  2:1152221/50‑1 (MQ=255)                                     cATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAg         >  1:272143/1‑50 (MQ=255)                                     cATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAg         <  1:472313/50‑1 (MQ=255)                                       tAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCg       <  2:153484/50‑1 (MQ=255)                                         aaGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAg     >  1:727714/1‑50 (MQ=255)                                          aGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAgg    >  1:1554264/1‑50 (MQ=255)                                          aGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAgg    >  2:1588762/1‑50 (MQ=255)                                            cGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGtt  <  1:898629/50‑1 (MQ=255)                                            cGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGtt  >  2:950879/1‑50 (MQ=255)                                                |                                              TGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCTTTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTT  >  NC_017635/1201412‑1201504

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