BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,858,325	A→G	D344D (GAT→GAC)	b1776 ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,858,325	0	A	G	99.3%	430.1 / ‑5.3	139	D344D (GAT→GAC)	b1776	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (83/55); total (84/55)
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.96e-01

Only 100 of 141 total aligned reads displayed.

GTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTATCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCT  >  NC_000913/1858235‑1858412

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                       >  1:87381/1‑93 (MQ=255)

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                       <  2:4554480/93‑1 (MQ=255)

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                       <  2:4185102/93‑1 (MQ=255)

gTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTc                                                                                       >  1:1213840/1‑93 (MQ=255)

  tttGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCgc                                                                                     >  1:2548596/1‑92 (MQ=255)

   ttGTTATCTGTTCCATGCGATAACGTTCCTTCCGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGAc                                                                                    <  1:508830/93‑1 (MQ=255)

     gTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCa                                                                                  >  1:2925511/1‑93 (MQ=255)

       tATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATc                                                                                <  1:3279295/93‑1 (MQ=255)

       tATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATc                                                                                <  2:1898627/93‑1 (MQ=255)

       tATCTGTTCCATGCAATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGG‑TTTGTCGACCATCg                                                                               >  2:2818777/1‑93 (MQ=255)

        aTCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCg                                                                               >  1:3996815/1‑93 (MQ=255)

         tCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGc                                                                              <  2:131203/93‑1 (MQ=255)

         tCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGc                                                                              >  2:3285851/1‑93 (MQ=255)

         tCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGCTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGc                                                                              >  1:2440807/1‑93 (MQ=255)

          cTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCa                                                                             >  2:1805501/1‑93 (MQ=255)

           tGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAt                                                                            <  1:828761/93‑1 (MQ=255)

           tGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAt                                                                            <  1:3708375/93‑1 (MQ=255)

            gTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATc                                                                           <  1:4572352/93‑1 (MQ=255)

            gTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATc                                                                           >  1:3680943/1‑93 (MQ=255)

             ttCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCAtta                                                                          >  2:4154779/1‑91 (MQ=255)

             ttCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCa                                                                          >  1:3060270/1‑93 (MQ=255)

               ccATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCaaa                                                                        >  1:3703756/1‑93 (MQ=255)

                   gCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCCTCAAACCCt                                                                    <  1:1004030/93‑1 (MQ=255)

                   gCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCt                                                                    <  1:447644/93‑1 (MQ=255)

                    cGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCtt                                                                   <  1:3111420/93‑1 (MQ=255)

                     gATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTc                                                                  <  2:3322812/93‑1 (MQ=255)

                     caTAACTTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTc                                                                  <  2:1154230/92‑1 (MQ=255)

                      aTAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTcg                                                                 <  2:901190/93‑1 (MQ=255)

                       tAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTcgc                                                                >  2:3812320/1‑93 (MQ=255)

                       tAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTcgc                                                                >  2:3910794/1‑93 (MQ=255)

                         aCGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCgcgc                                                              >  2:2091098/1‑93 (MQ=255)

                           gTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCa                                                            >  2:4100428/1‑93 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                           >  1:683482/1‑93 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                           >  1:856214/1‑93 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                           <  2:3482006/93‑1 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                           >  1:4307323/1‑93 (MQ=255)

                            ttCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAt                                                           <  1:1202131/93‑1 (MQ=255)

                             tCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCAtt                                                          <  1:2184850/93‑1 (MQ=255)

                               cTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGc                                                        <  2:1164956/93‑1 (MQ=255)

                               cTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGc                                                        >  2:4411363/1‑93 (MQ=255)

                               cTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGc                                                        <  2:2593287/93‑1 (MQ=255)

                               cTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGc                                                        <  1:2694904/93‑1 (MQ=255)

                                 tCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGa                                                      >  1:3666448/1‑93 (MQ=255)

                                   aGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACa                                                    <  2:4588570/93‑1 (MQ=255)

                                    ggAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAg                                                   >  2:888698/1‑93 (MQ=255)

                                    ggAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAg                                                   >  1:132015/1‑93 (MQ=255)

                                    ggAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTGTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAg                                                   <  1:3414044/93‑1 (MQ=255)

                                     gagaGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAgg                                                  >  1:1081487/1‑93 (MQ=255)

                                     gagaGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAgg                                                  >  1:1176765/1‑93 (MQ=255)

                                       gagTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGcc                                                >  2:4559120/1‑93 (MQ=255)

                                        agTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCg                                               >  1:1049749/1‑93 (MQ=255)

                                         gTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGa                                              <  2:2282256/93‑1 (MQ=255)

                                          tGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATGAAACACTTCGCGGCATTGCGACAGGCCGAt                                             >  2:2435871/1‑93 (MQ=255)

                                            gTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATAc                                           >  1:3294436/1‑93 (MQ=255)

                                              tATTCATCAGAGACGTAAGTCATGATCACTTAGATGGCGGTTTTGTCGACCATCACCTCAAACCCTTCGCGCCAGTGCGACAGGCCGATACGa                                         <  2:1363033/93‑1 (MQ=255)

                                              tATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGa                                         <  1:472665/93‑1 (MQ=255)

                                               aTTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGAt                                        <  2:22242/93‑1 (MQ=255)

                                                ttCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATg                                       <  2:2427213/93‑1 (MQ=255)

                                                ttCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATg                                       >  1:1617215/1‑93 (MQ=255)

                                                 tcatcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGc                                      <  1:2545886/93‑1 (MQ=255)

                                                 tcatcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGc                                      >  1:2790719/1‑93 (MQ=255)

                                                  catcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCg                                     >  1:1135017/1‑93 (MQ=255)

                                                   atcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGt                                    >  2:3023455/1‑93 (MQ=255)

                                                    tcaAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTg                                   >  1:4095428/1‑93 (MQ=255)

                                                     caAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                  >  2:2147643/1‑93 (MQ=255)

                                                     caAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                  >  1:396695/1‑93 (MQ=255)

                                                     caAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGa                                  <  2:3160633/93‑1 (MQ=255)

                                                        aGTCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatca                               <  2:607597/93‑1 (MQ=255)

                                                          tCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatcatc                             <  1:324474/93‑1 (MQ=255)

                                                          tCGTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGatcatc                             <  1:3668437/93‑1 (MQ=255)

                                                            gTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCgg                           >  2:1153909/1‑93 (MQ=255)

                                                            gTAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATAGCATCAAACCCTTCGCGCCATTGCGACAGGACGATACGATGCGTGATCATCgg                           >  2:1028016/1‑93 (MQ=255)

                                                             tAAGTCATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAACCCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGt                          <  2:4019438/93‑1 (MQ=255)

                                                                  cATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGAc                     >  2:574472/1‑93 (MQ=255)

                                                                  cATGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGAc                     <  1:4434214/93‑1 (MQ=255)

                                                                   aTGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACt                    <  1:2576574/93‑1 (MQ=255)

                                                                   aTGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACt                    >  2:2212510/1‑93 (MQ=255)

                                                                   aTGATCACTTTGATTGCGGTTTTATCGACCATCGCATCAAACCCTTCCCGCCATTGCGACAGGTCGATACGATGCGTGATCATCGGTGTGACt                    >  2:1141577/1‑93 (MQ=255)

                                                                    tGATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACtt                   >  2:1478018/1‑93 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  <  1:4657945/93‑1 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  >  2:4391829/1‑93 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  <  2:3107277/93‑1 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  <  1:2834855/93‑1 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  >  1:2366727/1‑93 (MQ=255)

                                                                     gATCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCCCCATTGCGAAAGGCCGATACGATGCGTGATCATCGGTTTGACttt                  <  2:2367219/93‑1 (MQ=255)

                                                                      aTCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTg                 <  1:1678163/93‑1 (MQ=255)

                                                                       tCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAGCCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGa                <  1:1384648/93‑1 (MQ=255)

                                                                       tCACTTTGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGa                >  2:3758670/1‑93 (MQ=255)

                                                                           tttGATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAgc            >  2:3286115/1‑93 (MQ=255)

                                                                           tttGATTGCGGTTTTGTCGAACATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATatc            >  1:1130729/1‑91 (MQ=255)

                                                                              gATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgcc         >  1:2163844/1‑93 (MQ=255)

                                                                              gATTGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgcc         <  2:1028859/93‑1 (MQ=255)

                                                                                ttGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGcgccgc       >  2:2335341/1‑93 (MQ=255)

                                                                                 tGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGGCGCt      >  1:2494400/1‑93 (MQ=255)

                                                                                 tGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCt      >  1:1620337/1‑93 (MQ=255)

                                                                                 tGCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCt      <  1:3868468/93‑1 (MQ=255)

                                                                                  gCGGTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTg     >  2:2184451/1‑93 (MQ=255)

                                                                                    ggTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGc   <  1:3611317/93‑1 (MQ=255)

                                                                                    ggTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATAGGTTTGACTTTGATAGCGCCGCTGGc   >  2:3795052/1‑93 (MQ=255)

                                                                                     gTTTTGTCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCt  >  2:1797874/1‑93 (MQ=255)

GTTTTGTTATCTGTTCCATGCGATAACGTTCCTTCAGGAGAGTGGTTATTCATCAAAGTCGTAAGTCATGATCACTTTGATTGCGGTTTTATCGACCATCGCATCAAACCCTTCGCGCCATTGCGACAGGCCGATACGATGCGTGATCATCGGTTTGACTTTGATAGCGCCGCTGGCT  >  NC_000913/1858235‑1858412

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