BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,586,923	C→T	Q437* (CAG→TAG)	acrD →	AcrD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,586,923	0	C	T	85.4%	763.9 / 86.7	296	Q437* (CAG→TAG)	acrD	AcrD
Reads supporting (aligned to +/- strand): ref base C (25/18); new base T (109/143); total (135/161)
Fisher's exact test for biased strand distribution p-value = 9.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-02

Only 100 of 296 total aligned reads displayed.

GATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTA  >  NC_000913/2586833‑2587014

 aTGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGa                                                                                          >  2:1364496/1‑93 (MQ=255)

 aTGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGa                                                                                          >  1:3208178/1‑93 (MQ=255)

     cGCCGTCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                      <  2:5081089/93‑1 (MQ=255)

     cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCa                                                                                      <  2:6740057/93‑1 (MQ=255)

     cGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCa                                                                                      >  2:446734/1‑93 (MQ=255)

      gCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAg                                                                                     <  1:3006250/93‑1 (MQ=255)

       ccATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAgg                                                                                    >  1:7452894/1‑93 (MQ=255)

        cATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAggg                                                                                   >  2:4849053/1‑93 (MQ=255)

         aTCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAgggg                                                                                  <  1:1242101/93‑1 (MQ=255)

         aTCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAgggg                                                                                  >  2:2482854/1‑93 (MQ=255)

         aTCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAggga                                                                                  <  1:1509680/93‑2 (MQ=255)

            gtggtggtgGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCa                                                                               <  2:2909529/93‑1 (MQ=255)

            gtggtggtgGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCa                                                                               >  1:3828210/1‑93 (MQ=255)

             tggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCAc                                                                              >  2:1947589/1‑93 (MQ=255)

              ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACt                                                                             >  1:2459854/1‑93 (MQ=255)

              ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACt                                                                             <  2:1796110/93‑1 (MQ=255)

              ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACt                                                                             <  1:882325/93‑1 (MQ=255)

              ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACt                                                                             >  1:6652261/1‑93 (MQ=255)

              ggtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACt                                                                             >  1:4987940/1‑93 (MQ=255)

               gtggtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTg                                                                            >  1:5823675/1‑93 (MQ=255)

                  gtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTc                                                                         <  2:3080525/93‑1 (MQ=255)

                  gtggAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTc                                                                         <  2:5370166/93‑1 (MQ=255)

                     gAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCggg                                                                      >  1:6722082/1‑93 (MQ=255)

                       aaaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGAt                                                                    >  1:6071152/1‑93 (MQ=255)

                        aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGAtt                                                                   >  2:524478/1‑93 (MQ=255)

                        aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGAtt                                                                   >  1:1759984/1‑93 (MQ=255)

                        aaCGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGAtt                                                                   >  2:6786245/1‑93 (MQ=255)

                          cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGc                                                                 <  2:1193502/93‑1 (MQ=255)

                          cGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGc                                                                 <  2:1695721/93‑1 (MQ=255)

                           gTCGAACGTATTATGAGTGAGGAAGGACTCGCTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCg                                                                <  2:6287039/93‑1 (MQ=255)

                           gTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCg                                                                >  2:4295082/1‑93 (MQ=255)

                           gTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCg                                                                >  2:1022101/1‑93 (MQ=255)

                            tCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                               >  2:1241237/1‑93 (MQ=255)

                            tCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                               <  2:4048393/93‑1 (MQ=255)

                            tCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                               <  2:4670315/93‑1 (MQ=255)

                            tCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGa                                                               <  1:5646358/93‑1 (MQ=255)

                             cGAACGTATTATGAGTGAGGAATGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGAt                                                              <  1:5332020/93‑1 (MQ=255)

                             cGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGAt                                                              >  1:4371994/1‑93 (MQ=255)

                             cGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGAt                                                              >  1:3395620/1‑93 (MQ=255)

                              gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATg                                                             <  1:6467854/93‑1 (MQ=255)

                              gAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATg                                                             >  1:227847/1‑93 (MQ=255)

                               aaCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATgg                                                            >  2:5140977/1‑93 (MQ=255)

                                aCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGt                                                           <  1:5156172/93‑1 (MQ=255)

                                aCGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGt                                                           <  1:7160805/93‑1 (MQ=255)

                                   tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCt                                                        <  2:7459108/93‑1 (MQ=255)

                                   tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCt                                                        >  1:5251495/1‑93 (MQ=255)

                                   tattatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCt                                                        >  2:7026884/1‑93 (MQ=255)

                                    attatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCtt                                                       <  1:5971680/93‑1 (MQ=255)

                                     ttatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCttt                                                      >  2:4390139/1‑93 (MQ=255)

                                      tatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTc                                                     >  1:3474195/1‑93 (MQ=255)

                                      tatGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTc                                                     >  1:4118317/1‑93 (MQ=255)

                                        tgagtgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTcgg                                                   >  2:6203038/1‑93 (MQ=255)

                                            tgagGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGt                                               >  2:3848393/1‑93 (MQ=255)

                                              aggaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAt                                             <  1:769247/93‑1 (MQ=255)

                                                gaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAttt                                           <  2:3508680/93‑1 (MQ=255)

                                                gaaggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTAttt                                           <  1:4693353/93‑1 (MQ=255)

                                                  aggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGt                                         >  1:4941684/1‑93 (MQ=255)

                                                  aggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGt                                         <  1:4483905/93‑1 (MQ=255)

                                                   ggtctctctcCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATTGTTCTTTCGGCGGTATTTGTa                                        <  2:1596281/90‑1 (MQ=255)

                                                   ggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                        >  1:3677255/1‑93 (MQ=255)

                                                   ggaCTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTa                                        >  2:369336/1‑93 (MQ=255)

                                                     actcactcCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTAcc                                      >  2:656221/1‑93 (MQ=255)

                                                       tcactcCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCaa                                    >  1:3557535/1‑93 (MQ=255)

                                                       tcactcCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCaa                                    <  2:2237721/93‑1 (MQ=255)

                                                         actcCTCGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATg                                  <  1:2876649/93‑1 (MQ=255)

                                                          ctcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTGCCAATgg                                 <  2:3586178/93‑1 (MQ=255)

                                                          ctcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                 >  2:1464075/1‑93 (MQ=255)

                                                          ctcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                 >  1:2951217/1‑93 (MQ=255)

                                                          ctcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                 >  1:102339/1‑93 (MQ=255)

                                                          ctcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATgg                                 >  2:6100087/1‑93 (MQ=255)

                                                           tcctcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGc                                >  2:5887758/1‑93 (MQ=255)

                                                              tcGCGAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                             <  2:2606992/93‑1 (MQ=255)

                                                              tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                             >  2:935055/1‑93 (MQ=255)

                                                              tcGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCctt                             >  1:5590377/1‑93 (MQ=255)

                                                               cgcgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttc                            <  1:7380680/93‑1 (MQ=255)

                                                                gcgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttct                           <  2:6146353/93‑1 (MQ=255)

                                                                gcgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttct                           <  1:1978568/93‑1 (MQ=255)

                                                                 gnaaGCCACACGT‑AATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCttcttc                         <  2:6932426/91‑1 (MQ=255)

                                                                 cgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttctt                          <  2:3458634/93‑1 (MQ=255)

                                                                 cgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttctt                          <  2:7160721/93‑1 (MQ=255)

                                                                 cgAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttctt                          <  1:4498208/93‑1 (MQ=255)

                                                                 cgAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCcttctt                          <  2:338680/93‑1 (MQ=255)

                                                                  gAAGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCCGCGGTATTTGTACCAATGGCCttcttc                         <  2:6423780/93‑1 (MQ=255)

                                                                    aGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcgg                       >  2:5583757/1‑93 (MQ=255)

                                                                    aGCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTcgg                       >  1:1424279/1‑93 (MQ=255)

                                                                     gCCACACGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGTTTGCGATGGGTCTTTCGGCGGTGTTTGTACCAATGGCCTTCTTcggc                      >  1:5004400/1‑93 (MQ=255)

                                                                        acacGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCggcggc                   >  2:757485/1‑93 (MQ=255)

                                                                         cacGTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGca                  >  2:7261064/1‑93 (MQ=255)

                                                                            gTAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcacca               <  2:4526831/93‑1 (MQ=255)

                                                                             tAAATCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGcaccac              <  2:2447582/93‑1 (MQ=255)

                                                                                aTCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCgg           <  1:6383148/93‑1 (MQ=255)

                                                                                 tCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt          <  1:6886518/93‑1 (MQ=255)

                                                                                 tCGATGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGt          <  1:3148843/93‑1 (MQ=255)

                                                                                    aTGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGcc       >  2:1989460/1‑93 (MQ=255)

                                                                                     tGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCa      <  1:3935924/93‑1 (MQ=255)

                                                                                     tGGGGTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCa      >  1:1600491/1‑93 (MQ=255)

                                                                                       gggTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc    <  1:3312285/93‑1 (MQ=255)

                                                                                       gggTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc    >  1:2763336/1‑93 (MQ=255)

                                                                                       gggTAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATc    <  2:5401262/93‑1 (MQ=255)

                                                                                         gCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTa  <  1:2541763/93‑1 (MQ=255)

GATGACGCCATCGTGGTGGTGGAAAACGTCGAACGTATTATGAGTGAGGAAGGACTCACTCCTCGCGAAGCCACACGTAAATCGATGGGGCAGATCCAGGGGGCACTGGTCGGGATTGCGATGGTTCTTTCGGCGGTATTTGTACCAATGGCCTTCTTCGGCGGCACCACCGGTGCCATCTA  >  NC_000913/2586833‑2587014

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