BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,520	A→G	S531P (TCC→CCC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,520	0	A	G	97.6%	386.2 / ‑1.1	124	S531P (TCC→CCC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/1); new base G (51/69); total (54/70)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 130 total aligned reads displayed.

ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609

aTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGgggg                                                                                         >  1:689587/1‑93 (MQ=255)

aTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGgggg                                                                                         >  2:3622675/1‑93 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       >  2:86692/1‑93 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       >  1:1958738/1‑93 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       >  2:1394291/1‑93 (MQ=255)

  acacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCa                                                                                       >  1:2425545/1‑93 (MQ=255)

   cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAt                                                                                      <  2:79172/93‑1 (MQ=255)

    acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAtt                                                                                     <  2:713652/93‑1 (MQ=255)

     ccGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCAttt                                                                                    <  1:4428107/93‑1 (MQ=255)

      cGACCTGTAGCAATTTTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtc                                                                                   <  2:1292432/93‑1 (MQ=255)

      cGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtc                                                                                   <  1:4147417/93‑1 (MQ=255)

          cTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTtcctcc                                                                               <  1:1337049/93‑1 (MQ=255)

                 aaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAAt                                                                        >  2:115252/1‑93 (MQ=255)

                 aaTTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAAt                                                                        >  2:3996049/1‑93 (MQ=255)

                   ttGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATcc                                                                      >  1:4714373/1‑93 (MQ=255)

                    tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCa                                                                     <  2:4010023/93‑1 (MQ=255)

                    tgtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCa                                                                     <  1:1968998/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  2:583315/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  2:850602/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  1:4419034/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  1:3388554/93‑1 (MQ=255)

                      tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATa                                                                   <  2:1137706/93‑1 (MQ=255)

                       gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAg                                                                  <  1:667160/93‑1 (MQ=255)

                       gCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAg                                                                  >  2:777653/1‑93 (MQ=255)

                         gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGcc                                                                >  1:3139205/1‑93 (MQ=255)

                         gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGcc                                                                <  2:1129327/93‑1 (MQ=255)

                          ttAGTAGAGGATGTCTTTAACTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTCCCGGGGGCATTTCCTCCGATAAATCCATAGccc                                                               <  2:48189/93‑1 (MQ=255)

                           tagtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTGCCGGGGGCATTTCCTCCGATAAATCCATAGCCCg                                                              <  1:758762/93‑1 (MQ=255)

                             gtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGtt                                                            >  2:131605/1‑93 (MQ=255)

                                gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         <  1:4099115/93‑1 (MQ=255)

                                gatgatGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGaa                                                         >  1:3355112/1‑93 (MQ=255)

                                 atgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAAt                                                        <  2:1656790/93‑1 (MQ=255)

                                   gatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTa                                                      >  2:3263800/1‑93 (MQ=255)

                                   gatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTa                                                      >  1:1975087/1‑93 (MQ=255)

                                    atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  2:41062/93‑1 (MQ=255)

                                    atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  2:133470/93‑1 (MQ=255)

                                    atgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTaa                                                     <  1:1451165/93‑1 (MQ=255)

                                     tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  2:2568763/93‑1 (MQ=255)

                                     tgTCTTTAAATTCAAAAGAATGTCCGTTAAAAAATTAACTATCCCCTTACCGGTGGCCTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                    <  2:4626282/93‑1 (MQ=255)

                                       tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAggg                                                  <  2:599324/93‑1 (MQ=255)

                                       tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAggg                                                  <  2:3674486/93‑1 (MQ=255)

                                        cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGt                                                 >  2:4177275/1‑93 (MQ=255)

                                         tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  1:2105568/1‑93 (MQ=255)

                                         tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  1:4706862/1‑93 (MQ=255)

                                         tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                >  1:4583516/1‑93 (MQ=255)

                                         tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                <  2:1503341/93‑1 (MQ=255)

                                         tttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTg                                                <  2:1503342/93‑1 (MQ=255)

                                          ttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTgg                                               >  2:4245979/1‑93 (MQ=255)

                                           tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              >  1:1262778/1‑93 (MQ=255)

                                           tAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGc                                              <  1:2861869/93‑1 (MQ=255)

                                            aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcg                                             >  1:1121892/1‑92 (MQ=255)

                                            aaaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcg                                             >  1:1943296/1‑92 (MQ=255)

                                             aaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcgc                                            >  1:2193193/1‑91 (MQ=255)

                                             aaTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGcgc                                            >  2:1809078/1‑91 (MQ=255)

                                              aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                           <  2:564418/93‑1 (MQ=255)

                                                tCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAt                                         <  2:168477/93‑1 (MQ=255)

                                                 ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                        <  2:118722/93‑1 (MQ=255)

                                                  cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                       >  1:2098333/1‑93 (MQ=255)

                                                      gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  1:1366243/1‑93 (MQ=255)

                                                      gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  2:3726567/1‑93 (MQ=255)

                                                      gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                   >  2:3726569/1‑93 (MQ=255)

                                                       aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                  >  1:2600323/1‑93 (MQ=255)

                                                       aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCAATTACTGAc                                  >  1:1594512/1‑93 (MQ=255)

                                                         tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATCGCCCGTTGCATTAAGGGTGGCGCAATTCCTGACAt                                >  1:2947616/1‑93 (MQ=255)

                                                         tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                <  2:104890/93‑1 (MQ=255)

                                                         tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                >  2:2155465/1‑93 (MQ=255)

                                                         tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                >  1:2947615/1‑93 (MQ=255)

                                                          gTCCGTTGTAATATTAACTATCCCCATACGGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               <  2:2501728/93‑1 (MQ=255)

                                                          gTCCGTTGAAAGATTAACTACCCCCTGACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                               <  1:1985363/93‑1 (MQ=255)

                                                           tCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAttt                              <  2:4393989/93‑1 (MQ=255)

                                                               ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                          >  2:2584698/1‑93 (MQ=255)

                                                                tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                         <  2:4657032/93‑1 (MQ=255)

                                                                tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                         <  2:4200546/93‑1 (MQ=255)

                                                                  aaaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAc                       <  2:3483046/93‑1 (MQ=255)

                                                                    atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                     >  1:2432424/1‑93 (MQ=255)

                                                                    atatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCg                     >  2:285233/1‑93 (MQ=255)

                                                                      atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCt                   <  1:2605307/93‑1 (MQ=255)

                                                                      atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGCGGCGCCATTACTGACATTTACGATACCGCt                   <  1:2605306/93‑1 (MQ=255)

                                                                       ttAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTa                  <  1:3835231/93‑1 (MQ=255)

                                                                       ttAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTa                  >  2:1150809/1‑93 (MQ=255)

                                                                         aaCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAcc                >  1:3481357/1‑93 (MQ=255)

                                                                          aCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa               >  2:3559486/1‑93 (MQ=255)

                                                                          aCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa               >  2:768685/1‑93 (MQ=255)

                                                                           cTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaa              >  2:2914461/1‑93 (MQ=255)

                                                                              tCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAt           <  1:3147934/93‑1 (MQ=255)

                                                                                cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa         <  1:3758896/93‑1 (MQ=255)

                                                                                cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa         <  1:4418314/93‑1 (MQ=255)

                                                                                cccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTa         <  2:1588466/93‑1 (MQ=255)

                                                                                 ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        <  1:4094287/93‑1 (MQ=255)

                                                                                 ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc        <  1:4250926/93‑1 (MQ=255)

                                                                                  cTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAcc       <  1:85332/93‑1 (MQ=255)

                                                                                   ttACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg      >  2:2789267/1‑93 (MQ=255)

                                                                                   ttACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCg      >  2:534990/1‑93 (MQ=255)

                                                                                    tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     <  2:2023309/93‑1 (MQ=255)

                                                                                    tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     <  2:693126/93‑1 (MQ=255)

                                                                                    tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     <  2:1687580/93‑1 (MQ=255)

                                                                                    tACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGAGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGa     >  1:3715612/1‑93 (MQ=255)

                                                                                      ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  2:2101743/93‑1 (MQ=255)

                                                                                      ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   <  2:912109/93‑1 (MQ=255)

                                                                                      ccGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa   >  2:4044077/1‑93 (MQ=255)

ATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAT  >  NC_000913/1594430‑1594609

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`