BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,562,936	(C)_5→6	coding (846/2424 nt)	b1489 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,562,931	1	.	C	84.2%	675.1 / 124.6	247	coding (851/2424 nt)	b1489	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (19/20); new base C (104/104); total (123/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 284 total aligned reads displayed.

ATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCG  >  NC_000913/1562839‑1563021

 tGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑c                                                                                           <  1:1014557/93‑1 (MQ=255)

 tGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑c                                                                                           >  1:2135513/1‑93 (MQ=255)

  gACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑cc                                                                                          >  2:8403748/1‑93 (MQ=255)

   aCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑ccc                                                                                         >  2:7487621/1‑93 (MQ=255)

    cGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑cccc                                                                                        <  1:4623576/93‑1 (MQ=255)

     gCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑ccccc                                                                                       >  1:6543177/1‑93 (MQ=255)

      cccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑cccccc                                                                                      >  2:10035154/1‑92 (MQ=255)

      cccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑cccccc                                                                                      >  1:3194848/1‑92 (MQ=255)

       ccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑cccccca                                                                                     >  2:4828225/1‑91 (MQ=255)

        cAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑ccccccat                                                                                    <  1:9160885/93‑4 (MQ=255)

        cAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑ccccccat                                                                                    <  1:5683124/93‑4 (MQ=255)

         aGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑ccccccatt                                                                                   <  1:7897806/93‑5 (MQ=255)

           tgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTc                                                                                  <  1:7148571/93‑1 (MQ=255)

           tgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTc                                                                                  >  1:6311351/1‑93 (MQ=255)

             tATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCAt                                                                                <  2:9378649/93‑1 (MQ=255)

             tATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCAt                                                                                >  2:3058563/1‑93 (MQ=255)

             tATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCAt                                                                                <  2:3792601/93‑1 (MQ=255)

              aTCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATg                                                                               >  1:6177865/1‑93 (MQ=255)

              aTCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGa                                                                              >  2:512184/1‑93 (MQ=255)

               tCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGa                                                                              >  2:1292925/1‑93 (MQ=255)

                cGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGCTGTTACGACAACTGATTTCCCCCCACTTCATGaa                                                                             >  1:6455963/1‑93 (MQ=255)

                 ggCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGaaa                                                                            <  1:2517639/93‑1 (MQ=255)

                 ggCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGaaa                                                                            <  1:7077503/93‑1 (MQ=255)

                  gCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAAt                                                                           >  2:3175329/1‑93 (MQ=255)

                     tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGc                                                                        >  2:1229284/1‑93 (MQ=255)

                     tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGc                                                                        >  2:6400504/1‑93 (MQ=255)

                     tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGc                                                                        >  2:5863297/1‑93 (MQ=255)

                     tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCg                                                                       <  2:3477613/93‑1 (MQ=255)

                     tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCg                                                                       >  1:8347628/1‑93 (MQ=255)

                      cccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCg                                                                       <  1:643124/93‑1 (MQ=255)

                       ccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCgg                                                                      >  2:7796300/1‑93 (MQ=255)

                        cGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCggg                                                                     <  2:252583/93‑1 (MQ=255)

                         gTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGgc                                                                    <  1:9280585/93‑1 (MQ=255)

                          ttGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGgct                                                                   <  1:1394734/93‑1 (MQ=255)

                           tGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGgctg                                                                  <  2:9580308/93‑1 (MQ=255)

                                 gTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGcac                                                            >  1:5497934/1‑93 (MQ=255)

                                  tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAACTGGCGGGCTGCTGccca                                                           >  2:2606729/1‑93 (MQ=255)

                                  tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGGTGCTGcaca                                                           >  2:2337077/1‑93 (MQ=255)

                                  tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGcaca                                                           >  1:9339607/1‑93 (MQ=255)

                                  tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACAt                                                          <  1:5116833/93‑1 (MQ=255)

                                   gCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACAt                                                          <  1:8627655/93‑1 (MQ=255)

                                     gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATgg                                                        >  2:2244908/1‑93 (MQ=255)

                                     gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATgg                                                        <  1:4098715/93‑1 (MQ=255)

                                     gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATgg                                                        <  1:2432333/93‑1 (MQ=255)

                                      aaCAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcg                                                      <  1:3560812/93‑1 (MQ=255)

                                        cAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcgg                                                     <  1:8793043/93‑1 (MQ=255)

                                         aGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcggc                                                    <  1:4338863/93‑1 (MQ=255)

                                          gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACACGgcggcg                                                   >  1:6908565/1‑93 (MQ=255)

                                          gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga                                                  <  2:5957030/93‑1 (MQ=255)

                                          gCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga                                                  <  2:615441/93‑1 (MQ=255)

                                           cGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCga                                                  <  1:9886204/93‑1 (MQ=255)

                                            gAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCgag                                                 >  1:5051574/1‑93 (MQ=255)

                                              agcATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCgagca                                               >  2:115451/3‑91 (MQ=255)

                                              agcATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGaa                                               >  1:1276327/3‑93 (MQ=255)

                                               gcATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAat                                              >  2:3016845/2‑93 (MQ=255)

                                                cATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAata                                             <  1:7446729/93‑1 (MQ=255)

                                                cATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAatat                                            <  1:6789137/93‑1 (MQ=255)

                                                cATCCGATGCGTAGAGTACAGATCCGATGTTACGACAAAGGATTTCCCCCCTTTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAata                                             <  1:1095270/93‑1 (MQ=255)

                                                  tGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGGCGCTTCACATAGGGGCGAGAAGAtt                                           >  2:7179198/1‑93 (MQ=255)

                                                   gCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTg                                          <  2:4282202/93‑1 (MQ=255)

                                                    cGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGGataattc                                         >  1:9645450/1‑89 (MQ=255)

                                                        tCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTc                                     <  2:1599728/93‑1 (MQ=255)

                                                        tCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTc                                     >  1:7782342/1‑93 (MQ=255)

                                                        tCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTc                                     <  1:2753880/93‑1 (MQ=255)

                                                          gTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGGCAATATTGGCTTcca                                   >  1:4887296/1‑91 (MQ=255)

                                                              aGTACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTg                               <  2:8759357/93‑1 (MQ=255)

                                                                tACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGc                             <  1:7741235/93‑1 (MQ=255)

                                                                tACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGc                             <  1:7512119/93‑1 (MQ=255)

                                                                tACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGc                             <  1:3890702/93‑1 (MQ=255)

                                                                tACAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGc                             >  2:5596651/1‑93 (MQ=255)

                                                                 aCAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCg                            <  2:8120514/93‑1 (MQ=255)

                                                                  cAGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGa                           <  2:6095399/93‑1 (MQ=255)

                                                                   aGATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGaa                          >  2:6993691/1‑93 (MQ=255)

                                                                    gATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAAt                         <  1:6234494/93‑1 (MQ=255)

                                                                    gATTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGATGGCGAAt                         >  2:4068192/1‑93 (MQ=255)

                                                                     aTTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATc                        <  1:2898186/93‑1 (MQ=255)

                                                                     aTTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATc                        <  2:9714420/93‑1 (MQ=255)

                                                                     aTTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATc                        <  2:1729988/93‑1 (MQ=255)

                                                                     aTTCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCAAATc                        >  1:1649429/1‑93 (MQ=255)

                                                                      ttCGATGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCt                       >  1:5277733/1‑93 (MQ=255)

                                                                           tGTTACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATTGCGGCGAGGATATTGCCTTCAAGCTAGCCAAACTGTCGt                  >  1:8251252/1‑93 (MQ=255)

                                                                             ttACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttc                <  1:4916246/93‑1 (MQ=255)

                                                                              tACGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttct               <  2:4724979/93‑1 (MQ=255)

                                                                               aCGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttctt              >  2:398950/1‑93 (MQ=255)

                                                                               aCGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttctt              <  1:3085714/93‑1 (MQ=255)

                                                                               aCGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttctt              <  2:9025786/93‑1 (MQ=255)

                                                                                cGACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGttcttc             <  2:4389637/93‑1 (MQ=255)

                                                                                aggcAAATGATTT‑TCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCg            <  1:1483560/90‑1 (MQ=255)

                                                                                 gACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCg            >  1:4713914/1‑93 (MQ=255)

                                                                                 gACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCg            <  2:7077792/93‑1 (MQ=255)

                                                                                 gACAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCg            <  1:7292957/93‑1 (MQ=255)

                                                                                  aCAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCgg           >  2:5743956/1‑93 (MQ=255)

                                                                                   cAAATGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGt          >  1:3703263/1‑93 (MQ=255)

                                                                                   cAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTg         >  2:7317816/1‑93 (MQ=255)

                                                                                    aaaTGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTg         >  2:7536899/1‑93 (MQ=255)

                                                                                     aaTGATTTCCCCCCATTTCATGAAATGGTGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGa        <  2:7605446/93‑1 (MQ=255)

                                                                                     aaTGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGa        <  1:5359914/93‑1 (MQ=255)

                                                                                      aTGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGat       <  1:2054043/93‑1 (MQ=255)

                                                                                      aTGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGat       <  2:2865772/93‑1 (MQ=255)

                                                                                       tGATTTCCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGata      <  2:1481775/93‑1 (MQ=255)

ATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTT‑CCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGAATATTGCCTTCAAGCTGGCGAATCTGCCGTTCTTCGGTGATATCCG  >  NC_000913/1562839‑1563021

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