BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,132,580	T→A	I389N (ATC→AAC)	katG →	catalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,132,580	0	T	A	98.4%	805.5 / ‑0.3	249	I389N (ATC→AAC)	katG	catalase
Reads supporting (aligned to +/- strand): ref base T (2/2); new base A (138/103); total (141/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 273 total aligned reads displayed.

TCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGA  >  NC_000913/4132490‑4132669

 cccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactc                                                                                        <  2:5005411/93‑1 (MQ=255)

  ccGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactct                                                                                       <  1:7264212/93‑1 (MQ=255)

   cGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAactctc                                                                                      <  2:7119767/93‑1 (MQ=255)

     gATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgt                                                                                    <  2:3745798/93‑1 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   >  1:7914435/1‑93 (MQ=255)

      aTCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtc                                                                                   <  1:4236432/93‑1 (MQ=255)

       tCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtcg                                                                                  >  1:104725/1‑93 (MQ=255)

       tCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCtcgtcg                                                                                  <  1:5036634/93‑1 (MQ=255)

         cGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGtt                                                                                >  1:931970/1‑93 (MQ=255)

         cGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGtt                                                                                <  2:3114183/93‑1 (MQ=255)

          gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGttt                                                                               <  1:7314170/93‑1 (MQ=255)

          gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGttt                                                                               <  1:7187571/93‑1 (MQ=255)

          gTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGttt                                                                               >  1:1868245/1‑93 (MQ=255)

            ttGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTcc                                                                             <  2:4073614/93‑1 (MQ=255)

              gATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTc                                                                           >  2:6596981/1‑93 (MQ=255)

              gATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTc                                                                           >  2:5478456/1‑93 (MQ=255)

              gATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTc                                                                           <  2:4778458/93‑1 (MQ=255)

              gATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTc                                                                           >  2:3980504/1‑93 (MQ=255)

               aTCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCa                                                                          >  1:6654736/1‑93 (MQ=255)

                 ccGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAAc                                                                        >  2:5551660/1‑93 (MQ=255)

                  cGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACg                                                                       <  1:5201467/93‑1 (MQ=255)

                   gTCGAAGAAACGTAAACCGACAATGCAGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGa                                                                      <  2:7899430/93‑1 (MQ=255)

                    tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGAt                                                                     <  1:8288187/93‑1 (MQ=255)

                    tCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGAt                                                                     <  2:8554430/93‑1 (MQ=255)

                     cGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATc                                                                    <  1:1706299/93‑1 (MQ=255)

                      gaagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                   >  2:273769/1‑93 (MQ=255)

                      gaagaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATcc                                                                   <  2:1936621/93‑1 (MQ=255)

                        agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGc                                                                 >  2:7205596/1‑93 (MQ=255)

                        agaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGc                                                                 >  1:7035596/1‑93 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGcg                                                                >  1:2937480/1‑92 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                >  1:4102132/1‑93 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                >  2:2418182/1‑93 (MQ=255)

                         gaaACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCa                                                                >  1:6233942/1‑93 (MQ=255)

                          aaaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAg                                                               >  1:6333999/1‑93 (MQ=255)

                          aaaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAg                                                               <  2:8425493/93‑1 (MQ=255)

                           aaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAgg                                                              >  2:6315429/1‑93 (MQ=255)

                           aaCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAgg                                                              >  2:2164013/1‑93 (MQ=255)

                            aCGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGc                                                             <  2:4325635/93‑1 (MQ=255)

                              ctAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGt                                                           <  2:1467707/92‑1 (MQ=255)

                               tAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGtt                                                          <  1:8705901/93‑1 (MQ=255)

                                aaaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATACTCAGGCGTTc                                                         >  2:8235655/1‑93 (MQ=255)

                                 aaCCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCa                                                        >  2:919085/1‑93 (MQ=255)

                                    cGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGTTCCTCAGGCGCTCAACg                                                     >  1:6053019/1‑93 (MQ=255)

                                     gACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGa                                                    <  2:3580032/93‑1 (MQ=255)

                                        aaTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGc                                                 >  1:3506652/1‑93 (MQ=255)

                                        aaTGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGc                                                 >  1:3860991/1‑93 (MQ=255)

                                          tGCTGTTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCt                                               <  2:3374599/93‑1 (MQ=255)

                                          tGCTGGTGACCGACCTGGCGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCt                                               <  1:5327405/93‑1 (MQ=255)

                                           gCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCtt                                              <  2:6839957/93‑1 (MQ=255)

                                             tGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTg                                            >  1:1276185/1‑93 (MQ=255)

                                              ggTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATTCGCAGGCGTTCAACGAATCCTTTGc                                           >  2:6996201/1‑93 (MQ=255)

                                                tGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGACCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGccc                                         <  2:3815998/93‑1 (MQ=255)

                                                 gactgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCg                                        <  1:317082/93‑1 (MQ=255)

                                                  accgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGt                                       >  1:5758026/1‑93 (MQ=255)

                                                   ccgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCGACGAAGCCTTTGCCCGTg                                      <  1:5974151/93‑1 (MQ=255)

                                                    cgaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGc                                     >  1:3573195/1‑93 (MQ=255)

                                                     gtccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    <  2:2103447/91‑1 (MQ=255)

                                                     gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    >  1:6869914/1‑93 (MQ=255)

                                                     gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    <  2:6091982/93‑1 (MQ=255)

                                                     gaccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGcc                                    <  2:2193682/93‑1 (MQ=255)

                                                      accTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCt                                   >  2:1430058/1‑93 (MQ=255)

                                                       ccTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTg                                  <  1:7978076/93‑1 (MQ=255)

                                                        ctaaCGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTgg                                 >  2:5772308/4‑93 (MQ=255)

                                                        cTGACGCTGCGTTTTGATCCTGAGTTCGAGCAGAACTCCCGTCGGTTCCTCAACGCCCCCCAGGCGTTCAACGAAGCCTTTGCCGTTGCCTgg                                 >  1:2890164/1‑93 (MQ=255)

                                                        cTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTgg                                 <  1:2542756/93‑1 (MQ=255)

                                                         tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCATCGCAGCCTTTGCCCGTGCCTGGt                                >  2:5986025/1‑93 (MQ=255)

                                                         tGACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGt                                >  2:1676810/1‑93 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               >  2:1805420/1‑93 (MQ=255)

                                                          gACGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGtt                               <  1:8662470/93‑1 (MQ=255)

                                                           aCGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTc                              <  2:3282413/93‑1 (MQ=255)

                                                            cGCTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCa                             >  2:3092889/1‑93 (MQ=255)

                                                              cTGCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCaaa                           >  2:4648720/1‑93 (MQ=255)

                                                                gCGTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACt                         <  1:4951603/93‑1 (MQ=255)

                                                                  gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTTGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGa                       >  2:2237014/1‑93 (MQ=255)

                                                                  gTTTTGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACtta                       >  1:5378474/1‑91 (MQ=255)

                                                                   ttttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAc                      >  1:3533017/1‑93 (MQ=255)

                                                                   ttttGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGAc                      >  2:7851294/1‑93 (MQ=255)

                                                                      tGATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGca                   <  2:3671144/93‑1 (MQ=255)

                                                                       gATCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcac                  >  2:1389880/1‑93 (MQ=255)

                                                                        aTCCTGAGTTCGAGAAGACCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                 <  2:6355060/93‑1 (MQ=255)

                                                                        aTCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGcaca                 <  1:2159889/93‑1 (MQ=255)

                                                                         tCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAg                <  1:3099868/93‑1 (MQ=255)

                                                                         tCCTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAg                <  2:7601725/93‑1 (MQ=255)

                                                                          ccTGAGTTCGAGAAGAACTCTCTTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTTCCCGTGCCTGGGTCAAACTGACGTACAgg               >  1:7819880/1‑93 (MQ=255)

                                                                          ccTGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAgg               >  2:80277/1‑93 (MQ=255)

                                                                            tGAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGGAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGa             >  1:1988849/1‑93 (MQ=255)

                                                                             gAGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGat            <  2:6517862/93‑1 (MQ=255)

                                                                              aggtCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGata           <  1:2724837/90‑1 (MQ=255)

                                                                              aGTTCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGata           <  1:752090/93‑1 (MQ=255)

                                                                                 tggagaAGATCCCCAGTCGTTGCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgg        <  2:4461885/91‑1 (MQ=255)

                                                                                 tCGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTCACGCACAGGGATATgg        >  1:1403506/1‑93 (MQ=255)

                                                                                  gtaaaagaGCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTGGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATggg       <  1:2095422/89‑1 (MQ=255)

                                                                                  cGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATggg       <  2:2796928/93‑1 (MQ=255)

                                                                                  cGAGAAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATggg       <  2:2948399/93‑1 (MQ=255)

                                                                                   gagaAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgggg      >  2:6529704/1‑93 (MQ=255)

                                                                                   gagaAGAACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATgggg      <  2:1195738/93‑1 (MQ=255)

                                                                                     gaagaaCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCACCGAAGCCTTGGCCCGGGCCTGGTCCAAACTGACGAACAGGAATATGGGGcg    >  1:4502125/1‑92 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCg   >  1:2151996/1‑93 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCg   <  2:5652493/93‑1 (MQ=255)

                                                                                      aagaACTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGCGATATGGGGCCg   >  1:5693931/1‑93 (MQ=255)

TCCCGGATCCGTTTGATCCGTCGAAGAAACGTAAACCGACAATGCTGGTGACCGACCTGACGCTGCGTTTTGATCCTGAGTTCGAGAAGATCTCTCGTCGTTTCCTCAACGATCCGCAGGCGTTCAACGAAGCCTTTGCCCGTGCCTGGTTCAAACTGACGCACAGGGATATGGGGCCGA  >  NC_000913/4132490‑4132669

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