D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,721,335	T→C	L198S (TTA→TCA)	SAUSA300_RS13975 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,721,335	0	T	C	99.3%	409.3 / ‑4.8	150	L198S (TTA→TCA)	SAUSA300_RS13975	hydrolase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (139/10); minor base A (0/1); total (139/11)
Fisher's exact test for biased strand distribution p-value = 7.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 150 total aligned reads displayed.

ACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTC  >  NC_007793/2721238‑2721431

     aaaTATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAAttt                                                                                           >  1:4957662/1‑100 (MQ=255)

             cACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAg                                                                                   >  1:1866912/1‑100 (MQ=255)

                      aTGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATAc                                                                          >  1:3782085/1‑100 (MQ=255)

                            tcttccgatctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTc                                                                    <  1:5740082/92‑1 (MQ=255)

                            aaaTGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTc                                                                    <  1:2723413/100‑1 (MQ=255)

                            aaaTGATGTCTCAACCTGCTGAGACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTc                                                                    <  1:3765896/100‑1 (MQ=255)

                              atgatgTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCag                                                                  >  1:224124/1‑99 (MQ=255)

                               tgatgTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCaga                                                                 >  1:176049/1‑98 (MQ=255)

                                   ctctcAACCTGCTGATACTGAAGAAGGCCGTATAGAACAGTACAAACGTACGCTGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATAtt                                                             >  1:30746/2‑100 (MQ=255)

                                   cgctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATAtt                                                             >  1:4679254/3‑100 (MQ=255)

                                    tctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTa                                                            >  1:1938523/1‑100 (MQ=255)

                                     ctcAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTAc                                                           >  1:2153843/1‑100 (MQ=255)

                                      acAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAGACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACt                                                          >  1:2316431/2‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGGCGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:3266852/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGTACAGTAAAAACGTCCGCTGTCCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTCCta                                                         >  1:2674248/1‑99 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGGTGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:2959678/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTCGACAATATACGCATTCAGATATTACtt                                                         >  1:5657433/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGATAATATACGCATTCAGATATTACtt                                                         >  1:1632618/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:92502/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:4016014/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:4086071/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:1169524/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:2772910/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:5644572/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:2980535/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:5540886/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:4184433/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:1616474/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:5494976/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:3115637/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGCACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:4430458/1‑100 (MQ=255)

                                       cAACCTGCTGATACTGAAGAAGGACGAATAGAACAGTACAAACGTACGATGTTCTGGTCAGAGTTTGAAATTAGACAATATACGCATTCAGATATTACtt                                                         >  1:5298589/1‑100 (MQ=255)

                                        tgccTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACttt                                                        >  1:397258/3‑100 (MQ=255)

                                        aaCCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACttt                                                        >  1:543415/1‑100 (MQ=255)

                                         cccTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTg                                                       >  1:367961/2‑100 (MQ=255)

                                         agctgctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTg                                                       >  1:4932594/3‑100 (MQ=255)

                                         aCCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTg                                                       >  1:1588245/1‑100 (MQ=255)

                                          ccTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTgg                                                      >  1:450242/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:4953367/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:763624/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:5611702/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:3669041/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:5465187/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:4644471/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:1665984/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:4967841/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:2441400/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:2270280/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:4887602/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatg                                                   >  1:2779750/1‑100 (MQ=255)

                                             gctgATACTGAAGAAGGACGTATAGAACAGTACAAAAGCACGATGATCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTATGgatg                                                   >  1:3926693/1‑100 (MQ=255)

                                              ctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatga                                                  >  1:1907036/1‑100 (MQ=255)

                                              ctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatga                                                  >  1:724403/1‑100 (MQ=255)

                                              ctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatga                                                  >  1:1405771/1‑100 (MQ=255)

                                              ctgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatga                                                  >  1:4620996/1‑100 (MQ=255)

                                               tgATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatgat                                                 >  1:4521960/1‑100 (MQ=255)

                                               actgccTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGgatgat                                                 >  1:1659509/6‑100 (MQ=255)

                                                gATACTGACGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:1840885/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:5323388/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:47272/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:5400382/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:3524866/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:1090134/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:3726597/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:1772045/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:2546097/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:2913870/1‑100 (MQ=255)

                                                gATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAtt                                                >  1:2918224/1‑100 (MQ=255)

                                                 gacgggaagaaGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAttt                                               >  1:4296358/6‑100 (MQ=255)

                                                 aTACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGAttt                                               >  1:642203/1‑100 (MQ=255)

                                                  tACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTc                                              >  1:2047865/1‑100 (MQ=255)

                                                   ccTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTca                                             >  1:3533119/2‑100 (MQ=255)

                                                   aCTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTca                                             >  1:5385011/1‑100 (MQ=255)

                                                    cTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTcac                                            >  1:2523237/1‑100 (MQ=255)

                                                         ttagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAAt                                       >  1:3645521/3‑100 (MQ=255)

                                                         taagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAAt                                       >  1:106235/2‑100 (MQ=255)

                                                         caagGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGGATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAAt                                       >  1:123143/2‑100 (MQ=255)

                                                          aagGACGTATAGAGCAGTACAAACGTAAGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATa                                      >  1:3060828/1‑100 (MQ=255)

                                                           agGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATAc                                     >  1:1090717/1‑100 (MQ=255)

                                                           agGACGTATAGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATAc                                     >  1:1583398/1‑100 (MQ=255)

                                                               cctataGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCg                                 >  1:1322006/3‑100 (MQ=255)

                                                                 tataGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGAt                               <  1:1373206/100‑1 (MQ=255)

                                                                  ataGAACAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATa                              <  1:5076116/100‑1 (MQ=255)

                                                                      ctcAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAAt                          >  1:1309544/3‑100 (MQ=255)

                                                                      caCAGTACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAAt                          >  1:4531784/2‑100 (MQ=255)

                                                                          ctACAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCAcc                      >  1:4899241/2‑100 (MQ=255)

                                                                            aCAAACGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCct                    >  1:5724522/1‑100 (MQ=255)

                                                                               aaCGTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCctact                 >  1:1626707/1‑100 (MQ=255)

                                                                                 ttgaCGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaa               >  1:1616872/4‑100 (MQ=255)

                                                                                 tgTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaa               >  1:3824966/2‑100 (MQ=255)

                                                                                  gTACGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaaa              >  1:3318239/1‑100 (MQ=255)

                                                                                  cgaCGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaaa              >  1:2043664/3‑100 (MQ=255)

                                                                                    aCGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaaacg            >  1:4532716/1‑98 (MQ=255)

                                                                                    aCGATGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTaaacg            >  1:2638694/1‑98 (MQ=255)

                                                                                        tGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGAc        <  1:4044179/100‑1 (MQ=255)

                                                                                        tGTTCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAACGGGAc        <  1:2646570/100‑1 (MQ=255)

                                                                                          atCTGGTAAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAg      <  1:3140342/99‑1 (MQ=255)

                                                                                           tCTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGa     >  1:2019906/1‑100 (MQ=255)

                                                                                            cTGGTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGAt    >  1:1057169/1‑100 (MQ=255)

                                                                                              ggTCAGAATTTGAAATTAGACAATATACGCATTCAGATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAAAGGGACAGAtcg  >  1:3807151/1‑98 (MQ=255)

ACATTAAATATCGCACCAATTGATGCAAAAATGATGTCTCAACCTGCTGATACTGAAGAAGGACGTATAGAACAGTACAAACGTACGATGTTCTGGTTAGAATTTGAAATTAGACAATATACGCATTCAAATATTACTTTGGATGATTTCACAAAATACAGCGATAAAATCACCCTACTAAATGGGACAGATTC  >  NC_007793/2721238‑2721431

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: `‑`