BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,721,938	A→G	V598V (GTA→GTG)	b1645 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,721,938	0	A	G	99.5%	605.4 / ‑5.2	196	V598V (GTA→GTG)	b1645	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (103/92); total (103/93)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 205 total aligned reads displayed.

ATTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACA  >  NC_000913/1721848‑1722027

attaCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCg                                                                                         >  1:3551346/1‑93 (MQ=255)

 ttaCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCgg                                                                                        <  2:6247187/93‑1 (MQ=255)

 ttaCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCgg                                                                                        >  2:5161369/1‑93 (MQ=255)

   aCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGa                                                                                      <  2:6663238/93‑1 (MQ=255)

    cGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGAt                                                                                     <  1:3412583/93‑1 (MQ=255)

    cGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGAt                                                                                     <  1:858696/93‑1 (MQ=255)

     gCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATa                                                                                    >  1:6570596/1‑93 (MQ=255)

      cTGGGGTTTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATaa                                                                                   <  2:1516472/93‑1 (MQ=255)

       tGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAAc                                                                                  >  1:1984868/1‑93 (MQ=255)

       tGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAAc                                                                                  >  2:38635/1‑93 (MQ=255)

         gggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTg                                                                                >  2:6446771/1‑93 (MQ=255)

         gggTGTAGAGCTGCTGAACTGGTCTCATGTTGTCTGGCAATTTCGCGACTGGGGATCGCGAGCCGATACGTTATCGCGAGTGCGGGATAACTg                                                                                >  2:5170411/1‑93 (MQ=255)

          ggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGt                                                                               <  2:6268911/93‑1 (MQ=255)

            tgtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAt                                                                             <  2:5326773/93‑1 (MQ=255)

             gtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAtt                                                                            <  2:4232624/93‑1 (MQ=255)

             gtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAtt                                                                            <  1:3253194/93‑1 (MQ=255)

              tAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAttt                                                                           >  1:6784296/1‑93 (MQ=255)

              tAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCCTTATCGCGAGTGCGGGATAACTGTAttt                                                                           >  1:7443884/1‑93 (MQ=255)

                gTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCa                                                                         >  2:1243032/1‑93 (MQ=255)

                gTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCa                                                                         <  2:1027416/93‑1 (MQ=255)

                 tgctgcTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCAc                                                                        <  1:7228922/93‑1 (MQ=255)

                 tgctgcTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCAc                                                                        <  1:592626/93‑1 (MQ=255)

                   ctgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTg                                                                      >  1:4831522/1‑93 (MQ=255)

                    tgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGt                                                                     >  2:822162/1‑93 (MQ=255)

                    tgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGt                                                                     >  1:3235480/1‑93 (MQ=255)

                     gctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGtt                                                                    >  1:1069194/1‑93 (MQ=255)

                       tgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGc                                                                  <  2:2771253/93‑1 (MQ=255)

                       tgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGc                                                                  >  1:4530812/1‑93 (MQ=255)

                        gAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCg                                                                 <  2:1446566/93‑1 (MQ=255)

                         aaCTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgg                                                                <  1:291312/93‑1 (MQ=255)

                         aaCTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgg                                                                <  1:2933205/93‑1 (MQ=255)

                         aaCTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgg                                                                >  1:1631390/1‑93 (MQ=255)

                         aaCTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCgg                                                                <  1:5481611/93‑1 (MQ=255)

                            tGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGa                                                             >  2:5707391/1‑93 (MQ=255)

                            tGTTCTAATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGGTCCGATCCGTTATCGCGAGTGCGGGCTAACTGTATTTCACTGCTGCGTGGa                                                             >  2:7804369/1‑93 (MQ=255)

                             gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAg                                                            >  2:3355972/1‑93 (MQ=255)

                               tctcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtg                                                          >  1:3048019/1‑93 (MQ=255)

                                ctcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtga                                                         <  1:3988785/93‑1 (MQ=255)

                                 tcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgat                                                        <  2:3578884/93‑1 (MQ=255)

                                  cATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgatg                                                       >  1:502821/1‑93 (MQ=255)

                                   aTGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgatga                                                      >  1:5782572/1‑93 (MQ=255)

                                   aTGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgatga                                                      <  1:5109445/93‑1 (MQ=255)

                                     gttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagt                                                    >  2:6726102/1‑93 (MQ=255)

                                     gttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagt                                                    <  1:4840336/93‑1 (MQ=255)

                                       tgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  <  1:3345523/93‑1 (MQ=255)

                                       tgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  >  2:1630217/1‑93 (MQ=255)

                                        gtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtgag                                                 >  2:2984021/1‑93 (MQ=255)

                                        gtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtgag                                                 >  2:2413143/1‑93 (MQ=255)

                                          cTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCg                                               >  1:221086/1‑93 (MQ=255)

                                           tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCgg                                              >  1:1828603/1‑92 (MQ=255)

                                           tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              >  1:2543686/1‑93 (MQ=255)

                                           tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              >  2:3768007/1‑93 (MQ=255)

                                           tGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              >  2:5010387/1‑93 (MQ=255)

                                           tGGCAATTGCACGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGt                                              <  1:3932114/93‑1 (MQ=255)

                                            ggCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTg                                             <  1:5742570/93‑1 (MQ=255)

                                            ggCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTg                                             >  1:7512616/1‑93 (MQ=255)

                                               aaTTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCg                                          >  1:2803197/1‑93 (MQ=255)

                                               aaTTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCg                                          <  2:724364/93‑1 (MQ=255)

                                                  tGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTc                                       <  1:55987/93‑1 (MQ=255)

                                                  tGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTc                                       >  2:4667119/1‑93 (MQ=255)

                                                   gcgcGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCa                                      <  2:3529420/93‑1 (MQ=255)

                                                     gcgACTGGGAAGCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGc                                    <  1:4108702/93‑1 (MQ=255)

                                                       gACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaa                                  >  2:6817081/1‑93 (MQ=255)

                                                        aCTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaa                                 >  2:4845080/1‑93 (MQ=255)

                                                        aCTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaa                                 >  1:3279159/1‑93 (MQ=255)

                                                        aCTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaa                                 >  1:2040328/1‑93 (MQ=255)

                                                         cTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaaa                                >  1:3622275/1‑93 (MQ=255)

                                                          tGGGAATCGCGTTGCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaaaa                               <  2:7635764/93‑1 (MQ=255)

                                                          tGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaaaa                               <  2:1777435/93‑1 (MQ=255)

                                                          tGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCaaaaa                               >  1:6893687/1‑93 (MQ=255)

                                                             gAATCGCGCTCCGATCCGTTATCGCGAGTGCGCGATAACTGTATTTCGCTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCa                            <  1:608958/93‑1 (MQ=255)

                                                               aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          <  2:2223418/93‑1 (MQ=255)

                                                               aTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          <  1:6338250/93‑1 (MQ=255)

                                                               aTCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                          <  1:2988966/93‑1 (MQ=255)

                                                                tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         <  1:5071810/93‑1 (MQ=255)

                                                                tCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTg                         >  2:7296911/1‑93 (MQ=255)

                                                                 cgcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTgg                        <  2:6621205/93‑1 (MQ=255)

                                                                  gcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTTTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggc                       <  2:6907996/93‑1 (MQ=255)

                                                                  gcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggc                       >  2:2923812/1‑93 (MQ=255)

                                                                  gcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggc                       >  2:7343834/1‑93 (MQ=255)

                                                                    gTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcgg                     >  1:3640264/1‑93 (MQ=255)

                                                                     ttCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcggc                    >  1:3569272/1‑93 (MQ=255)

                                                                     ttCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcggc                    >  1:4380503/1‑93 (MQ=255)

                                                                     ttCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcggc                    >  1:5008209/1‑93 (MQ=255)

                                                                      tCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGcc                   >  2:1220703/1‑93 (MQ=255)

                                                                      tCCGATCCGTTAGCGCGAGTGCGGGATAACTGGATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGcc                   <  2:5950093/93‑1 (MQ=255)

                                                                       ccGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCca                  >  1:3374465/1‑93 (MQ=255)

                                                                          aTCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCacac               >  1:2686226/1‑93 (MQ=255)

                                                                          aTCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCacac               <  2:4132281/93‑1 (MQ=255)

                                                                            ccGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACtt             <  2:1739806/93‑1 (MQ=255)

                                                                               ttATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  >  2:2858158/1‑53 (MQ=255)

                                                                               ttATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtga                                                  <  1:2858159/53‑1 (MQ=255)

                                                                               ttATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaa          <  1:6480455/93‑1 (MQ=255)

                                                                               ttATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaa          <  2:2154823/93‑1 (MQ=255)

                                                                                 aTCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGTCACACTTgaaga        <  2:7134562/93‑1 (MQ=255)

                                                                                 aTCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaaga        <  1:38635/93‑1 (MQ=255)

                                                                                 aTCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaaga        >  2:2632307/1‑93 (MQ=255)

                                                                                    gcgAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAAtt     >  1:2296391/1‑93 (MQ=255)

                                                                                       aGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACa  >  2:552264/1‑93 (MQ=255)

                                                                                       aGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACa  <  1:2585772/93‑1 (MQ=255)

ATTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACA  >  NC_000913/1721848‑1722027

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