Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	973,711	Δ1 bp	pseudogene (196/373 nt)	SAUSA300_RS04780 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	973,709	0	A	.	100.0%	98.0 / NA	26	pseudogene (198/373 nt)	SAUSA300_RS04780	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (4/22);  total (4/22)

TACCTCTTGTGGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTTAAATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCTAATAAAAATAATTAATGAAACAT  >  NC_007793/973636‑973800                                                                          |                                                                                            tACCTCTTGTGGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCa                                                                  >  1:2117511/1‑100 (MQ=255)  aCCTCTTGTGGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCag                                                                 <  1:3847726/100‑2 (MQ=255)     tctTGTGGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTGCTAAAACCATAACTTGTCGCATCaa                                                              <  1:4497574/100‑1 (MQ=255)     tctTGTGGATTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCaa                                                              <  1:2770666/100‑1 (MQ=255)       gtgtgGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCaaaa                                                            >  1:6218103/2‑100 (MQ=255)              ttGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCctt                                                     <  1:2692044/100‑1 (MQ=255)                   aaaTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGACACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGt                                                <  1:6244437/100‑1 (MQ=255)                    aaTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCGTAACTTGTCGCATCAAAAATGCCTTCTTggg                                               <  1:5592943/100‑3 (MQ=255)                     aTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGGCATTTGATACTTT‑AATCTTCTAAATCCATAACTTGTCGCATCAAAAATGCCTTCTTGTAc                                              <  1:5139241/100‑1 (MQ=255)                     aTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATTCTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTAc                                              <  1:403739/100‑1 (MQ=255)                     aTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCTTCAAAAATGCCTTCTTGTAc                                              <  1:3404396/100‑1 (MQ=255)                     aTGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTAc                                              <  1:5301112/100‑1 (MQ=255)                      tGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACa                                             <  1:5213478/100‑1 (MQ=255)                              caGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAAt                                     >  1:6263207/1‑100 (MQ=255)                                 agaTAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCaa                                  <  1:2338059/100‑1 (MQ=255)                                          tcttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCt                         <  1:2166458/100‑1 (MQ=255)                                          tcttcatTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCt                         >  1:974947/1‑100 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATGACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCta                        <  1:3360236/100‑1 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAGTATGCta                        <  1:6189018/100‑1 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGcct                        <  1:337521/100‑3 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCta                        <  1:1541456/100‑1 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCta                        <  1:3539209/100‑1 (MQ=255)                                           cttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCta                        <  1:810489/100‑1 (MQ=255)                                            ttcttTTTCGAAGATGACATTTGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCtaa                       <  1:4284786/100‑1 (MQ=255)                                                               tttGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCTAATAAAAATAATTAATGAAAc    <  1:961474/100‑1 (MQ=255)                                                                 tGATACTTT‑AATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCTAATAAAAATAATTAATGAAACAt  <  1:1764749/100‑1 (MQ=255)                                                                          |                                                                                            TACCTCTTGTGGGTTGAAAAATGGATCATCAGAGATAGACTTCTTCTTTTTCGAAGATGACATTTGATACTTTAAATCTTCTAAAACCATAACTTGTCGCATCAAAAATGCCTTCTTGTACAAGTAAAATCAAAAATATGCTAATAAAAATAATTAATGAAACAT  >  NC_007793/973636‑973800

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