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	96.8%	473.8 / 3.4	157	V598V (GTA→GTG)	b1645	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/2); new base G (66/86); total (69/88)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 162 total aligned reads displayed.

TTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACAG  >  NC_000913/1721849‑1722028

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

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

ttaCGCGGGGGTGTGGTGCTGCTGAACTGTTTTCATGTTGTTTAGCTATAGGGACACTGGGAAGAGCGGTCCGATCAGTTTTCGCGTGTGCgg                                                                                         >  2:3431881/1‑93 (MQ=255)

 taCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCggg                                                                                        <  1:1716272/93‑1 (MQ=255)

 taCGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCggg                                                                                        <  2:3637505/93‑1 (MQ=255)

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

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

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

   cGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGAt                                                                                      >  2:852586/1‑93 (MQ=255)

    gCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGaca                                                                                     <  1:1332042/93‑3 (MQ=255)

    gCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATa                                                                                     <  2:2121892/93‑1 (MQ=255)

      tGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAAc                                                                                   <  1:3316485/93‑1 (MQ=255)

       ggggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACt                                                                                  <  2:1814297/93‑1 (MQ=255)

         ggTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGt                                                                                >  2:78676/1‑93 (MQ=255)

          gtgtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTa                                                                               <  2:322175/93‑1 (MQ=255)

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

           tgtAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTAt                                                                              >  1:2364241/1‑93 (MQ=255)

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

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

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

                 gctgctGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACt                                                                        <  1:533156/93‑1 (MQ=255)

                  ctgctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTg                                                                       >  1:4003558/1‑93 (MQ=255)

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

                    gctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGtt                                                                     <  1:4262098/93‑1 (MQ=255)

                    gctgAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGtt                                                                     <  2:4013395/93‑1 (MQ=255)

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

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

                       gAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCg                                                                  >  2:114843/1‑93 (MQ=255)

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

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

                           tGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGa                                                              <  1:977802/93‑1 (MQ=255)

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

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

                           tGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGa                                                              <  1:4346578/93‑1 (MQ=255)

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

                           tGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGCATAACTGTATTTCACTGTTGCGGGGa                                                              <  2:3051256/93‑1 (MQ=255)

                            gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAg                                                             <  1:4434149/93‑1 (MQ=255)

                            gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAg                                                             <  2:2047310/93‑1 (MQ=255)

                            gTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAg                                                             >  2:4529888/1‑93 (MQ=255)

                             ttCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGt                                                            <  1:2157572/93‑1 (MQ=255)

                               ctcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtga                                                          >  2:3678470/1‑93 (MQ=255)

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

                                tcATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGtgat                                                         >  1:1777590/1‑93 (MQ=255)

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

                                   tgttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgag                                                      <  1:2404615/93‑1 (MQ=255)

                                   tgttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgag                                                      >  2:1601657/1‑93 (MQ=255)

                                   tgttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgag                                                      <  2:261375/93‑1 (MQ=255)

                                   tgttgtCTGGCAATTGCGCGACTGGGAATCGCGGTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgag                                                      <  1:4566839/93‑1 (MQ=255)

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

                                     ttgtCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGAtgagtg                                                    >  2:4446596/1‑93 (MQ=255)

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

                                        tCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGc                                                 <  2:1292344/93‑1 (MQ=255)

                                        tCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGc                                                 <  2:4176839/93‑1 (MQ=255)

                                         cTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCAAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCg                                                <  1:1649754/93‑1 (MQ=255)

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

                                             gaaTTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGc                                            <  2:1242153/92‑1 (MQ=255)

                                             cAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGc                                            >  1:2374479/1‑93 (MQ=255)

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

                                                ttGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGtt                                         >  1:3284668/1‑93 (MQ=255)

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

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

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

                                                    gcgACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCCGc                                     >  2:3249103/1‑93 (MQ=255)

                                                     cgACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAGCTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGGGTGGCGTTCAGCa                                    >  1:1089962/1‑93 (MQ=255)

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

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

                                                          gggAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAAt                               >  2:539809/1‑93 (MQ=255)

                                                          gggAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAAt                               <  2:4668391/93‑1 (MQ=255)

                                                          gggAATCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAAt                               >  2:2180639/1‑93 (MQ=255)

                                                             aaTCGCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCAc                            <  1:1654204/93‑1 (MQ=255)

                                                              aTCTCGTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                           <  1:100294/93‑1 (MQ=255)

                                                              aTCGGGTTCGGATCCTTTATCGCGAGTGCGGGATAACTGTATTTCACTGGTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACt                           <  1:3673183/93‑1 (MQ=255)

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

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

                                                                 gcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggg                        >  1:2634147/1‑92 (MQ=255)

                                                                 gcgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggc                        <  1:886839/93‑1 (MQ=255)

                                                                  cgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcg                       <  1:3369180/93‑1 (MQ=255)

                                                                  cgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcg                       <  1:3369179/93‑1 (MQ=255)

                                                                  cgTTCCGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTggcg                       <  2:334711/93‑1 (MQ=255)

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

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

                                                                    ttCCGATCCGTTATCGCGAGAGCGGGATAACTGTATTTCACTGCTGCGGGGAGCGATGAGTGAGCGTGGCGGTCAGCAAAAATGACTggcggc                     >  2:4234291/1‑93 (MQ=255)

                                                                     tcggATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGcc                    <  1:4342577/90‑1 (MQ=255)

                                                                       cGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCcac                  <  1:4286580/93‑1 (MQ=255)

                                                                       cGATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCcac                  >  1:1184515/1‑93 (MQ=255)

                                                                        gATCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCcaca                 >  1:2871896/1‑93 (MQ=255)

                                                                          tCCGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACt               <  1:1177142/93‑1 (MQ=255)

                                                                            cGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTg             <  2:1460676/93‑1 (MQ=255)

                                                                            cGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTg             >  1:4729297/1‑93 (MQ=255)

                                                                            cGTTATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTg             >  2:2889711/1‑93 (MQ=255)

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

                                                                               tATCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaag          >  2:3378090/1‑93 (MQ=255)

                                                                                aTCGCGAGTGCGGGATAACTGAATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaaga         >  2:4672624/1‑93 (MQ=255)

                                                                                 tCGCGAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTgaagaa        <  1:32394/93‑1 (MQ=255)

                                                                                   gcgAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAAtt      <  2:1116478/93‑1 (MQ=255)

                                                                                     gAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTaa    >  1:2111184/1‑92 (MQ=255)

                                                                                     gAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTaa    >  2:3445645/1‑92 (MQ=255)

                                                                                     gAGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTAc    <  1:4151661/93‑1 (MQ=255)

                                                                                      aGTGCGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAATAATCACTGGCGGCCACACTTGAAGAATTACa   >  2:3213874/1‑93 (MQ=255)

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

TTACGCTGGGGTGTAGTGCTGCTGAACTGTTCTCATGTTGTCTGGCAATTGCGCGACTGGGAATCGCGTTCCGATCCGTTATCGCGAGTACGGGATAACTGTATTTCACTGTTGCGGGGAGTGATGAGTGAGCGTGGCGTTCAGCAAAAATCACTGGCGGCCACACTTGAAGAATTACAG  >  NC_000913/1721849‑1722028

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