BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,698,875	A→G	V499V (GTA→GTG)	malX →	MalX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,698,875	0	A	G	98.6%	472.2 / ‑4.7	147	V499V (GTA→GTG)	malX	MalX
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (71/74); total (72/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 153 total aligned reads displayed.

TAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTACAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTGTCC  >  NC_000913/1698785‑1698963

  aCTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAAc                                                                                      <  2:3952308/93‑1 (MQ=255)

  aCTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAAc                                                                                      >  1:232948/1‑93 (MQ=255)

   cTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACt                                                                                     >  2:1670345/1‑93 (MQ=255)

   cTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACt                                                                                     >  2:4727183/1‑93 (MQ=255)

    tGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGTAGGACAATAAGGCACTTGGCGTAGTGCAACtt                                                                                    >  2:3725916/1‑93 (MQ=255)

    tGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACtt                                                                                    >  1:622848/1‑93 (MQ=255)

    tGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTACAAcgt                                                                                    <  2:4103300/93‑3 (MQ=255)

      cATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTaa                                                                                  >  2:607767/1‑93 (MQ=255)

       aTTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAAt                                                                                 <  2:735061/93‑1 (MQ=255)

        ttACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATc                                                                                <  1:1940396/93‑1 (MQ=255)

        ttACCCGTCTGCGTTTGACTGTGAAAGATATGTCGCTTGTTGATGTGCAGGCACTGAAGGCCACTCGGGCAATAGGCGTAGTGCAACTTAATc                                                                                <  1:4198437/93‑1 (MQ=255)

           cccGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAAc                                                                             >  2:1980710/1‑93 (MQ=255)

           cccGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAAc                                                                             <  2:3233209/93‑1 (MQ=255)

              gTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATa                                                                          >  1:4631412/1‑93 (MQ=255)

               tCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATaa                                                                         <  2:4208255/93‑1 (MQ=255)

               tCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGAACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATaa                                                                         <  1:1194806/93‑1 (MQ=255)

                cTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAAc                                                                        <  1:4451093/93‑1 (MQ=255)

                cTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAAc                                                                        <  1:1571659/93‑1 (MQ=255)

                 tGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAAcc                                                                       <  2:4269585/93‑1 (MQ=255)

                 tGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAAcc                                                                       <  2:243609/93‑1 (MQ=255)

                  gCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCt                                                                      >  1:3168958/1‑93 (MQ=255)

                   cGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTg                                                                     <  1:2380023/93‑1 (MQ=255)

                    gTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGc                                                                    >  2:17021/1‑93 (MQ=255)

                    gTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGc                                                                    >  1:1353019/1‑93 (MQ=255)

                    gTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGc                                                                    >  2:4349933/1‑93 (MQ=255)

                     tttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   <  1:583803/93‑1 (MQ=255)

                     tttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   >  2:2886729/1‑93 (MQ=255)

                     tttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   <  2:1331931/93‑1 (MQ=255)

                     tttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   <  2:1331932/93‑1 (MQ=255)

                     tttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   >  1:1535799/1‑93 (MQ=255)

                     gttGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCa                                                                   <  1:2653997/92‑1 (MQ=255)

                       tGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAgg                                                                 <  1:1164714/93‑1 (MQ=255)

                       tGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAgg                                                                 <  2:1985861/93‑1 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                >  1:3364051/1‑93 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                >  1:2608286/1‑93 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                <  1:3726535/93‑1 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                >  1:1130781/1‑93 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                >  1:1050918/1‑93 (MQ=255)

                        gTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgt                                                                <  1:608928/93‑1 (MQ=255)

                         tCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgtt                                                               >  1:151172/1‑93 (MQ=255)

                           tgtgAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgttgt                                                             <  1:1933335/93‑1 (MQ=255)

                           tgtgAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgttgt                                                             <  1:801551/93‑1 (MQ=255)

                            gtgAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGgttgtt                                                            <  1:3410367/93‑1 (MQ=255)

                             tgAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTa                                                           >  2:869638/1‑93 (MQ=255)

                              gAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTAt                                                          >  1:2486795/1‑93 (MQ=255)

                               aaaGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATc                                                         <  1:1864955/93‑1 (MQ=255)

                               aaaGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATc                                                         >  2:1177119/1‑93 (MQ=255)

                               aaaGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCAGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATc                                                         <  1:2948374/93‑1 (MQ=255)

                                aaGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCg                                                        >  1:1857269/1‑93 (MQ=255)

                                aaGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGGTGTTATCg                                                        >  2:3128827/1‑93 (MQ=255)

                                   atatGTCGCTTGTTAATGTGCAGGCNCNGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGc                                                     <  2:2323800/93‑1 (MQ=255)

                                    tatGTCGGTTGTTAATGTGCAGGCACTGAAGGGCAATCGGGCAATTAGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGTcc                                                    >  2:4147155/1‑93 (MQ=255)

                                    tatGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGcc                                                    <  2:4651212/93‑1 (MQ=255)

                                     atGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCca                                                   >  2:1782943/1‑93 (MQ=255)

                                      tGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCcac                                                  <  2:4158644/93‑1 (MQ=255)

                                        tCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACaa                                                <  2:3022366/93‑1 (MQ=255)

                                        tCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACaa                                                >  1:2153729/1‑93 (MQ=255)

                                          gCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGt                                              <  2:2465626/93‑1 (MQ=255)

                                            ttGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtc                                            >  2:3812131/1‑93 (MQ=255)

                                             tGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtca                                           >  2:2678554/1‑93 (MQ=255)

                                             tGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtca                                           <  1:1669993/93‑1 (MQ=255)

                                             tGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtca                                           <  1:81958/93‑1 (MQ=255)

                                              gTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtcag                                          <  2:2383520/93‑1 (MQ=255)

                                                 aaTGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtcagtca                                       <  2:1902676/93‑1 (MQ=255)

                                                  aTGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtcagtcag                                      >  2:4702378/1‑93 (MQ=255)

                                                  aTGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTtcagtcag                                      >  2:4502793/1‑93 (MQ=255)

                                                     tgCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTaa                                   <  1:526164/93‑1 (MQ=255)

                                                     tgCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGAAACTTAATCAACATAACCTGCAGGTTATTATCGGGCCACAAGTTCAGTCAGTaa                                   <  1:2419829/93‑1 (MQ=255)

                                                     tgCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTACAACTTAATCAACATAACCTGCAGGTTGTTATCGTGCCACAAGTTCAGTCAGTaa                                   >  1:3819258/1‑93 (MQ=255)

                                                       cAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTaaaa                                 >  1:2978611/1‑93 (MQ=255)

                                                         ggCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTAAGTCAGTAAAAGa                               >  2:1435453/1‑93 (MQ=255)

                                                          gCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGAt                              <  2:2516079/93‑1 (MQ=255)

                                                          gCACTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGAt                              <  2:2892005/93‑1 (MQ=255)

                                                             cTGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGaa                           <  2:2587404/93‑1 (MQ=255)

                                                              tGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGaag                          <  1:2901294/93‑2 (MQ=255)

                                                              tGAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGaaa                          <  2:865561/93‑1 (MQ=255)

                                                               gAAGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAAt                         >  1:2373227/1‑93 (MQ=255)

                                                                aaGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATg                        <  2:1249628/93‑1 (MQ=255)

                                                                aaGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATg                        <  1:3746666/93‑1 (MQ=255)

                                                                aaGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATg                        >  1:4428980/1‑93 (MQ=255)

                                                                 aGGACAATCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATgg                       >  2:1520434/1‑93 (MQ=255)

                                                                      aaTCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGt                  >  1:2753183/1‑93 (MQ=255)

                                                                      aaTCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGt                  >  2:2178782/1‑93 (MQ=255)

                                                                       aTCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTc                 <  1:1017530/93‑1 (MQ=255)

                                                                       aTCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTc                 >  2:4089940/1‑93 (MQ=255)

                                                                        tCGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCt                >  2:222159/1‑93 (MQ=255)

                                                                         cGGGCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTg               <  2:2339234/93‑1 (MQ=255)

                                                                           ggCAATTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGAt             >  2:308754/1‑93 (MQ=255)

                                                                              aaTTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCa          >  1:3141508/1‑93 (MQ=255)

                                                                               aTTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCAt         >  1:3551716/1‑93 (MQ=255)

                                                                               aTTGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGACGCAt         >  1:3570085/1‑93 (MQ=255)

                                                                                ttGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATa        <  2:2761288/93‑1 (MQ=255)

                                                                                ttGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATa        <  1:90265/93‑1 (MQ=255)

                                                                                 tGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATAc       <  2:3634682/93‑1 (MQ=255)

                                                                                 tGGCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATAc       <  1:2002292/93‑1 (MQ=255)

                                                                                   gCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTg     >  2:3934739/1‑93 (MQ=255)

                                                                                   gCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTg     <  1:3349166/93‑1 (MQ=255)

                                                                                   gCGTAGTGCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTg     <  1:1354853/93‑1 (MQ=255)

                                                                                     gtagtgCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTGTc   >  1:1255200/1‑93 (MQ=255)

                                                                                     gtagtgCAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTGTc   >  1:4577838/1‑93 (MQ=255)

TAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTAATGTGCAGGCACTGAAGGACAATCGGGCAATTGGCGTAGTACAACTTAATCAACATAACCTGCAGGTTGTTATCGGGCCACAAGTTCAGTCAGTAAAAGATGAAATGGCCGGTCTGATGCATACTGTCC  >  NC_000913/1698785‑1698963

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