BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,923,857	T→C	Q76R (CAG→CGG)	mioC ←	MioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,923,857	0	T	C	98.9%	578.5 / ‑2.9	180	Q76R (CAG→CGG)	mioC	MioC
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (79/99); total (79/101)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

Only 100 of 183 total aligned reads displayed.

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAG  >  NC_000913/3923767‑3923946

tCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGt                                                                                         <  1:3105887/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                      <  1:1024199/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                      <  2:2845817/93‑1 (MQ=255)

   aGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTcc                                                                                      <  1:3926731/93‑1 (MQ=255)

     tttATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTg                                                                                    >  2:4280318/1‑93 (MQ=255)

       tATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCa                                                                                  <  1:613457/93‑1 (MQ=255)

        aTCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCaa                                                                                 <  1:3843269/93‑1 (MQ=255)

        aTCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCaa                                                                                 <  2:3839882/93‑1 (MQ=255)

         tCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAAt                                                                                <  2:900017/93‑1 (MQ=255)

         tCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAAt                                                                                <  2:304173/93‑1 (MQ=255)

          cGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATg                                                                               <  2:443556/93‑1 (MQ=255)

           gATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATCCCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGc                                                                              <  1:2824245/93‑1 (MQ=255)

             tAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCtt                                                                            >  2:1764788/1‑93 (MQ=255)

              aGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTc                                                                           <  2:2501635/93‑1 (MQ=255)

              aGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTc                                                                           >  1:3820748/1‑93 (MQ=255)

              aGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTc                                                                           <  1:3161571/93‑1 (MQ=255)

              aGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTc                                                                           <  2:2209895/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          <  1:4727335/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          <  1:3146693/93‑1 (MQ=255)

               gCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCa                                                                          <  2:802980/93‑1 (MQ=255)

                ccccACAAAAGGTGTCATATTCACGCCTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCAt                                                                         <  2:4456756/93‑1 (MQ=255)

                 cccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATa                                                                        >  2:2781422/1‑93 (MQ=255)

                 cccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATa                                                                        <  1:1038302/93‑1 (MQ=255)

                  ccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGTCTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATag                                                                       >  2:946990/1‑93 (MQ=255)

                  ccACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATag                                                                       >  1:806152/1‑93 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAATGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      >  1:2171423/1‑93 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      >  1:2355043/1‑93 (MQ=255)

                   cacaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATaga                                                                      >  1:4364414/1‑93 (MQ=255)

                    acaAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaa                                                                     >  1:702139/1‑93 (MQ=255)

                     caAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaa                                                                    <  1:3749662/93‑1 (MQ=255)

                     caAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaa                                                                    <  1:1957048/93‑1 (MQ=255)

                      aaaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaag                                                                   >  1:2555126/1‑93 (MQ=255)

                      aaaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaag                                                                   >  1:4557570/1‑93 (MQ=255)

                      aaaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaag                                                                   <  2:4584628/93‑1 (MQ=255)

                       aaaGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaagg                                                                  <  2:555938/93‑1 (MQ=255)

                         aGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTTCTGCAATGCTTCATagaaaggag                                                                >  2:3397609/1‑93 (MQ=255)

                           gtgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATagaaaggagaa                                                              <  1:1566182/93‑1 (MQ=255)

                            tgtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaa                                                             <  2:3487077/93‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                            <  2:1909286/93‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                            <  2:1262563/93‑1 (MQ=255)

                             gtCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATagaaaggagaaag                                                            <  2:511846/93‑1 (MQ=255)

                               cATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGt                                                          >  2:4549451/1‑93 (MQ=255)

                                atatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGtt                                                         <  2:1978995/93‑1 (MQ=255)

                                 tatTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAGAGGTTg                                                        >  1:1967234/1‑93 (MQ=255)

                                    tCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTcc                                                     >  1:4050/1‑93 (MQ=255)

                                     cACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCg                                                    <  1:4238819/93‑1 (MQ=255)

                                      aCGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCgg                                                   >  1:2853830/1‑93 (MQ=255)

                                         aCTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAat                                                >  2:3491071/1‑93 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                               <  2:4657996/93‑1 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                               <  2:3631531/93‑1 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                               <  2:2233/93‑1 (MQ=255)

                                          cTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAata                                               <  2:1523905/93‑1 (MQ=255)

                                           tGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAatat                                              >  2:2685462/1‑93 (MQ=255)

                                             ccAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtct                                            >  1:3006209/1‑93 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                           <  2:1740470/93‑1 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                           >  2:1176799/1‑93 (MQ=255)

                                              cAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATAtctc                                           <  1:4318256/93‑1 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                          <  2:2055395/93‑1 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                          <  2:1618671/93‑1 (MQ=255)

                                               aaTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTcc                                          >  1:1610300/1‑93 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         >  1:2993895/1‑93 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         >  2:3046347/1‑93 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         <  1:4270061/93‑1 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         >  2:4150201/1‑93 (MQ=255)

                                                aTACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCg                                         >  2:4002857/1‑93 (MQ=255)

                                                    cGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAc                                     >  2:2460071/1‑93 (MQ=255)

                                                     gATTGCGCCAAAGCGGACTGCAGAAAGATCTGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAcc                                    <  1:2511614/93‑1 (MQ=255)

                                                     gATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAcc                                    >  2:2187401/1‑93 (MQ=255)

                                                     gATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCAcc                                    >  2:2322026/1‑93 (MQ=255)

                                                      aTTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCg                                   >  1:1718165/1‑93 (MQ=255)

                                                        tGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtg                                 >  2:4133015/1‑93 (MQ=255)

                                                          cgcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtggg                               <  1:4580147/93‑1 (MQ=255)

                                                           gcCAAAGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCgtgggt                              <  2:1653010/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAACGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           <  2:276142/93‑1 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           >  2:1436363/1‑93 (MQ=255)

                                                              aaaGCGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGa                           >  2:1145424/1‑93 (MQ=255)

                                                                  cGGACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGGGGAGCTg                       >  2:2387457/1‑93 (MQ=255)

                                                                   ggACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGa                      >  2:56467/1‑93 (MQ=255)

                                                                   ggACTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGa                      <  2:4064788/93‑1 (MQ=255)

                                                                      cTGCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATaa                   >  2:1669587/1‑93 (MQ=255)

                                                                        gCAGAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAcc                 <  1:2315732/93‑1 (MQ=255)

                                                                           gAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagc              >  2:256605/1‑93 (MQ=255)

                                                                           gAAAGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagc              >  2:4485890/1‑93 (MQ=255)

                                                                             aaGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagcca            >  2:2429466/1‑93 (MQ=255)

                                                                              aGATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAAccagccag           >  1:4728841/1‑93 (MQ=255)

                                                                               gATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGa          <  2:3756157/93‑1 (MQ=255)

                                                                               gATCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGa          >  2:3319361/1‑93 (MQ=255)

                                                                                aTCGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGAt         >  2:3583230/1‑93 (MQ=255)

                                                                                  cGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATcc       <  1:2355394/93‑1 (MQ=255)

                                                                                  cGGGCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATcc       >  2:4513202/1‑93 (MQ=255)

                                                                                   gggCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATccc      <  1:1232875/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    <  1:3604888/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTg    <  1:580998/93‑1 (MQ=255)

                                                                                     gCTTCCGTTCCTGCAATGCTTCATAGAAAGGAG‑AAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa   <  1:681960/93‑1 (MQ=255)

                                                                                      cTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa   <  1:4197219/93‑1 (MQ=255)

                                                                                      cTTCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGa   <  1:1454224/93‑1 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  <  2:3248698/93‑1 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  >  2:3183765/1‑93 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  >  1:3336230/1‑93 (MQ=255)

                                                                                       ttCCGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAg  >  1:3529986/1‑93 (MQ=255)

TCGAGTTTATCGATAGCCCCACAAAAGGTGTCATATTCACGACTGCCAATACCGATTGCGCCAAAGCGGACTGCAGAAAGATCGGGCTTCTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAG  >  NC_000913/3923767‑3923946

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