Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	525,622	(A)8→7	pseudogene (63/157 nt)	SAUSA300_RS02505 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	525,615	0	A	.	100.0%	181.9 / NA	46	pseudogene (70/157 nt)	SAUSA300_RS02505	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (46/0);  total (46/0)

AAACGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGATAAAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACGA  >  NC_007793/525530‑525710                                                                                      |                                                                                                aaaCGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATCTTGACTATGAT‑AAAAAAATGATGACa                                                                                  >  1:4580110/1‑100 (MQ=255)             taaCGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGtttttt                                                                      >  1:1767375/2‑100 (MQ=255)             aaaCGCAAGCTGTAATTTTAAATGGTGAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGtttttt                                                                      >  1:4120114/1‑100 (MQ=255)              aaCGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTa                                                                     >  1:3570617/1‑100 (MQ=255)                 gCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATg                                                                  >  1:1941868/1‑100 (MQ=255)                  cAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGt                                                                 >  1:2643171/1‑100 (MQ=255)                   taGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTa                                                                >  1:5386462/2‑100 (MQ=255)                   caGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTa                                                                >  1:723355/2‑100 (MQ=255)                   aaGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTa                                                                >  1:4470675/1‑100 (MQ=255)                    aGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTaa                                                               >  1:4306901/1‑100 (MQ=255)                     tcTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTaaa                                                              >  1:1217118/2‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:2246458/1‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:994989/1‑100 (MQ=255)                                aaaTGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACAtt                                                   >  1:794571/1‑100 (MQ=255)                                               cttttGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATAc                                    >  1:4414272/2‑100 (MQ=255)                                                 cttGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCt                                  >  1:1687605/2‑100 (MQ=255)                                                  ttGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCtt                                 >  1:2956271/1‑100 (MQ=255)                                                   tGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTa                                >  1:5238754/1‑100 (MQ=255)                                                    gCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTaa                               >  1:443995/1‑100 (MQ=255)                                                        tGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAGGTAATAAATACCTTAACAtt                           >  1:4516637/1‑100 (MQ=255)                                                        tGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACAtt                           >  1:4735507/1‑100 (MQ=255)                                                        cgATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACAtt                           >  1:86547/2‑100 (MQ=255)                                                         taTACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTg                          >  1:4666367/2‑100 (MQ=255)                                                         gATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTg                          >  1:5466500/1‑100 (MQ=255)                                                         gATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTg                          >  1:2654230/1‑100 (MQ=255)                                                            aCTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAAt                       >  1:1323348/1‑100 (MQ=255)                                                             ttGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATa                      >  1:1692879/2‑100 (MQ=255)                                                               gATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAg                    >  1:4132174/1‑100 (MQ=255)                                                                      tgatgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGaaaa             >  1:1440482/1‑100 (MQ=255)                                                                      tgatgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGaaaa             >  1:2022179/1‑100 (MQ=255)                                                                       aatgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAAt            >  1:5310544/2‑100 (MQ=255)                                                                         tgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATga          >  1:4040026/1‑100 (MQ=255)                                                                           ttGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATG‑AAATGAGAt       >  1:5495107/1‑100 (MQ=255)                                                                           ctGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATgaga        >  1:2671597/2‑100 (MQ=255)                                                                            tGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGAt       >  1:2216062/1‑100 (MQ=255)                                                                            tGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGAt       >  1:2021876/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATGAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:3333383/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGGTg      >  1:3920223/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGGTg      >  1:4841204/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:338761/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:2797542/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:5524851/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:1941809/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAACgagaca      >  1:4030572/1‑98 (MQ=255)                                                                                 aTGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAGTGagaagatcg  >  1:3082105/1‑97 (MQ=255)                                                                                 aTGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACAGTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGGTGACga  >  1:4225394/1‑100 (MQ=255)                                                                                      |                                                                                                AAACGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGATAAAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACGA  >  NC_007793/525530‑525710

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑