BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	384,546	A→G	T50A (ACC→GCC)	tauA →	TauA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	384,546	0	A	G	99.5%	639.3 / ‑4.6	195	T50A (ACC→GCC)	tauA	TauA
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (97/97); minor base T (1/0); total (98/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

Only 100 of 204 total aligned reads displayed.

ATGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAACCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGC  >  NC_000913/384456‑384635

aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                         <  2:207243/93‑1 (MQ=255)

aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                         >  2:802506/1‑93 (MQ=255)

aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                         <  2:3411626/93‑1 (MQ=255)

 tGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCt                                                                                        >  2:6274017/1‑93 (MQ=255)

  ggCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTc                                                                                       <  2:4490735/93‑1 (MQ=255)

  ggCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTc                                                                                       <  2:5014656/93‑1 (MQ=255)

  ggCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTc                                                                                       >  1:2458275/1‑93 (MQ=255)

     aaTTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGc                                                                                    >  2:417037/1‑93 (MQ=255)

     aaTTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGc                                                                                    >  1:3713552/1‑93 (MQ=255)

      aTTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGcc                                                                                   >  1:2515042/1‑93 (MQ=255)

       tttCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCg                                                                                  >  1:5178239/1‑93 (MQ=255)

       tttCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCg                                                                                  >  1:1394483/1‑93 (MQ=255)

        ttCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGa                                                                                 <  2:389470/93‑1 (MQ=255)

          cATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAc                                                                               <  1:5747712/93‑1 (MQ=255)

           aTCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAcc                                                                              <  1:424807/93‑1 (MQ=255)

           aTCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAcc                                                                              <  1:4026898/93‑1 (MQ=255)

           aTCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAcc                                                                              <  2:183414/93‑1 (MQ=255)

           aTCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAcc                                                                              >  2:5688540/1‑93 (MQ=255)

            tCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCg                                                                             <  2:638349/93‑1 (MQ=255)

             cgcgTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCgg                                                                            >  1:532761/1‑93 (MQ=255)

                gTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGa                                                                         >  1:1867610/1‑93 (MQ=255)

                  aaCACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGaaa                                                                       <  1:4792430/93‑1 (MQ=255)

                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                      <  1:2765733/93‑1 (MQ=255)

                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                      <  1:4940401/93‑1 (MQ=255)

                     acacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTCTCAAGCCTCAGCCGAACCGGCGAAAGTg                                                                    >  2:1151410/1‑93 (MQ=255)

                     acacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTg                                                                    <  1:630845/93‑1 (MQ=255)

                     acacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTg                                                                    <  2:46290/93‑1 (MQ=255)

                      cacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTgg                                                                   <  1:1797454/93‑1 (MQ=255)

                      cacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTgg                                                                   <  1:1035551/93‑1 (MQ=255)

                      cacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTgg                                                                   >  1:3763644/1‑93 (MQ=255)

                      cacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTgg                                                                   <  1:6165586/93‑1 (MQ=255)

                        cttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCt                                                                 >  2:5504078/1‑93 (MQ=255)

                        cttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCt                                                                 <  2:4001340/93‑1 (MQ=255)

                        cttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCt                                                                 >  1:6113119/1‑93 (MQ=255)

                         ttctcGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTc                                                                <  1:5641454/93‑1 (MQ=255)

                          tcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCa                                                               >  2:1425012/1‑93 (MQ=255)

                           cttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAg                                                              >  2:1888235/1‑93 (MQ=255)

                            ttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAgg                                                             >  1:2864695/1‑93 (MQ=255)

                               ccGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCg                                                          <  2:735878/93‑1 (MQ=255)

                                cGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGa                                                         <  1:1792028/93‑1 (MQ=255)

                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGGCGac                                                        >  1:3967614/1‑93 (MQ=255)

                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGac                                                        <  1:4149123/93‑1 (MQ=255)

                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGac                                                        <  1:3534689/93‑1 (MQ=255)

                                  cACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGaca                                                       <  2:2521291/93‑1 (MQ=255)

                                   aCTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaa                                                      >  2:6439681/1‑93 (MQ=255)

                                   aCTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaa                                                      >  1:3321573/1‑93 (MQ=255)

                                   aCTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaa                                                      >  2:6550079/1‑93 (MQ=255)

                                    cTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaac                                                     <  2:6393816/93‑1 (MQ=255)

                                     tGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGCGGCTCAGGCCGacaaca                                                    >  2:361294/1‑93 (MQ=255)

                                      ggCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAacac                                                   <  2:5468048/93‑1 (MQ=255)

                                      ggCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAacac                                                   >  2:2826384/1‑93 (MQ=255)

                                        cATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCt                                                 <  1:3164927/93‑1 (MQ=255)

                                         aTTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCtt                                                >  1:3395854/1‑93 (MQ=255)

                                         aTTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCtt                                                <  1:1480733/93‑1 (MQ=255)

                                            cATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTTGCTCAGGCCGACAACACCTTTGc                                             <  1:5374971/93‑1 (MQ=255)

                                            cATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGc                                             <  2:2948921/93‑1 (MQ=255)

                                            cATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGc                                             >  1:6683468/1‑93 (MQ=255)

                                            cATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGc                                             >  1:1116586/1‑93 (MQ=255)

                                                ccTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaa                                         <  1:1743723/92‑1 (MQ=255)

                                                   tttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaa                                      <  1:4503623/93‑1 (MQ=255)

                                                   tttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaa                                      >  1:4151761/1‑93 (MQ=255)

                                                   tttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaa                                      <  1:2123324/93‑1 (MQ=255)

                                                   tttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaa                                      <  1:1999740/93‑1 (MQ=255)

                                                    ttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaaa                                     <  1:6357567/93‑1 (MQ=255)

                                                    ttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaaa                                     <  2:5755879/93‑1 (MQ=255)

                                                    ttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaaa                                     <  2:2352820/93‑1 (MQ=255)

                                                     tCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaaag                                    <  1:505551/93‑1 (MQ=255)

                                                      cAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGc                                   <  2:5269506/93‑1 (MQ=255)

                                                       aGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCg                                  >  2:4314809/1‑93 (MQ=255)

                                                       aGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCg                                  >  2:3550998/1‑93 (MQ=255)

                                                        ggCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCgg                                 >  2:98267/1‑93 (MQ=255)

                                                        ggCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCgg                                 >  1:131104/1‑93 (MQ=255)

                                                        ggCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCgg                                 <  1:3731705/93‑1 (MQ=255)

                                                        ggCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCgg                                 >  1:5513663/1‑93 (MQ=255)

                                                          cacaGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAg                               <  2:5059946/93‑1 (MQ=255)

                                                             aGGCGGTGAACGTCACCGTGGCGTATCAATCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCaa                            >  2:3757359/1‑93 (MQ=255)

                                                             aGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCac                            >  2:1278907/1‑92 (MQ=255)

                                                             aGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGCCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCaa                            <  1:5961793/93‑1 (MQ=255)

                                                               gCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCCAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCCAcc                          >  2:610118/1‑93 (MQ=255)

                                                                  gTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTg                       <  1:5156953/93‑1 (MQ=255)

                                                                     aaCGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGAc                    >  1:2709165/1‑93 (MQ=255)

                                                                      aCGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACt                   >  1:896434/1‑93 (MQ=255)

                                                                        gTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTgg                 >  1:4987863/1‑93 (MQ=255)

                                                                         tCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGc                >  1:5277321/1‑93 (MQ=255)

                                                                         tCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGc                >  2:2324011/1‑93 (MQ=255)

                                                                         tCACCGTGCCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGc                <  2:3106023/93‑1 (MQ=255)

                                                                           aCCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGt              >  2:170228/1‑93 (MQ=255)

                                                                            ccGTGGCGTATCAAGCCACAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTa             <  1:1419758/93‑1 (MQ=255)

                                                                              gTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAg           <  2:4887127/93‑1 (MQ=255)

                                                                              gTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAg           >  1:3843272/1‑93 (MQ=255)

                                                                               tGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGt          <  1:6552545/93‑1 (MQ=255)

                                                                               tGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGt          >  2:12490/1‑93 (MQ=255)

                                                                                ggCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCGTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGtt         <  2:6043071/93‑1 (MQ=255)

                                                                                   gTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGa      <  2:5881294/93‑1 (MQ=255)

                                                                                   gTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGa      <  2:228534/93‑1 (MQ=255)

                                                                                    tATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGAc     <  2:5678249/93‑1 (MQ=255)

                                                                                    tATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGAc     <  1:5462658/93‑1 (MQ=255)

                                                                                     aTCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACa    <  2:3286990/93‑1 (MQ=255)

                                                                                      tCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAg   <  2:5645455/93‑1 (MQ=255)

                                                                                      tCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAg   >  1:1651628/1‑93 (MQ=255)

ATGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAACCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGC  >  NC_000913/384456‑384635

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