Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	621,998	T→G	S1035A (TCA→GCA)	SAUSA300_RS02925 →	MSCRAMM family adhesin SdrE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	621,998	0	T	G	81.5%	10.3 / 39.0	94	S1035A (TCA→GCA)	SAUSA300_RS02925	MSCRAMM family adhesin SdrE
Reads supporting (aligned to +/- strand):  ref base T (16/1);  new base G (73/2);  total (91/3)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.21e-01

CAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGACTCAGATA  >  NC_007793/621903‑622094                                                                                                |                                                                                                 cagattcagatAGCGATTCGGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGACAGCGACTCGGACTCAGa                                                                                              >  1:771374/1‑100 (MQ=255) cagacgcagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGa                                                                                              >  1:1970806/7‑100 (MQ=0)  cgacgcagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGAt                                                                                             >  1:4356651/6‑100 (MQ=2)  agacgcagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGAt                                                                                             >  1:2708598/6‑100 (MQ=2)     cgcagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGc                                                                                          >  1:3236714/3‑100 (MQ=11)     cgcagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGc                                                                                          >  1:5633087/3‑100 (MQ=11)       ctgatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGa                                                                                        >  1:1762285/3‑100 (MQ=11)       cagatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGa                                                                                        >  1:4370790/1‑100 (MQ=11)       cagaaAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCGGACGCAGATAGCGa                                                                                        >  1:4786048/1‑100 (MQ=11)        agatAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGTGAc                                                                                       >  1:5549173/1‑100 (MQ=11)             tcaccagtgACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGAATCAGATTCAGATAGAGATTCGGACCCAGATAGCGACTCAGa                                                                                  >  1:4479657/9‑100 (MQ=17)                    tgatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGAc                                                                           >  1:3318293/8‑100 (MQ=11)                    tgatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGAc                                                                           >  1:5333587/8‑100 (MQ=11)                    tgatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGAc                                                                           >  1:1946525/8‑100 (MQ=11)                    agatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGAc                                                                           >  1:5376545/8‑100 (MQ=11)                    acgtgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGAc                                                                           >  1:371752/8‑100 (MQ=11)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACa                                                                          >  1:2538671/7‑100 (MQ=2)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACa                                                                          >  1:3108364/7‑100 (MQ=11)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACa                                                                          >  1:5058041/7‑100 (MQ=11)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACa                                                                          >  1:3415496/7‑100 (MQ=11)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACa                                                                          >  1:2446083/7‑100 (MQ=11)                     gatgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACa                                                                          >  1:2668679/7‑100 (MQ=11)                      ctgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAg                                                                         >  1:5866427/6‑100 (MQ=1)                      atgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATGGCGACTCTGATGCGGACAg                                                                         >  1:6500706/6‑100 (MQ=0)                       tgcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGACAGCGACTCTGATGCGGACAGc                                                                        >  1:588482/5‑100 (MQ=11)                        gcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAgtg                                                                       >  1:1941969/4‑98 (MQ=1)                        gcggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAACg                                                                       >  1:5306453/4‑100 (MQ=11)                        ctggACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCg                                                                       >  1:228122/4‑100 (MQ=11)                         tgtgcAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGa                                                                      >  1:1666763/5‑100 (MQ=11)                         cAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGGCAGCGa                                                                      >  1:1444034/1‑100 (MQ=255)                         cAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGa                                                                      >  1:6416007/1‑100 (MQ=255)                         cAGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGa                                                                      >  1:731881/1‑100 (MQ=11)                          tgACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGac                                                                     >  1:4451957/2‑99 (MQ=11)                          gcccAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGac                                                                     >  1:5386500/4‑99 (MQ=11)                          acaCAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAgatcg                                                                     >  1:2927435/3‑96 (MQ=1)                          acaCAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAgatcg                                                                     >  1:3733524/3‑96 (MQ=1)                          aagcAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGac                                                                     >  1:2264759/4‑99 (MQ=11)                          aGACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGac                                                                     >  1:353538/1‑99 (MQ=11)                           taCAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGact                                                                    >  1:4649390/2‑98 (MQ=11)                           gACAGCGATTCGGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGact                                                                    >  1:2226404/1‑98 (MQ=11)                           gACAGCGATTCAGACGCAGATAGCGATTCAGATTCAGACAGCGACTCAGACGCAGATAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAgatcgg                                                                    >  1:3892437/1‑95 (MQ=1)                            ataGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTc                                                                   >  1:5049745/3‑100 (MQ=2)                            aCAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCAATTc                                                                   <  1:2498155/100‑1 (MQ=255)                              aGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAg                                                                 >  1:4757380/1‑100 (MQ=11)                              aGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAg                                                                 >  1:2241666/1‑100 (MQ=11)                               gCGATTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGa                                                                >  1:1072787/1‑100 (MQ=11)                                actaccAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAc                                                               >  1:479611/6‑100 (MQ=2)                                  aTTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcgc                                                             <  1:4093712/100‑3 (MQ=2)                                  aTTCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGAcgc                                                             <  1:1487784/100‑3 (MQ=2)                                   ttCAGACGCGGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCa                                                            >  1:3891958/1‑100 (MQ=2)                                    tCAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAg                                                           >  1:6218680/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGCTTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:6553633/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGCTTCAGACGCAGa                                                          >  1:2006529/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:5578876/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:5729543/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:6079432/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:6626753/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:4773021/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:4316075/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:1058672/1‑100 (MQ=2)                                     cAGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGa                                                          >  1:3557434/1‑100 (MQ=2)                                      cgACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGAt                                                         >  1:5559445/2‑100 (MQ=2)                                      aGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGAt                                                         >  1:6063745/1‑100 (MQ=2)                                      aGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGAt                                                         >  1:5467323/1‑100 (MQ=2)                                      aGACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCAATTCAGACGCAGAt                                                         >  1:4887830/1‑100 (MQ=2)                                       gACGCAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATa                                                        >  1:4232293/1‑100 (MQ=2)                                       gACGCAGATAGCGACTCTGATGCGGACAGAGACTAAGATTCAGACAGAGATTCTGACGCCGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGaaa                                                        >  1:4470891/1‑98 (MQ=2)                                        acgcAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAg                                                       >  1:2127051/4‑100 (MQ=2)                                         tgcAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGTTGCGGACAGCGATTCAGACGCAGATAGc                                                      >  1:2472605/3‑100 (MQ=2)                                         gccAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGc                                                      >  1:4436071/3‑100 (MQ=2)                                         cTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGc                                                      >  1:4847407/1‑100 (MQ=255)                                          gcAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCg                                                     >  1:3510038/2‑100 (MQ=11)                                          gcAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCg                                                     >  1:435516/2‑100 (MQ=11)                                          gcAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCg                                                     >  1:6487560/2‑100 (MQ=11)                                          ccAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCg                                                     >  1:2854025/2‑100 (MQ=11)                                           cAGATAGCGACTCTGCTGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCga                                                    >  1:3691951/1‑100 (MQ=2)                                           cAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCga                                                    >  1:4225112/1‑100 (MQ=11)                                           cAGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCga                                                    >  1:1768578/1‑100 (MQ=11)                                           cAGATAGCGACTCAGACTCAGACAGCGACTCAGCTTCAGATAGCGATTCGGACACAGATAGCGAATAAGATTCAGACAGCGATTCAGACTCAGATAGCga                                                    >  1:4573792/1‑100 (MQ=14)                                           cAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCga                                                    >  1:3173489/1‑100 (MQ=255)                                           cAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCga                                                    >  1:879034/1‑100 (MQ=255)                                            aGATAGCGACTCTGATGCGGACAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCGGACGCAGAtcggaag                                                   >  1:1422355/1‑94 (MQ=2)                                                           agcggacAGCGACTCAGATTCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCGGACAGCGATTCAGACGCAGATAGCGACTCTGATGCGGACAgc                                    >  1:3906094/5‑99 (MQ=2)                                                             cagacAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTga                                  >  1:5374881/1‑100 (MQ=255)                                                             cagacAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTga                                  >  1:1952207/1‑100 (MQ=255)                                                             cagacAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTga                                  >  1:1547360/1‑100 (MQ=255)                                                             cagacAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTga                                  >  1:1238746/1‑100 (MQ=255)                                                             cagacAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCga                                  >  1:1867806/1‑100 (MQ=255)                                                                         cagattcagatAGCGATTCGGACTCACATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCGGACTCAGACAGTGAttcagattcaga                      >  1:4714165/1‑100 (MQ=255)                                                                                agatAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCgac               >  1:5128846/1‑100 (MQ=255)                                                                                agatAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCgac               >  1:836328/1‑100 (MQ=255)                                                                                    aGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCgactcag           >  1:2168418/1‑100 (MQ=255)                                                                                           cGGACGCAGATAGCGATTCAGATGCGGACAGCGATTCAGACGCAGACAGCGATTCAGACGCAGATAGCGATTCTGATGCAGACAGCGActcagacgcaga    >  1:5110261/1‑100 (MQ=2)                                                                                             gACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGACTCAGATa  >  1:6011579/1‑100 (MQ=40)                                                                                                |                                                                                                 CAGATTCAGATAGCGATTCGGACTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGATAGCGACTCAGATTCAGACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGACTCAGATA  >  NC_007793/621903‑622094

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