Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	781,031	A→G	Y557C (TAT→TGT)	SAUSA300_RS03785 →	heme ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	781,031	0	A	G	100.0%	273.8 / NA	91	Y557C (TAT→TGT)	SAUSA300_RS03785	heme ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (83/8);  total (83/8)

TATGAAAGTTATAAAAAATCATTAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATATAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTGC  >  NC_007793/780934‑781127                                                                                                  |                                                                                                  tATGAAAGTTATAAAAAATCATTAGATAAAAATAAATCCACATTGAGGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCAtgta                                                                                                 >  1:2569093/1‑100 (MQ=255)  aTGAAAGTTATAAAAAATCATTAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTaa                                                                                                >  1:3131794/1‑100 (MQ=255)         gtataAAAAATCATTAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGaaaaa                                                                                         >  1:4549887/2‑100 (MQ=255)                 aaTCATTAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAata                                                                                 >  1:4357303/1‑100 (MQ=255)                  tagggtAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAatat                                                                                >  1:2184704/6‑100 (MQ=255)                                   aatCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACg                                                               >  1:1929484/1‑100 (MQ=255)                                    ccacacaTTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGc                                                              >  1:3074182/4‑100 (MQ=255)                                    atCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGc                                                              >  1:344246/1‑100 (MQ=255)                                     accACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCa                                                             >  1:769371/2‑100 (MQ=255)                                     accACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCa                                                             >  1:1471497/2‑100 (MQ=255)                                     accACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCa                                                             >  1:4341802/2‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCag                                                            >  1:3591742/1‑99 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCAt                                                            >  1:1053785/1‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCAt                                                            >  1:1518378/1‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCAt                                                            >  1:5251168/1‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAAAGCAt                                                            >  1:2458292/1‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGCAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCAt                                                            >  1:1884471/1‑100 (MQ=255)                                      ccACATTGAAGCAACAATCTAAATCTTCTACAACTGAACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCAt                                                            >  1:4136681/1‑100 (MQ=255)                                              caGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑Ag                                                   >  1:625940/2‑100 (MQ=255)                                              caGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑Ag                                                   >  1:3284371/2‑100 (MQ=255)                                               aGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGc                                                  >  1:1437292/1‑100 (MQ=255)                                                 caacaaTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACG‑AGCgg                                                >  1:1259042/1‑100 (MQ=255)                                                  cacaaTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAGCGCACAGatcggaagag                                                >  1:2856957/2‑91 (MQ=255)                                                    gaaTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATGGAACA‑Agatcgga                                             >  1:1772995/2‑94 (MQ=255)                                                    caaTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGaga                                             <  1:4851016/100‑3 (MQ=255)                                                     aaTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGGACAATTGATGAAACGCATAGAACA‑AGCGGagta                                            <  1:4090139/100‑4 (MQ=255)                                                      aTCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTa                                           <  1:3355645/100‑1 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGGAGTAAGat                                       >  1:455782/1‑99 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:2409714/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:4039566/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:2406412/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:4343787/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:4847617/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGaa                                       >  1:5052129/1‑100 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGG‑AGTAAGatc                                      >  1:5166960/1‑98 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCACAGAACAGATCGGAAgagcac                                        >  1:546249/1‑95 (MQ=255)                                                          aaaTCTTCTACAACTGTACGTAAGAAAAATGGTTTATCACGAAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGat                                       >  1:4740291/1‑99 (MQ=255)                                                              tcttccgatcTGGACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGga                                   <  1:3711015/91‑1 (MQ=255)                                                              gctctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTagagatcgg                                   >  1:1257193/3‑92 (MQ=255)                                                               agcaaCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaa                                  <  1:2795451/96‑1 (MQ=255)                                                                tctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaag                                 >  1:4600578/1‑100 (MQ=255)                                                                tctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaag                                 >  1:4842183/1‑100 (MQ=255)                                                                tctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaag                                 >  1:597113/1‑100 (MQ=255)                                                                gctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaag                                 >  1:3694716/2‑100 (MQ=255)                                                                 ctACAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaaga                                >  1:1445160/1‑100 (MQ=255)                                                                   aCAACTGTACGTAGGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGgaagatc                              >  1:1088546/1‑98 (MQ=255)                                                                   aCAACTGTACGTAAGAAATATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:316361/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGGAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:4180992/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAGGTAAGATCGgaagagc                              >  1:4631358/1‑98 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGATCGgaagagc                              >  1:5478637/1‑98 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:4037385/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:4511980/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:966387/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3922939/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3170270/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:2892868/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:2899398/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:2953953/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3018610/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3103991/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3114242/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:2104150/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3130093/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:1821518/1‑100 (MQ=255)                                                                   aCAACTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGaaa                              >  1:3783898/1‑100 (MQ=255)                                                                     aaCTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAAtt                            >  1:1880053/1‑100 (MQ=255)                                                                      ccTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTg                           >  1:3170589/2‑100 (MQ=255)                                                                      aCTGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑TGCGGAAGTAAGAATGGAAGAAATTg                           >  1:1522190/1‑100 (MQ=255)                                                                       atGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGa                          >  1:3396039/2‑100 (MQ=255)                                                                       atGTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGa                          >  1:1433258/2‑100 (MQ=255)                                                                         gTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtg                        >  1:2401521/1‑100 (MQ=255)                                                                         gTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtg                        >  1:1914538/1‑100 (MQ=255)                                                                         gTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtg                        >  1:859995/1‑100 (MQ=255)                                                                         gTACGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtg                        >  1:4462459/1‑100 (MQ=255)                                                                          gaCGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtgt                       >  1:3221680/2‑100 (MQ=255)                                                                           tcGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtgtg                      >  1:342324/2‑100 (MQ=255)                                                                           aCGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtgtg                      >  1:5393922/1‑100 (MQ=255)                                                                           aCGTAAGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGAtgtg                      >  1:4932816/1‑100 (MQ=255)                                                                               aaGAAAAATGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCa                  <  1:5385458/100‑1 (MQ=255)                                                                                    ccaTGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGAtt             >  1:4407739/3‑100 (MQ=255)                                                                                      aTGGTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGa           >  1:1988821/1‑100 (MQ=255)                                                                                        ggTTTATCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACACATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAgg         >  1:142585/1‑100 (MQ=255)                                                                                             ctCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGt    >  1:3246845/2‑100 (MQ=255)                                                                                             aTCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCGTGATTGAGGCAAGt    >  1:3936545/1‑100 (MQ=255)                                                                                             aTCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGt    >  1:569717/1‑100 (MQ=255)                                                                                             aTCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGt    >  1:1399320/1‑100 (MQ=255)                                                                                             aTCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTAAGGCAAga    >  1:5448663/1‑99 (MQ=255)                                                                                             aTCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGAAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGt    >  1:5175293/1‑100 (MQ=255)                                                                                              tCATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTg   <  1:387837/100‑1 (MQ=255)                                                                                               cATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTGc  >  1:1661893/1‑100 (MQ=255)                                                                                               cATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTGc  <  1:3416331/100‑1 (MQ=255)                                                                                               cATGTAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGGAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTGc  >  1:2039405/1‑100 (MQ=255)                                                                                                  |                                                                                                  TATGAAAGTTATAAAAAATCATTAGATAAAAATAAATCCACATTGAAGCAACAATCTAAATCTTCTACAACTGTACGTAAGAAAAATGGTTTATCATATAAAGAAAAATTAGAATATGAACAATTGATGAAACGCATAGAACA‑AGCGGAAGTAAGAATGGAAGAAATTGATGTGCTCATGATTGAGGCAAGTGC  >  NC_007793/780934‑781127

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