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%	776.9 / ‑2.9	253	V548A (GTG→GCG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (133/117); total (133/120)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 267 total aligned reads displayed.

TGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGTGC  >  NC_000913/423291‑423470

tGATTGCGCTGCGTTAGAAAAGTCCGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGTGTGc                                                                                         >  2:5532814/1‑93 (MQ=255)

tGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCgcgcgc                                                                                         <  1:7712573/93‑1 (MQ=255)

tGATTGCACTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCgcgcgc                                                                                         <  2:3813940/93‑1 (MQ=255)

  aTTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGGATTTGTCCGCGCGCTg                                                                                       >  2:1027297/1‑93 (MQ=255)

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

   ttGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGa                                                                                      >  2:4341012/1‑93 (MQ=255)

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

     gcgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAAt                                                                                    >  2:7591298/1‑93 (MQ=255)

     gcgcTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAAt                                                                                    <  2:6521651/93‑1 (MQ=255)

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

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

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

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

        cTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAg                                                                                 >  1:4327961/1‑93 (MQ=255)

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGc                                                                                >  1:5208241/1‑93 (MQ=255)

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

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

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGc                                                                                >  1:6684063/1‑93 (MQ=255)

         tGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGc                                                                                >  1:4161849/1‑93 (MQ=255)

          gCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCa                                                                               <  1:2075554/93‑1 (MQ=255)

          gCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCa                                                                               >  2:5046316/1‑93 (MQ=255)

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

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

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

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

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

             tAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACg                                                                            >  2:2195574/1‑93 (MQ=255)

             tAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACg                                                                            >  1:5569411/1‑93 (MQ=255)

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

                gAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACgtgt                                                                         <  2:7359140/93‑1 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                        <  2:2853501/93‑1 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                        <  1:7019437/93‑1 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                        <  1:3565557/93‑1 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                        >  2:628772/1‑93 (MQ=255)

                 aaaaGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTa                                                                        <  1:6071031/93‑1 (MQ=255)

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

                     gacAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACtgg                                                                    >  2:7549993/3‑93 (MQ=255)

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

                       cAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACtggtg                                                                  <  2:2684247/93‑1 (MQ=255)

                       cAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACtggtg                                                                  >  2:7711568/1‑93 (MQ=255)

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

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

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

                              tACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCACCGTGTACTGGTGGCAATCa                                                           <  1:6675356/93‑1 (MQ=255)

                              tACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCa                                                           >  1:4821446/1‑93 (MQ=255)

                              tACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCa                                                           >  2:1124335/1‑93 (MQ=255)

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

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

                               aCGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCACGCGCTGAATCAGCAACGTGTACTGGTGGCAATCaa                                                          >  2:1829426/1‑93 (MQ=255)

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

                                   cATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGt                                                      >  1:1709473/1‑93 (MQ=255)

                                    aTGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTg                                                     >  2:160678/1‑93 (MQ=255)

                                     tGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTgg                                                    >  1:7872906/1‑93 (MQ=255)

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

                                      ggcggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGGGGc                                                   >  2:1991491/1‑93 (MQ=255)

                                       gcggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCg                                                  <  1:2129989/93‑1 (MQ=255)

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

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

                                         ggcTGTCAGGTGCTGTATGCGGAAGATAACGTGGGGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAg                                                >  1:6717098/1‑93 (MQ=255)

                                          gcTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAgg                                               >  2:7139258/1‑93 (MQ=255)

                                             gTCAGGTGCTGTATGCGTAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCt                                            >  1:926931/1‑93 (MQ=255)

                                             gTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCt                                            >  1:3040513/1‑93 (MQ=255)

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

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

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

                                               cAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtgt                                          >  2:773850/1‑93 (MQ=255)

                                                aGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtgtg                                         >  2:3510644/1‑93 (MQ=255)

                                                aGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtgtg                                         >  1:4916856/1‑93 (MQ=255)

                                                aGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCtgtg                                         >  1:3371273/1‑93 (MQ=255)

                                                  gcgCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGaa                                       >  2:4755575/3‑93 (MQ=255)

                                                    gCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgg                                     >  1:1920871/1‑92 (MQ=255)

                                                     cTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtg                                    <  1:6626252/93‑1 (MQ=255)

                                                      tGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtgg                                   >  1:2234016/1‑93 (MQ=255)

                                                       gTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtggt                                  >  2:558904/1‑93 (MQ=255)

                                                        tATGCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAgtggtg                                 <  1:3778923/93‑1 (MQ=255)

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

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

                                                           gCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCtc                              >  2:5866043/1‑92 (MQ=255)

                                                           gCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTa                              <  2:536299/93‑1 (MQ=255)

                                                           gCGGAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTa                              <  2:2443975/93‑1 (MQ=255)

                                                             ggAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTAcc                            <  1:5698811/93‑1 (MQ=255)

                                                              gAAGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTAccc                           <  2:5173853/93‑1 (MQ=255)

                                                               aaGATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCcg                          <  1:5946555/93‑1 (MQ=255)

                                                                 gATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCcgcg                        >  1:3942607/1‑93 (MQ=255)

                                                                 gATAACGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCcgcg                        >  1:254513/1‑93 (MQ=255)

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

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

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

                                                                    aaCGTGGTGGTATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCa                     >  2:4061139/1‑93 (MQ=255)

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

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

                                                                            ggtATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTtct             <  1:6341418/93‑1 (MQ=255)

                                                                            ggtATTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAAGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTtct             >  2:6747842/1‑93 (MQ=255)

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

                                                                               aTTTGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCaa          <  1:3896903/93‑1 (MQ=255)

                                                                                tttGTCCGCGCGCTGTAACAGCAACGTGTACTGGTGGCAATCAGCCGAGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAAt         <  2:7025551/93‑1 (MQ=255)

                                                                                tttGTCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAAt         >  1:145174/1‑93 (MQ=255)

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

                                                                                    tCCGCGCGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCg     >  1:5496183/1‑93 (MQ=255)

                                                                                       gcgcgcTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGTgc  >  2:786715/1‑93 (MQ=255)

TGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAACGTGGTGGTATTTGTCCGCGTGCTGAATCAGCAACGTGTACTGGTGGCAATCAACCGTGGCGAGGCCTGTGAAGTGGTGCTACCCGCGTCACCGTTTCTCAATGCCGTGC  >  NC_000913/423291‑423470

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: `‑`