BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	591,675	T→C	S243G (AGC→GGC)	nfrB ←	NfrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	591,675	0	T	C	97.7%	500.0 / ‑0.0	171	S243G (AGC→GGC)	nfrB	NfrB
Reads supporting (aligned to +/- strand): ref base T (3/1); new base C (84/82); total (87/84)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 182 total aligned reads displayed.

GAAGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAACTCATC  >  NC_000913/591586‑591765

gAAGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGaa                                                                                         >  2:53772/1‑93 (MQ=255)

 aaGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGaaa                                                                                        <  1:642638/93‑1 (MQ=255)

 aaGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGaaa                                                                                        <  1:323724/93‑1 (MQ=255)

  aGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAAc                                                                                       >  2:3028667/1‑93 (MQ=255)

  aGCCAATGTCGTAATCTTCAGGAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAAc                                                                                       >  2:4470581/1‑93 (MQ=255)

   gCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACa                                                                                      >  2:4225018/1‑93 (MQ=255)

    ccAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGACACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAg                                                                                     >  1:2937370/1‑93 (MQ=255)

     cAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAgg                                                                                    >  1:4261212/1‑93 (MQ=255)

      aaTGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAACCAggg                                                                                   >  1:1359802/1‑92 (MQ=255)

      aaTGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAggg                                                                                   >  1:3480154/1‑92 (MQ=255)

       aTGTCGTAATCTTCAGTAAGCCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTg                                                                                  <  2:2865732/93‑1 (MQ=255)

        tGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGc                                                                                 >  2:4178180/1‑93 (MQ=255)

         gTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGcc                                                                                >  1:330646/1‑93 (MQ=255)

          tCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCg                                                                               >  1:3170691/1‑93 (MQ=255)

           cctAATCTTGAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGa                                                                              <  1:1421595/91‑1 (MQ=255)

           cGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGa                                                                              <  1:1347924/93‑1 (MQ=255)

             tAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCCCGGCGCGGCGGCCGAAACAGGTGCCGACg                                                                            >  1:2289849/1‑93 (MQ=255)

             tAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACg                                                                            >  2:3788688/1‑93 (MQ=255)

              aaTCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGc                                                                           >  1:389505/1‑93 (MQ=255)

               aTCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGcc                                                                          >  1:1148095/1‑93 (MQ=255)

                     aGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCt                                                                    >  1:364221/1‑93 (MQ=255)

                      gTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTg                                                                   >  1:595383/1‑93 (MQ=255)

                      gTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTg                                                                   <  1:2415519/93‑1 (MQ=255)

                       tAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTgg                                                                  <  1:4155321/93‑1 (MQ=255)

                       tAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTgg                                                                  >  1:1574347/1‑93 (MQ=255)

                       tAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTgg                                                                  <  2:666804/93‑1 (MQ=255)

                       tAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTgg                                                                  >  1:1782523/1‑93 (MQ=255)

                         aGACTCTTCACGTCGGCAGCAATACCCTCACCGGCAGCTAACAGTGCGGTCACGGCGCGGCCGGCGCAAAACGCGCCGACGCCCGCGCCGGGc                                                                >  2:3534346/1‑93 (MQ=255)

                         aGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGc                                                                >  2:77035/1‑93 (MQ=255)

                         aGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGc                                                                <  1:4364771/93‑1 (MQ=255)

                         aGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGc                                                                >  1:1560263/1‑93 (MQ=255)

                          gACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGGCGACGCCTGCGCTGGGCa                                                               >  2:917659/1‑93 (MQ=255)

                          gACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCAGAAACAGGTGCCGCCGCCTGCGCCgggac                                                               >  1:1324119/1‑91 (MQ=255)

                           aCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCAc                                                              >  1:3067509/1‑93 (MQ=255)

                           aCTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCAc                                                              <  2:2683017/93‑1 (MQ=255)

                            ctctGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACt                                                             >  2:2581417/1‑93 (MQ=255)

                             tctGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACtt                                                            <  1:4167598/93‑1 (MQ=255)

                               tGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAATAGGTGCCGACGTCTGCGCTGGGCACTTGt                                                          >  2:4681741/1‑93 (MQ=255)

                               tGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGt                                                          <  1:1305729/93‑1 (MQ=255)

                                 cACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTcc                                                        <  1:1691602/93‑1 (MQ=255)

                                  aCGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCg                                                       <  2:1801225/93‑1 (MQ=255)

                                  aCGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCg                                                       >  1:525662/1‑93 (MQ=255)

                                  aCGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCg                                                       >  1:165730/1‑93 (MQ=255)

                                  aCGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGAGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCg                                                       >  2:1940751/1‑93 (MQ=255)

                                   cGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCgg                                                      <  2:1619513/93‑1 (MQ=255)

                                   cGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCgg                                                      >  2:3868927/1‑93 (MQ=255)

                                    gTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGc                                                     >  2:2889933/1‑93 (MQ=255)

                                    gTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGc                                                     <  2:3693291/93‑1 (MQ=255)

                                       gAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAg                                                  <  2:1061814/93‑1 (MQ=255)

                                        aaaGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAgg                                                 <  1:1613250/93‑1 (MQ=255)

                                         aaGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAggg                                                >  2:2387754/1‑93 (MQ=255)

                                         aaGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAggg                                                >  2:105771/1‑93 (MQ=255)

                                           gCAATACCGTCACCGTCAGCTAACAGTGCGGCCACGCGGCGGCGGCCGAAACAGGTGCCGACCCCTGCGCTGGGCACTTGTCCGGCGAGGGCt                                              <  1:1860968/93‑1 (MQ=255)

                                            cAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCtt                                             <  1:924696/93‑1 (MQ=255)

                                             aaTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTc                                            <  1:1826535/93‑1 (MQ=255)

                                             aaTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTc                                            <  2:376075/93‑1 (MQ=255)

                                             aaTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTc                                            <  1:3829488/93‑1 (MQ=255)

                                             aaTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTc                                            >  2:1882395/1‑93 (MQ=255)

                                              aTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTcc                                           >  2:2235701/1‑92 (MQ=255)

                                              aTACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCa                                           >  1:3972877/1‑93 (MQ=255)

                                               tACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCAc                                          <  2:1275501/93‑1 (MQ=255)

                                               tACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCAc                                          >  1:1126136/1‑93 (MQ=255)

                                                accgtcaccgtcAGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACg                                         >  1:4047409/1‑93 (MQ=255)

                                                   gtcaccgtcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCAc                                      >  1:2975481/1‑93 (MQ=255)

                                                    tcaccgtcaGCTCACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCAcc                                     >  2:4164131/1‑93 (MQ=255)

                                                    tcaccgtcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCAcc                                     <  2:3723511/93‑1 (MQ=255)

                                                     caccgtcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCg                                    <  1:1931180/93‑1 (MQ=255)

                                                       ccgtcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGa                                  >  1:3485144/1‑93 (MQ=255)

                                                        cgtcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGaa                                 <  2:1391262/93‑1 (MQ=255)

                                                         gtcaGCTAACAGTGCGGTCACGGCGCGGCGGCTGACACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCCGCGAGGGCTTCACGCCCCGGCAc                                >  2:2315141/1‑93 (MQ=255)

                                                          tcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACa                               <  2:1374827/93‑1 (MQ=255)

                                                          tcaGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACa                               <  2:4096971/93‑1 (MQ=255)

                                                            aGCTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATc                             >  1:3482631/1‑93 (MQ=255)

                                                             gCTAACAGTGCGGTCACGCCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCt                            <  2:1289472/93‑1 (MQ=255)

                                                              cTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCtt                           >  1:1791252/1‑93 (MQ=255)

                                                              cTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCtt                           >  2:2923827/1‑93 (MQ=255)

                                                              cTAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCtt                           <  1:1065315/93‑1 (MQ=255)

                                                              cTAACAGTGCGGGCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCtt                           <  2:1989319/93‑1 (MQ=255)

                                                               tAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCttt                          <  1:2786305/93‑1 (MQ=255)

                                                               tAACAGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGAGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACCTCttt                          >  2:3041223/1‑93 (MQ=255)

                                                                   aGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCa                      <  1:4574885/93‑1 (MQ=255)

                                                                   aGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCa                      <  2:1529135/93‑1 (MQ=255)

                                                                   aGTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCa                      >  1:3312548/1‑93 (MQ=255)

                                                                    gTGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCAt                     >  2:2474755/1‑93 (MQ=255)

                                                                     tGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCAtc                    >  1:2087604/1‑92 (MQ=255)

                                                                     tGCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATg                    <  1:4651874/93‑1 (MQ=255)

                                                                      gCGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGc                   <  1:2260233/93‑1 (MQ=255)

                                                                       cGGTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCa                  <  1:4169904/93‑1 (MQ=255)

                                                                         gTCACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGc                <  2:3294815/93‑1 (MQ=255)

                                                                          gcACGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGct               <  2:4440748/92‑1 (MQ=255)

                                                                            aCGGCGCGGCGGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGctct             >  1:2016470/1‑93 (MQ=255)

                                                                             cggcgcggcgGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTg            >  1:2227355/1‑93 (MQ=255)

                                                                               gcgcggcgGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGaa          >  2:829118/1‑93 (MQ=255)

                                                                                cgcggcgGCCGAAACAGGTGCCGACTCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGaaa         <  2:4595816/93‑1 (MQ=255)

                                                                                 gcggcgGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGaaaa        <  1:2349356/93‑1 (MQ=255)

                                                                                 gcggcgGCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGaaaa        >  2:1478419/1‑93 (MQ=255)

                                                                                  cggcggCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAAc       <  1:202399/93‑1 (MQ=255)

                                                                                  cggcggCCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAAc       <  1:753869/93‑1 (MQ=255)

                                                                                     cggcCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAACtca    <  2:3443118/93‑1 (MQ=255)

                                                                                      ggcCGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAACtcat   >  1:3105191/1‑93 (MQ=255)

GAAGCCAATGTCGTAATCTTCAGTAAGACTCTGCACGTCGAAAGCAATACCGTCACCGTCAGCTAACAGTGCGGTCACGGCGCGGCGGCTGAAACAGGTGCCGACGCCTGCGCTGGGCACTTGTCCGGCGAGGGCTTCACGCACCGGAACATCTTTGCCATGCAGCTCTGAAAACTCATC  >  NC_000913/591586‑591765

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