BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,281,741	T→C	D76G (GAC→GGC)	yejK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,281,741	0	T	C	99.1%	352.0 / ‑5.1	113	D76G (GAC→GGC)	yejK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (63/49); total (63/50)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 122 total aligned reads displayed.

AAAACAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGTCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCTAT  >  NC_000913/2281652‑2281828

aaaaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCct                                                                                      >  2:793494/1‑93 (MQ=255)

 aaaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCctc                                                                                     >  2:3507313/1‑93 (MQ=255)

 aaaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCctc                                                                                     >  2:1174649/1‑93 (MQ=255)

  aaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCctct                                                                                    >  1:2937091/1‑93 (MQ=255)

  aaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCctct                                                                                    <  1:2658477/93‑1 (MQ=255)

  aaCAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGAGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCctct                                                                                    >  1:4320549/1‑93 (MQ=255)

    aaGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTc                                                                                  <  2:3462096/92‑1 (MQ=255)

     aGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGTGGCCGGTTGCCGCGAGGCTAAATGCCAGGAAGCCCTCTTct                                                                                 >  2:945222/1‑92 (MQ=255)

      gCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGc                                                                                <  2:4697379/93‑1 (MQ=255)

      gCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGc                                                                                <  2:4689070/93‑1 (MQ=255)

         cGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCt                                                                             <  1:2924337/93‑1 (MQ=255)

          gAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTg                                                                            >  2:2001293/1‑93 (MQ=255)

          gAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTg                                                                            >  1:297736/1‑93 (MQ=255)

           aaTCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGa                                                                           >  2:1514017/1‑93 (MQ=255)

           aaTCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGa                                                                           >  2:3342468/1‑93 (MQ=255)

           aaTCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGa                                                                           <  1:4488604/93‑1 (MQ=255)

            aTCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGAc                                                                          >  2:3908560/1‑93 (MQ=255)

              ccACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGc                                                                        <  1:693633/93‑1 (MQ=255)

               cACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCt                                                                       >  2:2560760/1‑93 (MQ=255)

               cACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCt                                                                       >  2:917570/1‑93 (MQ=255)

                 ccGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGt                                                                     <  2:3061923/93‑1 (MQ=255)

                  cGTGAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTa                                                                    <  1:1907049/93‑1 (MQ=255)

                   gTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTaa                                                                   <  2:2765451/93‑1 (MQ=255)

                     cAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAgc                                                                 >  2:450438/1‑93 (MQ=255)

                     cAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAgc                                                                 >  2:1210123/1‑93 (MQ=255)

                       gCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAgcgc                                                               >  1:1818716/1‑93 (MQ=255)

                        cAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCa                                                              <  2:676030/93‑1 (MQ=255)

                        cAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCa                                                              <  2:676157/93‑1 (MQ=255)

                         aaaaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAg                                                             >  2:598819/1‑93 (MQ=255)

                         aaaaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGATGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAg                                                             >  2:3005088/1‑93 (MQ=255)

                          aaaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGc                                                            <  1:715260/93‑1 (MQ=255)

                          aaaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGc                                                            <  1:2606292/93‑1 (MQ=255)

                           aaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCg                                                           >  2:3877951/1‑93 (MQ=255)

                           aaGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCg                                                           <  1:3461236/93‑1 (MQ=255)

                            aGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGt                                                          >  1:3860023/1‑93 (MQ=255)

                            aGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGt                                                          >  2:4503126/1‑93 (MQ=255)

                             ggATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGACCTGACGCTGTAAGCGCAGCGTc                                                         <  1:3426946/93‑1 (MQ=255)

                              gATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCt                                                        >  1:2608634/1‑93 (MQ=255)

                                tatTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTgc                                                      >  2:714141/1‑93 (MQ=255)

                                tatTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTgc                                                      >  2:714142/1‑93 (MQ=255)

                                    tCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCa                                                  <  1:2315391/93‑1 (MQ=255)

                                      gCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGt                                                >  2:3993285/1‑93 (MQ=255)

                                        tAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGGAGTAATCTTCGCCCTGACGCTGTGAGCGCAGCGTCTGCGTCAGTTc                                              >  2:2693791/1‑93 (MQ=255)

                                         aaTTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCa                                             >  1:1025591/1‑93 (MQ=255)

                                         aaTTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCa                                             >  1:1025589/1‑93 (MQ=255)

                                          aTTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGTCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCAc                                            <  1:2140896/93‑1 (MQ=255)

                                          aTTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCAc                                            <  1:2809813/93‑1 (MQ=255)

                                              gTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctc                                        >  1:1805334/1‑93 (MQ=255)

                                               tCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       <  1:743060/93‑1 (MQ=255)

                                               tCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       <  1:122134/93‑1 (MQ=255)

                                               tCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       <  2:1453579/93‑1 (MQ=255)

                                               tCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       >  2:3433597/1‑93 (MQ=255)

                                               tCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       <  1:391488/93‑1 (MQ=255)

                                               tCACGCAGGAGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCActctct                                       >  1:909580/1‑93 (MQ=255)

                                                cACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCtt                                      <  1:941841/93‑1 (MQ=255)

                                                 aCGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTc                                     >  2:3334852/1‑93 (MQ=255)

                                                 aCGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTc                                     <  1:175632/93‑1 (MQ=255)

                                                 aCGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTc                                     <  1:95615/93‑1 (MQ=255)

                                                    cAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCt                                  >  1:2287210/1‑93 (MQ=255)

                                                     aGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTa                                 >  1:2651971/1‑93 (MQ=255)

                                                     aGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTa                                 >  2:2721254/1‑93 (MQ=255)

                                                         ggcCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACa                             >  2:239923/1‑93 (MQ=255)

                                                          gcCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACac                            >  1:527662/1‑92 (MQ=255)

                                                          gcCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACaa                            <  1:2063494/93‑1 (MQ=255)

                                                          gcCGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACaa                            <  2:40265/93‑1 (MQ=255)

                                                            cGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGc                          <  1:1115494/93‑1 (MQ=255)

                                                            cGGTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGc                          <  1:1228043/93‑1 (MQ=255)

                                                             tgTTGCCGCGGGGCTAAATGCCTGGAAGCACTCTTCGCCATGACGCTGTAAGCGCAGCGTCTGGGCCAGTTCACGCTCTTCGTTAAACTAGcc                         >  2:916938/2‑93 (MQ=255)

                                                             ggTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGcc                         >  2:26364/1‑93 (MQ=255)

                                                             ggTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGcc                         <  1:4560919/93‑1 (MQ=255)

                                                             ggTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGcc                         >  2:3993497/1‑93 (MQ=255)

                                                             ggTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGcc                         <  2:3965538/93‑1 (MQ=255)

                                                              gTTGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCg                        <  2:4731316/93‑1 (MQ=255)

                                                               ttGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGt                       >  2:3853802/1‑93 (MQ=255)

                                                               ttGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGt                       <  2:4454213/93‑1 (MQ=255)

                                                                tGCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTa                      >  2:2385253/1‑93 (MQ=255)

                                                                 gCCGCGCGGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTaa                     <  1:4562130/93‑1 (MQ=255)

                                                                   cgcgcgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGc                   <  1:376174/93‑1 (MQ=255)

                                                                      gcgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCttt                >  1:147786/1‑93 (MQ=255)

                                                                      gcgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCttt                >  2:3647723/1‑93 (MQ=255)

                                                                       gggCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTa               <  2:3865805/92‑1 (MQ=255)

                                                                       gggCTAAATGCCAGGAAGCCCTCTTCGCCCAGTGACTGGAAGCGCAGCGTGTGCGCCACTCCACCCCCTCCCCAAACCAAGCCGTAAGCTTTa               <  2:1350236/92‑1 (MQ=255)

                                                                       cgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTa               <  2:4567579/93‑1 (MQ=255)

                                                                       cgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTa               >  1:2619760/1‑93 (MQ=255)

                                                                       cgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTa               <  1:4619638/93‑1 (MQ=255)

                                                                       cgGCTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTa               <  2:3813011/93‑1 (MQ=255)

                                                                        ggCAAAATCCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCGGCGCCAGTTCACTCTCTTCGCAAAACAAGCCGTAAGCTTTAt              <  2:2918058/93‑1 (MQ=255)

                                                                          cTAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTAttt            >  1:3302459/1‑93 (MQ=255)

                                                                           tAAATGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTAtttt           >  1:4496579/1‑93 (MQ=255)

                                                                            aaaTGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTAttttt          >  1:88451/1‑93 (MQ=255)

                                                                            aaaTGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTAttttt          >  1:88452/1‑93 (MQ=255)

                                                                             aaTGCCAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTg         >  2:1856353/1‑93 (MQ=255)

                                                                              aTGCCAGGAAGCCCTCTTCGCCCTGACGCTGAAAGCGCAGCGTCTCCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTgg        <  2:3383699/93‑1 (MQ=255)

                                                                               tGCCAGGAAGCCCTCTCCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGgc       <  1:3363149/93‑1 (MQ=255)

                                                                               tGCCAGGAACCCCTCTTCGCCCGACCGCTGTAAGCGCAGCGTCTGCGCCATTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGgc       <  2:329362/93‑1 (MQ=255)

                                                                               tGCCAGAAACCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCCCTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGgc       <  2:754296/93‑1 (MQ=255)

                                                                                  cAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCt    >  2:4286391/1‑93 (MQ=255)

                                                                                  cAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCt    >  2:2090678/1‑93 (MQ=255)

                                                                                  cAGGAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGCTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCt    >  1:3116502/1‑93 (MQ=255)

                                                                                    ggAAGCCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCtat  >  2:2017722/1‑93 (MQ=255)

AAAACAGCACGAATCCACCGTCAGCAAAAGGATATTTCGCTAATTCGTCACGCAGGCGGCCGGTTGCCGCGCGGCTAAATGCCAGGAAGTCCTCTTCGCCCTGACGCTGTAAGCGCAGCGTCTGCGCCAGTTCACTCTCTTCGCTAAACAAGCCGTAAGCTTTATTTTTGGCGCTAT  >  NC_000913/2281652‑2281828

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: `‑`