D66 :: BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,608,681	T→C	R59R (CGT→CGC)	SAUSA300_RS13420 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,608,681	0	T	C	99.1%	295.2 / ‑5.3	114	R59R (CGT→CGC)	SAUSA300_RS13420	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (46/67); total (47/67)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 115 total aligned reads displayed.

ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACA  >  NC_007793/2608585‑2608774

atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGc                                                                                            >  1:4039079/1‑100 (MQ=255)

atataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGc                                                                                            >  1:5081052/1‑100 (MQ=255)

  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtt                                                                                          >  1:4378216/1‑100 (MQ=255)

  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtt                                                                                          >  1:5687870/1‑100 (MQ=255)

  ataGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtt                                                                                          >  1:1559491/1‑100 (MQ=255)

   taGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCttt                                                                                         >  1:3029740/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCACTATAGGACGCGGCtttt                                                                                        >  1:4781435/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtttt                                                                                        >  1:90773/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtttt                                                                                        >  1:3362695/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtttt                                                                                        >  1:2019792/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCtttt                                                                                        >  1:1733710/1‑100 (MQ=255)

    aGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACACCCAATATAGGACGCGGCtttt                                                                                        >  1:821640/1‑100 (MQ=255)

     tcAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTa                                                                                       >  1:1745815/2‑100 (MQ=255)

     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTa                                                                                       >  1:1250004/1‑100 (MQ=255)

     gCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTa                                                                                       >  1:3269026/1‑100 (MQ=255)

      cAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAg                                                                                      >  1:871666/1‑100 (MQ=255)

       aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGt                                                                                     >  1:4830700/1‑100 (MQ=255)

       aGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGt                                                                                     >  1:4279882/1‑100 (MQ=255)

               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTcc                                                                             <  1:5136403/100‑1 (MQ=255)

               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTcc                                                                             <  1:1849938/100‑1 (MQ=255)

               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTcc                                                                             <  1:4164787/100‑1 (MQ=255)

               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTcc                                                                             <  1:3456405/100‑1 (MQ=255)

               tAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTcc                                                                             <  1:2864065/100‑1 (MQ=255)

                  aaTTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATACAGGACGCGGCTTTTAGTATCGTTCCAGa                                                                          <  1:648230/100‑1 (MQ=255)

                       tataTGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTtgt                                                                     <  1:1253386/100‑1 (MQ=255)

                       gctcttccgaTCTTGTAGATAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTtgt                                                                     <  1:5698716/91‑1 (MQ=255)

                         tcttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGGCTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgg                                                                   <  1:1144764/93‑1 (MQ=255)

                         tcttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTtgtgc                                                                   <  1:4855488/93‑2 (MQ=255)

                         tcttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTtgtgc                                                                   <  1:3094166/93‑2 (MQ=255)

                         tcttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgg                                                                   <  1:3361921/93‑1 (MQ=255)

                         tcttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgg                                                                   <  1:143455/93‑1 (MQ=255)

                         tatGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgg                                                                   <  1:1778409/100‑1 (MQ=255)

                          cttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggg                                                                  <  1:317473/94‑1 (MQ=255)

                          atGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggg                                                                  <  1:2440653/100‑1 (MQ=255)

                           ttccgattttGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:4586031/92‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTGTCGTTCCAGGTTTGTgggg                                                                 <  1:4815252/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:641448/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:2565364/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:3983022/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:4746687/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:5062290/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:5353585/95‑1 (MQ=255)

                           ttccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggga                                                                 <  1:4132139/95‑2 (MQ=255)

                           ttccgaTCTTGTAGATAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:3043106/95‑1 (MQ=255)

                           ttccgaTCTCGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:654133/95‑1 (MQ=255)

                           ttccgaTCTAGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggg                                                                 <  1:742923/95‑1 (MQ=255)

                           tGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggca                                                                 <  1:6762/100‑3 (MQ=255)

                            tccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTtggcacac                                                                <  1:2804450/96‑7 (MQ=255)

                            tccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggggt                                                                <  1:4949376/96‑2 (MQ=255)

                            tccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggggg                                                                <  1:301732/96‑2 (MQ=255)

                            tccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGa                                                                <  1:4674463/96‑1 (MQ=255)

                             tAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                               <  1:5265496/100‑1 (MQ=255)

                             ccgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAg                                                               <  1:5448244/97‑1 (MQ=255)

                              cgaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGt                                                              <  1:6508/98‑1 (MQ=255)

                               gaTCTTGTGGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTggggggag                                                             <  1:836187/99‑5 (MQ=255)

                               gaTCTTGTGGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                             <  1:3696849/99‑1 (MQ=255)

                               gaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCTAGATTTGTGGGGAGtt                                                             <  1:1616101/99‑1 (MQ=255)

                               gaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                             <  1:1147274/99‑1 (MQ=255)

                               gaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                             <  1:3018356/99‑1 (MQ=255)

                               gaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                             <  1:120428/99‑1 (MQ=255)

                               gaTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGtt                                                             <  1:5439942/99‑1 (MQ=255)

                                aTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                                            <  1:4932234/100‑1 (MQ=255)

                                aTCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTa                                                            <  1:4774889/100‑1 (MQ=255)

                                   ttGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAggg                                                         <  1:2837696/100‑1 (MQ=255)

                                    tggaGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:2256728/4‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATCATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATATGTGgagatcggaaga                                                        >  1:5553223/1‑92 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTgggagatcggaag                                                        >  1:1404396/1‑90 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:2851424/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:4421341/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:1345173/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:329868/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaa                                                        >  1:2800199/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                                        >  1:3782711/1‑100 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                                        >  1:1448763/1‑100 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTATgggagatcggaag                                                        >  1:5330204/1‑90 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATCTATGGGGAgatcgaag                                                        >  1:3594331/1‑93 (MQ=255)

                                    tGTAGGTAATCGTTTTAAAAATAATATAAGCATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGt                                                        >  1:3822792/1‑100 (MQ=255)

                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGTGGAgatcggaag                                                       >  1:272218/1‑92 (MQ=255)

                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaag                                                       >  1:3358214/1‑92 (MQ=255)

                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaag                                                       >  1:2747231/1‑92 (MQ=255)

                                     gTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGtt                                                       >  1:655697/1‑100 (MQ=255)

                                      tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                                      <  1:1111647/100‑1 (MQ=255)

                                      tAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                                      <  1:2610526/100‑1 (MQ=255)

                                      gaGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGttt                                                      <  1:476717/99‑1 (MQ=255)

                                       aGGTAATCGTTTTAAAAATAATCTAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAgatcggaagag                                                     >  1:5570350/1‑90 (MQ=255)

                                                       aataatATAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTc                                     <  1:1663568/100‑1 (MQ=255)

                                                        ataataTAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCt                                    <  1:3674595/100‑1 (MQ=255)

                                                          aataTAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTgt                                  <  1:4185778/100‑2 (MQ=255)

                                                          aataTAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                                  <  1:2679429/100‑1 (MQ=255)

                                                          aataTAAGTATGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGc                                  <  1:3384865/100‑1 (MQ=255)

                                                                   aTGACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTg                         >  1:2442878/1‑100 (MQ=255)

                                                                     gACTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCt                       <  1:3600515/100‑1 (MQ=255)

                                                                      aCTAAAGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTc                      <  1:5246998/100‑1 (MQ=255)

                                                                         aaaGCCACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCaga                   <  1:1694815/100‑3 (MQ=255)

                                                                             ccACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg               <  1:4442451/100‑1 (MQ=255)

                                                                             ccACATCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACg               <  1:2572247/100‑1 (MQ=255)

                                                                              cacaTCCAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGa              >  1:2457666/1‑100 (MQ=255)

                                                                                   ccAATATAGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCa         >  1:4147770/1‑100 (MQ=255)

                                                                                          aGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACa  >  1:4361779/1‑100 (MQ=255)

                                                                                          aGGACGCGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACa  >  1:554312/1‑100 (MQ=255)

ATATAGCAGACATGATAGAATTTTATATGTAAATCTTGTAGGTAATCGTTTTAAAAATAATATAAGTATGACTAAAGCCACATCCAATATAGGACGTGGCTTTTAGTATCGTTCCAGATTTGTGGGGAGTTAGGGTTTACTTAATTAAAGGGGTCTGCCCTCAATTGCTCACCAACGAGGGCAAGTTACA  >  NC_007793/2608585‑2608774

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