BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	907,798	A→G	Y240Y (TAT→TAC)	ybjU ←	putative arylsulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,798	0	A	G	99.2%	418.4 / NA	133	Y240Y (TAT→TAC)	ybjU	putative arylsulfatase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (74/58); minor base T (0/1); total (74/59)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

CACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCAC  >  NC_000913/907708‑907886

cacaTCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCgta                                                                                        <  2:1955741/93‑1 (MQ=255)

 ccaTCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCgtat                                                                                       <  1:1376275/92‑1 (MQ=255)

  caTCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCgtata                                                                                      <  2:1486573/93‑1 (MQ=255)

    tCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATc                                                                                    >  2:884621/1‑93 (MQ=255)

    tCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATc                                                                                    >  2:3967230/1‑93 (MQ=255)

    tCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATc                                                                                    >  2:683446/1‑93 (MQ=255)

    tCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATc                                                                                    <  2:559647/93‑1 (MQ=255)

     ccGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATcc                                                                                   <  1:288318/93‑1 (MQ=255)

      cgcgCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATccc                                                                                  >  2:3216829/1‑93 (MQ=255)

       gcgcCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcg                                                                                 >  1:3536321/1‑93 (MQ=255)

        cgcCTGCTTCACGCAGCTGCTCCTCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcgg                                                                                <  2:731016/93‑1 (MQ=255)

          ccTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcggcg                                                                              >  2:2015756/1‑93 (MQ=255)

           cTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcggcgg                                                                             <  1:2346852/93‑1 (MQ=255)

           cTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcggcgg                                                                             >  2:2722630/1‑93 (MQ=255)

           cTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCcggcgg                                                                             >  2:4102922/1‑93 (MQ=255)

            tGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCggcggc                                                                            <  2:3175876/93‑1 (MQ=255)

            tGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCggcggc                                                                            >  2:76642/1‑93 (MQ=255)

             gCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCa                                                                           >  2:2839699/1‑93 (MQ=255)

             gCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCa                                                                           >  2:3931672/1‑93 (MQ=255)

                 cACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCa                                                                       <  2:2479495/93‑1 (MQ=255)

                 cACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCa                                                                       <  2:108775/93‑1 (MQ=255)

                    gCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGaa                                                                    <  1:636242/93‑1 (MQ=255)

                      aGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATg                                                                  <  2:3433103/93‑1 (MQ=255)

                          gctCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCgg                                                              >  2:1849228/1‑93 (MQ=255)

                          gctCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCgg                                                              >  2:3663238/1‑93 (MQ=255)

                          gctCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCgg                                                              >  2:1115184/1‑93 (MQ=255)

                           ctCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGa                                                             <  1:2355101/93‑1 (MQ=255)

                            cccGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGAc                                                            <  1:1703064/92‑1 (MQ=255)

                             ccGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACt                                                           <  1:3870344/93‑1 (MQ=255)

                             ccGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACt                                                           >  1:3506189/1‑93 (MQ=255)

                               gCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTgg                                                         >  2:3832526/1‑93 (MQ=255)

                               gCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTgg                                                         >  1:2612761/1‑93 (MQ=255)

                               gCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTgg                                                         <  1:2543083/93‑1 (MQ=255)

                               gCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTgg                                                         >  2:230074/1‑93 (MQ=255)

                                ccATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGgc                                                        <  2:2612406/93‑1 (MQ=255)

                                 cATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGgcg                                                       >  1:4021021/1‑93 (MQ=255)

                                  aTCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGgcgc                                                      >  2:2641456/1‑93 (MQ=255)

                                    ccAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTTGCGCAt                                                    <  2:2986232/93‑1 (MQ=255)

                                    ccAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCAt                                                    <  2:2041298/93‑1 (MQ=255)

                                     cAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATc                                                   >  2:1444872/1‑93 (MQ=255)

                                      aGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATcc                                                  >  1:2737608/1‑93 (MQ=255)

                                        gCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCg                                                >  2:3425078/1‑93 (MQ=255)

                                         cAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGc                                               <  1:4537662/93‑1 (MQ=255)

                                         cAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGc                                               <  1:2505383/93‑1 (MQ=255)

                                         cAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGc                                               >  2:2057003/1‑93 (MQ=255)

                                          aGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGcc                                              >  1:4189728/1‑93 (MQ=255)

                                          aGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGcc                                              >  2:3342971/1‑93 (MQ=255)

                                          aGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGcc                                              >  1:2986799/1‑93 (MQ=255)

                                           gCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCa                                             >  2:2765260/1‑93 (MQ=255)

                                           gCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCa                                             >  2:2169693/1‑93 (MQ=255)

                                           gCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCa                                             <  1:4361948/93‑1 (MQ=255)

                                             gTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCAcc                                           >  1:1394785/1‑93 (MQ=255)

                                             gTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCAcc                                           <  1:3833586/93‑1 (MQ=255)

                                                   gTGGTCTTCCTGCAAGCGCGCACCGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCAt                                     <  1:3472570/93‑1 (MQ=255)

                                                      gTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCAtttt                                  >  2:3135450/1‑93 (MQ=255)

                                                      gTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGCATATCCCGGCGGCAGCCAGCATGCCGGTCTGGCGCATCCCGCCACCTGTTAtttt                                  >  2:968253/1‑93 (MQ=255)

                                                       tCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTc                                 >  1:1173526/1‑93 (MQ=255)

                                                       tCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTc                                 <  1:2902908/93‑1 (MQ=255)

                                                       tCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTc                                 <  1:3883762/93‑1 (MQ=255)

                                                       tCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTc                                 >  2:1086688/1‑93 (MQ=255)

                                                        cTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTcc                                >  1:4059828/1‑93 (MQ=255)

                                                            cTGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCa                            >  1:313135/1‑93 (MQ=255)

                                                             tGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAg                           <  2:3659544/93‑1 (MQ=255)

                                                             tGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAg                           >  2:417210/1‑93 (MQ=255)

                                                             tGCAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAg                           >  1:4392464/1‑93 (MQ=255)

                                                               cAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGTCAGCg                         >  1:3146715/1‑93 (MQ=255)

                                                               cAAGCGCGCAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCg                         <  1:1096609/93‑1 (MQ=255)

                                                                  gcgcgcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAAt                      <  2:4514033/93‑1 (MQ=255)

                                                                  gcgcgcAACCTTATTTTGCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGCCCTGCCGCATCCCGCCCCCTGTCATTTTCCGCCATCGAAt                      <  2:3518650/93‑1 (MQ=255)

                                                                   cgcgcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATg                     <  1:221619/93‑1 (MQ=255)

                                                                    gcgcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATgg                    >  2:3960237/1‑93 (MQ=255)

                                                                    gcgcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATgg                    >  1:1647059/1‑93 (MQ=255)

                                                                     cgcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGc                   <  2:776697/93‑1 (MQ=255)

                                                                      gcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCa                  <  1:4016287/93‑1 (MQ=255)

                                                                      gcAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCa                  >  1:3680688/1‑93 (MQ=255)

                                                                       cAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCAc                 <  2:4402375/93‑1 (MQ=255)

                                                                       cAACGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCAc                 <  2:485367/93‑1 (MQ=255)

                                                                        aaCGTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACg                <  1:3405405/93‑1 (MQ=255)

                                                                          cGTTATTTTGCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGtt              <  1:931020/93‑1 (MQ=255)

                                                                           gTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGttt             >  1:2188113/1‑93 (MQ=255)

                                                                           gTTATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGttt             >  2:3853357/1‑93 (MQ=255)

                                                                           gTTATTTTTCAGGGAGTATATCCCGGCCGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCAcactt             >  2:2000933/1‑89 (MQ=255)

                                                                            ttATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTa            >  2:3323072/1‑93 (MQ=255)

                                                                            ttATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTa            <  1:2727062/93‑1 (MQ=255)

                                                                            ttATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTa            <  2:4491573/93‑1 (MQ=255)

                                                                            ttATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTa            >  1:786227/1‑93 (MQ=255)

                                                                            ttATTTTTCAGGGCGTATATCCCGGCCGCAGCCAGAATGCCGGACTGGCGCATCCCGCCCCCTGTCATTTTCCGCCAGCGAATGGCACGTTTa            <  1:3600718/93‑1 (MQ=255)

                                                                             tATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTaa           <  2:2921256/93‑1 (MQ=255)

                                                                             tATTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTaa           >  2:2319522/1‑93 (MQ=255)

                                                                              aTTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAAt          >  1:3230360/1‑93 (MQ=255)

                                                                              aTTTTTCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAAt          <  2:1200245/93‑1 (MQ=255)

                                                                               tttttCAGGGCTTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATg         <  1:704558/93‑1 (MQ=255)

                                                                                ttttCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCCGCGAATGGCACGTTTAATGt        >  2:2050588/1‑93 (MQ=255)

                                                                                ttttCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGt        <  2:641598/93‑1 (MQ=255)

                                                                                  ttCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTaa      <  2:2405745/93‑1 (MQ=255)

                                                                                   tCAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAAt     <  1:1011888/93‑1 (MQ=255)

                                                                                    cAGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATc    <  1:1969741/93‑1 (MQ=255)

                                                                                     aGGGCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCa   >  2:1983433/1‑93 (MQ=255)

                                                                                      gggCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCCCCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCAc  <  2:2230183/93‑1 (MQ=255)

                                                                                      gggCGTATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCAc  <  1:4157760/93‑1 (MQ=255)

CACATCCGCGCCTGCTTCACGCAGCTGCTCCGCCATCCAGGCAGCGTTGTCGTGGTCTTCCTGCAAGCGCGCAACGTTATTTTTCAGGGCATATATCCCGGCGGCAGCCAGAATGCCGGACTGGCGCATCCCGCCACCTGTCATTTTCCGCCAGCGAATGGCACGTTTAATGTAATCAC  >  NC_000913/907708‑907886

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