BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	423,381	T→C	V548A (GTG→GCG)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	423,381	0	T	C	98.8%	488.5 / ‑5.0	161	V548A (GTG→GCG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (73/85); total (73/88)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 175 total aligned reads displayed.

TGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGTG  >  NC_000913/423291‑423469

tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGGAGATAACGTGGTGGTATTTGTCCgcgcgc                                                                                        >  1:3382698/1‑93 (MQ=255)

 gATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCt                                                                                       <  1:104727/93‑1 (MQ=255)

 gATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCt                                                                                       >  1:1000502/1‑93 (MQ=255)

   ttGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGa                                                                                     <  2:4419227/93‑1 (MQ=255)

   ttGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGa                                                                                     <  2:615334/93‑1 (MQ=255)

    tGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGaa                                                                                    <  2:721010/93‑1 (MQ=255)

     tcgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAAt                                                                                   <  2:4331214/92‑1 (MQ=255)

      cgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATc                                                                                  <  2:786042/93‑1 (MQ=255)

      cgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATc                                                                                  <  1:4152783/93‑1 (MQ=255)

      cgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATc                                                                                  <  1:2430678/93‑1 (MQ=255)

       gcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCa                                                                                 <  2:758801/93‑1 (MQ=255)

       gcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCa                                                                                 >  2:1165981/1‑93 (MQ=255)

       gcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCa                                                                                 <  2:2199773/93‑1 (MQ=255)

       gcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCa                                                                                 >  2:2326929/1‑93 (MQ=255)

        cTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAg                                                                                <  1:1098818/93‑1 (MQ=255)

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTTCGCGCGCTGAATCAGc                                                                               >  2:3591685/1‑93 (MQ=255)

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGc                                                                               <  1:1820112/93‑1 (MQ=255)

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGc                                                                               <  1:2299426/93‑1 (MQ=255)

           cGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCaa                                                                             <  1:3980159/93‑1 (MQ=255)

           cGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCaa                                                                             <  2:131886/93‑1 (MQ=255)

           cGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCaa                                                                             <  2:2460753/93‑1 (MQ=255)

           cGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCaa                                                                             <  2:3523179/93‑1 (MQ=255)

           cGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCaa                                                                             >  2:4105004/1‑93 (MQ=255)

            gTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGGAAc                                                                            >  2:2785380/1‑93 (MQ=255)

            gTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAAc                                                                            <  2:2548905/93‑1 (MQ=255)

            gTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAAc                                                                            >  2:947316/1‑93 (MQ=255)

            gTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAAc                                                                            <  1:3514586/93‑1 (MQ=255)

             tAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACg                                                                           <  2:1722815/93‑1 (MQ=255)

             tAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGGATTTGTCCGCGCGCTGAATCAGCAACg                                                                           >  2:1450474/1‑93 (MQ=255)

              aaGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACgt                                                                          <  2:600139/93‑1 (MQ=255)

              aaGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGGATTTGTCCGCGCGCTGAATCAGCAACgt                                                                          >  2:4645066/1‑93 (MQ=255)

               aGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACgtg                                                                         >  1:3561873/1‑93 (MQ=255)

                gAAAAGACAGGCGCTACGTCATGGCGGCTGTCAGGGGCTGTATGCGGAAGATAACGTGGTGGTACTTGTCCGCGCGCCGAATCAGTAACgtgt                                                                        >  2:2410441/1‑93 (MQ=255)

                 aaaaGTCAGGCTCTACGTCATGGCGGCGGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                       <  1:4562250/93‑1 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACgtgtc                                                                       >  1:458106/1‑92 (MQ=255)

                  aaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTAc                                                                      >  1:2781121/1‑93 (MQ=255)

                  aaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTAc                                                                      <  2:2131940/93‑1 (MQ=255)

                      tCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACtggt                                                                  >  2:3949513/1‑93 (MQ=255)

                         ggCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGc                                                               <  1:3985383/93‑1 (MQ=255)

                         ggCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGc                                                               >  1:314101/1‑93 (MQ=255)

                         ggCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGc                                                               >  2:1213069/1‑93 (MQ=255)

                            gcTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAAt                                                            <  1:3712277/93‑1 (MQ=255)

                            gcTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAAt                                                            >  2:3761628/1‑93 (MQ=255)

                             cTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATc                                                           >  1:52392/1‑93 (MQ=255)

                             cTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATc                                                           >  1:2174150/1‑93 (MQ=255)

                             cTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATc                                                           <  1:3470019/93‑1 (MQ=255)

                             cTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCCCGCTGAATCAGCAACGTGTACTGGTGGCAATc                                                           >  2:2048586/1‑93 (MQ=255)

                               aCGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCaa                                                         <  1:1419584/93‑1 (MQ=255)

                               aCGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCaa                                                         >  1:941733/1‑93 (MQ=255)

                               aCGTCATGGCGCCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGGGCTGAATCAGCAACGTGTACTGGTGGCAATCat                                                         <  2:3201685/93‑2 (MQ=255)

                                cGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAAc                                                        <  1:3846839/93‑1 (MQ=255)

                                cGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAAc                                                        >  1:4255445/1‑93 (MQ=255)

                                  tCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCg                                                      <  1:1342332/93‑1 (MQ=255)

                                   cATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGt                                                     <  1:1566389/93‑1 (MQ=255)

                                    aTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTg                                                    <  1:4227973/93‑1 (MQ=255)

                                      ggcggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGc                                                  <  2:3388205/93‑1 (MQ=255)

                                        gggcTGTGAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGa                                                <  1:2502274/92‑1 (MQ=255)

                                        cggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGa                                                <  2:1018934/93‑1 (MQ=255)

                                        cggcTGTATGGTGCGGTATTCGAAAGATTACGTGGTGGTATTTGTCCGCTCTATGAATGAGCAACGGTTACTGTTGGCAATCAGCCGTGGcta                                                <  2:3664274/93‑3 (MQ=255)

                                         ggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAg                                               <  2:4347160/93‑1 (MQ=255)

                                         ggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAg                                               <  1:2639554/93‑1 (MQ=255)

                                          gcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCCATCAACCGTGGCGAgg                                              >  2:553446/1‑93 (MQ=255)

                                          gcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAgg                                              <  1:3930776/93‑1 (MQ=255)

                                           cTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGc                                             <  1:752123/93‑1 (MQ=255)

                                            tGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGcc                                            <  2:358773/93‑1 (MQ=255)

                                              tCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtg                                          >  1:2966732/1‑93 (MQ=255)

                                               cAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtgt                                         <  1:3676716/93‑1 (MQ=255)

                                                 ggTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGa                                       >  2:2851841/1‑93 (MQ=255)

                                                 ggTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGa                                       >  2:1666433/1‑93 (MQ=255)

                                                 ggTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGa                                       >  1:1877558/1‑93 (MQ=255)

                                                  gTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGaa                                      <  2:843694/93‑1 (MQ=255)

                                                  gTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGaa                                      >  1:4308885/1‑93 (MQ=255)

                                                    gCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgt                                    >  1:680312/1‑93 (MQ=255)

                                                    gCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgt                                    <  2:774792/93‑1 (MQ=255)

                                                      tGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtgg                                  <  1:3222050/93‑1 (MQ=255)

                                                      tGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtgg                                  <  1:2202971/93‑1 (MQ=255)

                                                         aTGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGc                               <  2:1759262/93‑1 (MQ=255)

                                                         aTGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGc                               <  1:1788856/93‑1 (MQ=255)

                                                            cGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTAc                            <  1:471313/93‑1 (MQ=255)

                                                               aaGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCcg                         >  2:3812571/1‑93 (MQ=255)

                                                                  aTAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGt                      <  2:1066319/93‑1 (MQ=255)

                                                                  aTAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGt                      >  2:1486710/1‑93 (MQ=255)

                                                                  aTAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCAGCGt                      <  1:3680267/93‑1 (MQ=255)

                                                                   tAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTc                     >  2:306006/1‑93 (MQ=255)

                                                                     aCGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCAc                   >  1:715101/1‑93 (MQ=255)

                                                                     aCGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCAc                   <  2:1572540/93‑1 (MQ=255)

                                                                     aCGTGGTGGTATTTGTCCGCGCGCTCAATCAGCAACGTGTACTGGTGGCAATCAACCGCGGCGAGGCCTGGGCAGGGGTGCTACCCGCGTCAc                   >  2:3221384/1‑93 (MQ=255)

                                                                       gtggtgGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCg                 >  2:1952636/1‑93 (MQ=255)

                                                                         ggtggtATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGtt               >  2:2055284/1‑93 (MQ=255)

                                                                           tggtATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGCGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTtc             >  1:780939/1‑93 (MQ=255)

                                                                            ggtATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTtct            >  1:2106614/1‑93 (MQ=255)

                                                                              tATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCa          >  1:3833171/1‑93 (MQ=255)

                                                                               atgtGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGGGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCaa         <  2:3955844/90‑1 (MQ=255)

                                                                               aTTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCaa         >  2:4633914/1‑93 (MQ=255)

                                                                               aTTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGGGCTACCCGCGTCACCGTTTCTCaa         <  2:4517444/93‑1 (MQ=255)

                                                                               aTTTGTCCGCGCGCTGAAACAGCCGCGTGAAATGGTGACCCTCAACCGTGGAAAGGCCTATGAAGTGATGCTGCACACGTCACCGTTTCTCaa         <  2:4685550/93‑1 (MQ=255)

                                                                                tttGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAAt        <  1:1361300/93‑1 (MQ=255)

                                                                                   gTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGcc     <  1:2201724/93‑1 (MQ=255)

                                                                                     ccGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGt   <  2:3434157/93‑1 (MQ=255)

                                                                                      cgcgcgCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGTCGTg  >  2:1278968/1‑93 (MQ=255)

TGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGTG  >  NC_000913/423291‑423469

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