Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	529,026	A→C	E98D (GAA→GAC)	SAUSA300_RS02520 →	ribonuclease M5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	529,026	0	A	C	100.0%	157.1 / NA	60	E98D (GAA→GAC)	SAUSA300_RS02520	ribonuclease M5
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (53/7);  total (53/7)

ATTACTGAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGAATTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTA  >  NC_007793/528934‑529122                                                                                             |                                                                                                 aTTACTGAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAGGGGAAAATTGGTGTTGAACATGCCGACTTAATTg                                                                                           <  1:1875673/100‑1 (MQ=255)   gtgtGAACAAGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGat                                                                                         >  1:5561384/4‑100 (MQ=255)       gAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAGAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTa                                                                                     >  1:37187/1‑100 (MQ=255)       gAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTa                                                                                     >  1:928189/1‑100 (MQ=255)       gAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTa                                                                                     >  1:458946/1‑100 (MQ=255)        gaCATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTaa                                                                                    >  1:3595105/2‑100 (MQ=255)         aCATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTaaa                                                                                   >  1:3395515/1‑100 (MQ=255)            tGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAaagaa                                                                                >  1:1691183/1‑100 (MQ=255)               cAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCg                                                                             >  1:4123714/1‑100 (MQ=255)                 taGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGtt                                                                           >  1:4963264/2‑100 (MQ=255)                 aaGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGtt                                                                           >  1:346349/1‑100 (MQ=255)                  aGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGttg                                                                          <  1:2728993/100‑2 (MQ=255)                          aaCATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGt                                                                  >  1:1149618/1‑100 (MQ=255)                            cATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTa                                                                >  1:2529190/1‑100 (MQ=255)                            aaTGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTa                                                                >  1:3390016/2‑100 (MQ=255)                             gtGCGTATATTGATAGAGAAAAAGCTAACAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAg                                                               >  1:4438072/2‑100 (MQ=255)                             ctGCGTATATTGATAGAGAAAAAGCTAAAAATAAACAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAg                                                               >  1:604869/2‑100 (MQ=255)                              tGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGt                                                              >  1:4541102/1‑100 (MQ=255)                               gCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTCAATGCATGTTAGtt                                                             >  1:4671751/1‑100 (MQ=255)                                         aTAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTtga                                                   <  1:1660475/100‑1 (MQ=255)                                         aTAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTtga                                                   >  1:5529437/1‑100 (MQ=255)                                                 ccaGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCtt                                           >  1:1595247/3‑100 (MQ=255)                                                   aGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTAt                                         >  1:4839462/1‑100 (MQ=255)                                                   aGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCATTTGATGAAGCTTAt                                         >  1:2265664/1‑100 (MQ=255)                                                    ccTAAAAATAAAAGAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATg                                        >  1:4832576/2‑100 (MQ=255)                                                       aaaaataaaaaAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAAt                                     >  1:2751333/1‑100 (MQ=255)                                                       aaaaataaaaaAGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAAt                                     >  1:4297541/1‑100 (MQ=255)                                                              ttgcaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGa                              >  1:3163743/5‑100 (MQ=255)                                                                ccaGGGAAAATTGGTGTTGAACACGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATa                            >  1:2993927/3‑100 (MQ=255)                                                                aaaGGGAAAATTGGTGTCGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATa                            >  1:817701/1‑100 (MQ=255)                                                                 caGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATaa                           >  1:2720479/2‑100 (MQ=255)                                                                 aaGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATaa                           >  1:2352580/1‑100 (MQ=255)                                                                  aGGGAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATaaa                          >  1:3734231/1‑100 (MQ=255)                                                                   gggAAAATTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAAt                         >  1:3850205/1‑100 (MQ=255)                                                                      aaaaTTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCtg                      >  1:3940525/1‑100 (MQ=255)                                                                         cttGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGc                   >  1:3003326/2‑100 (MQ=255)                                                                         aTTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCtgagg                   <  1:4165752/100‑4 (MQ=255)                                                                         aTTGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATCGATAAATCTGTGc                   >  1:4055328/1‑100 (MQ=255)                                                                          ttGGTGTTGAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCt                  >  1:3869297/1‑100 (MQ=255)                                                                                  gAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGt          >  1:4076261/1‑100 (MQ=255)                                                                                  gAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGt          >  1:3903308/1‑100 (MQ=255)                                                                                  gAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGt          >  1:3093127/1‑100 (MQ=255)                                                                                  gAACATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGt          >  1:2503122/1‑100 (MQ=255)                                                                                    aCATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATTGAGTTg        >  1:5631048/1‑100 (MQ=255)                                                                                     gaTGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTgg       >  1:2207974/2‑100 (MQ=255)                                                                                     cATGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTgg       >  1:2319121/1‑100 (MQ=255)                                                                                      gtGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTggg      >  1:3699439/2‑100 (MQ=255)                                                                                      gtGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTggg      >  1:3229595/2‑100 (MQ=255)                                                                                      aTGCCGACTTAATTGATATTAAAGAAGCGTTAATGCGTGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTggg      >  1:3471951/1‑100 (MQ=255)                                                                                      aTGCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTggg      >  1:233313/1‑100 (MQ=255)                                                                                        gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt    >  1:4592721/1‑100 (MQ=255)                                                                                        gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt    >  1:4848463/1‑100 (MQ=255)                                                                                        gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt    >  1:3043617/1‑100 (MQ=255)                                                                                        gCCGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt    >  1:5674836/1‑100 (MQ=255)                                                                                        accGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGt    >  1:2382270/2‑100 (MQ=255)                                                                                         gcGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGtt   >  1:2720562/2‑100 (MQ=255)                                                                                         ccGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGtt   <  1:5342371/100‑1 (MQ=255)                                                                                          cGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTa  <  1:2524440/100‑1 (MQ=255)                                                                                          cGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTa  >  1:4739124/1‑100 (MQ=255)                                                                                          cGACTTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTa  <  1:5715946/100‑1 (MQ=255)                                                                                             |                                                                                                 ATTACTGAACATGTCAAAGGTGTTAAACATGCGTATATTGATAGAGAAAAAGCTAAAAATAAAAAAGGGAAAATTGGTGTTGAACATGCCGAATTAATTGATATTAAAGAAGCGTTAATGCATGTTAGTTCACCCTTTGATGAAGCTTATGAATCAATTGATAAATCTGTGCTAATAGAGTTGGGGTTA  >  NC_007793/528934‑529122

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