Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	1,830,677	(A)6→7	pseudogene (137/929 nt)	tus →	DNA replication terminus site‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	1,830,671	1	.	A	100.0%	102.2 / NA	29	pseudogene (131/929 nt)	tus	DNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (12/17);  total (12/17)

GAGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑AAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAGGTAAA  >  NC_017635/1830625‑1830715                                                |                                             gAGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aa                                            <  2:1372806/49‑1 (MQ=255) gAGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aa                                            <  2:770653/49‑1 (MQ=255) gAGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑a                                             <  2:1198299/48‑1 (MQ=255)  aGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaa                                          <  2:1917949/50‑1 (MQ=255)     aaCACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaaaaa                                       <  1:1076652/50‑2 (MQ=255)     aaCACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaaa                                         >  1:1515644/1‑48 (MQ=255)     aaCACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaaa                                         >  2:75070/1‑48 (MQ=255)      acacAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaaaaag                                      <  2:1464521/50‑3 (MQ=255)      acacAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑aaaaaaag                                      <  2:1509924/50‑3 (MQ=255)            gCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATg                                 <  2:574683/50‑1 (MQ=39)              tATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAg                               >  2:1781663/1‑50 (MQ=39)              tATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAg                               >  1:901099/1‑50 (MQ=39)              tATTGGTTGCCCGCGTGTTCTCCTTGCCGGAGGTAAAAAAAGAGGATGAg                               >  2:131334/1‑50 (MQ=38)               aTTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAg                               >  1:625182/1‑49 (MQ=38)               aTTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAg                               >  2:103464/1‑49 (MQ=38)                ttGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGc                              <  1:1580728/49‑1 (MQ=38)                ttGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCa                             <  2:793468/50‑1 (MQ=39)                ttGGTTGCCCGCGTGTTCTCTTTGACGGAGGTAAAAAAAGAGGATGAGc                              <  2:1925036/49‑1 (MQ=38)                 tGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCa                             >  2:913643/1‑49 (MQ=38)                 tGGTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCa                             >  1:1267190/1‑49 (MQ=38)                  ggTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCAt                            >  2:1849280/1‑49 (MQ=38)                  ggTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCAt                            <  1:648211/49‑1 (MQ=38)                   gTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATaa                          <  2:103785/50‑1 (MQ=39)                   gTTGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATa                           <  2:1721596/49‑1 (MQ=38)                    ttGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAAt                         <  1:1439745/50‑1 (MQ=39)                    ttGCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAAt                         <  2:301336/50‑1 (MQ=39)                      gCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATcc                       >  2:33027/1‑50 (MQ=39)                      gCCCGCGTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATc                        <  2:712314/49‑1 (MQ=38)                          gcgTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCtt                   <  2:209645/50‑1 (MQ=39)                          gcgTGTTCTCTTCGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCtt                   <  1:1873425/50‑1 (MQ=38)                           cgTGTTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTa                  <  2:1569048/50‑1 (MQ=39)                            gtgtTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTaa                 >  1:730541/1‑50 (MQ=39)                              gtTCTCTTTGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATc               >  2:1538592/1‑50 (MQ=39)                                      tGCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTg         <  2:461065/48‑1 (MQ=38)                                       gCCGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAg       <  1:1985846/49‑1 (MQ=38)                                        ccGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAgg      >  2:1969943/1‑49 (MQ=38)                                        ccGGAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAGGt     <  1:1560368/50‑1 (MQ=39)                                           gAGGTAAAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAGGTaaa  <  1:7040/50‑1 (MQ=39)                                                |                                             GAGCAACACAAGCTATTGGTTGCCCGCGTGTTCTCTTTGCCGGAGGT‑AAAAAAGAGGATGAGCATAATCCGCTTAATCGTATTGAGGTAAA  >  NC_017635/1830625‑1830715

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