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%	141.8 / NA	36	pseudogene (70/157 nt)	SAUSA300_RS02505	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (36/0);  total (36/0)

AAACGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGATAAAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACGA  >  NC_007793/525530‑525710                                                                                      |                                                                                                aaaCGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACa                                                                                  >  1:389628/1‑100 (MQ=255)          tataAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGttt                                                                         >  1:4947792/1‑100 (MQ=255)          tataAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGttt                                                                         >  1:2951929/1‑100 (MQ=255)          caccgaCGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGttt                                                                         >  1:3356175/6‑100 (MQ=255)             aaaCGCAAGCTGTAATTTTAAATGGAAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGtttttt                                                                      >  1:5938929/1‑100 (MQ=255)                 gCAAGCTGTAATTTTAAATGGTAAGTTATCGTTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATg                                                                  >  1:4813613/1‑100 (MQ=255)                 gCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATg                                                                  >  1:16534/1‑100 (MQ=255)                   gcccTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTGTAATGTa                                                                >  1:4610075/4‑100 (MQ=255)                   aaGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTa                                                                >  1:4648547/1‑100 (MQ=255)                   aaGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTa                                                                >  1:2545760/1‑100 (MQ=255)                    aGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTaa                                                               >  1:4387560/1‑100 (MQ=255)                      gtGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTaaaa                                                             >  1:5206412/2‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:2800622/1‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:4563652/1‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:3589485/1‑100 (MQ=255)                        gTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAaaata                                                           >  1:2245271/1‑100 (MQ=255)                                 aaTGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTg                                                  >  1:5608026/1‑100 (MQ=255)                                    cttAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGaaa                                               >  1:629022/3‑100 (MQ=255)                                                ttttGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATAcc                                   >  1:2091667/1‑100 (MQ=255)                                                  ttGCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCtt                                 >  1:600261/1‑100 (MQ=255)                                                  ttGCATTGATACTGATAAAATGATGTTGACTATGAT‑‑AAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTa                                >  1:1127653/1‑100 (MQ=255)                                                    gCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTaa                               >  1:3972336/1‑100 (MQ=255)                                                    gCATTGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTaa                               >  1:3923179/1‑100 (MQ=255)                                                       ttGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACAt                            >  1:1769467/1‑100 (MQ=255)                                                        tGATACTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACAtt                           >  1:4232010/1‑100 (MQ=255)                                                            aCTGATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAAt                       >  1:3407381/1‑100 (MQ=255)                                                               gATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAg                    >  1:687192/1‑100 (MQ=255)                                                               gATAAAATGATGTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAg                    >  1:1225082/1‑100 (MQ=255)                                                                     ttgatgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGaaa              >  1:146386/2‑100 (MQ=255)                                                                         tgTTGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATga          >  1:2846472/1‑100 (MQ=255)                                                                            tGACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGAt       >  1:4641320/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:5470800/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATg      >  1:2828461/1‑100 (MQ=255)                                                                             gACTATGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGAT‑TAAATGAGAtcg     >  1:5491168/1‑98 (MQ=255)                                                                                 aTGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACga  >  1:542634/1‑100 (MQ=255)                                                                                 aTGAT‑AAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACga  >  1:879403/1‑100 (MQ=255)                                                                                      |                                                                                                AAACGATTTTATAAACGCAAGCTGTAATTTTAAATGGTAAGTTATCATTTTGCATTGATACTGATAAAATGATGTTGACTATGATAAAAAAAATGATGACATCGACGTTTTTTAATGTAAAATAAATACATTGAAAGTAATAAATACCTTAACATTGAATAAGATGAAAATGAGATGACGA  >  NC_007793/525530‑525710

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