BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,616,692	C→T	G81S (GGC→AGC)	ydeB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,616,692	0	C	T	99.5%	752.7 / ‑5.6	203	G81S (GGC→AGC)	ydeB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (117/85); total (118/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 207 total aligned reads displayed.

TCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616604‑1616781

tctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCt                                                                                       >  2:3909153/1‑93 (MQ=255)

 ctTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTg                                                                                      >  2:1798203/1‑93 (MQ=255)

  ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGctgcctgc                                                                                     <  2:2511956/93‑1 (MQ=255)

  ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGctgcctgc                                                                                     <  2:5845501/93‑1 (MQ=255)

  ttCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGc                                                                                     >  2:2812674/1‑93 (MQ=255)

   tCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCa                                                                                    >  1:2127698/1‑93 (MQ=255)

    ccGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCaa                                                                                   <  1:3784047/93‑1 (MQ=255)

     cGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAt                                                                                  >  1:2371623/1‑93 (MQ=255)

     cGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTAGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAt                                                                                  >  2:5136572/1‑93 (MQ=255)

      gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  2:3736768/1‑93 (MQ=255)

      gACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAAtt                                                                                 >  1:2968213/1‑93 (MQ=255)

         ttGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGc                                                                              >  2:4675160/1‑93 (MQ=255)

          tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             >  2:2761047/1‑93 (MQ=255)

          tGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCa                                                                             <  1:1208270/93‑1 (MQ=255)

            cTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGa                                                                           >  1:3366785/1‑93 (MQ=255)

             ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  1:2239569/93‑1 (MQ=255)

             ttttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGAc                                                                          <  2:1505709/93‑1 (MQ=255)

               ttCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCg                                                                        >  1:2435238/1‑93 (MQ=255)

                tCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCgg                                                                       <  1:4517620/93‑1 (MQ=255)

                 cGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCggg                                                                      <  1:3236618/93‑1 (MQ=255)

                  gCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGa                                                                     <  1:2851098/93‑1 (MQ=255)

                  gCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGa                                                                     >  1:901905/1‑93 (MQ=255)

                    cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCCGACCGGGAAt                                                                   >  1:766516/1‑93 (MQ=255)

                    cTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAAt                                                                   >  1:5486025/1‑93 (MQ=255)

                     tCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATc                                                                  >  2:279653/1‑93 (MQ=255)

                      ccGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCg                                                                 >  2:4124119/1‑93 (MQ=255)

                      ccGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCg                                                                 >  1:3068741/1‑93 (MQ=255)

                      ccGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCg                                                                 >  1:5299078/1‑93 (MQ=255)

                       cGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGa                                                                >  2:4161602/1‑93 (MQ=255)

                       cGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGa                                                                >  1:5773210/1‑93 (MQ=255)

                         gTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGaaa                                                              <  2:2433931/93‑1 (MQ=255)

                          tGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAAt                                                             >  2:4925399/1‑93 (MQ=255)

                          tGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAAt                                                             >  1:5866468/1‑93 (MQ=255)

                           gAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATg                                                            <  2:5215532/93‑1 (MQ=255)

                              tcaatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCa                                                         >  2:3199442/1‑93 (MQ=255)

                               caatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaa                                                        <  2:1691265/93‑1 (MQ=255)

                                aatcGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCaaa                                                       >  1:1662639/1‑93 (MQ=255)

                                  tcGCTTTCTGTTGCGGAAAGAGCATCCGAAACCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATg                                                     <  2:2214823/93‑1 (MQ=255)

                                   cGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt                                                    <  1:1043140/93‑1 (MQ=255)

                                   cGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGt                                                    <  2:3151486/93‑1 (MQ=255)

                                    gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  1:166915/1‑93 (MQ=255)

                                    gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  1:3895395/1‑93 (MQ=255)

                                    gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  1:4151303/1‑93 (MQ=255)

                                    gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   >  2:1025988/1‑93 (MQ=255)

                                    gCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTa                                                   <  1:654309/93‑1 (MQ=255)

                                      tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 <  2:3721467/93‑1 (MQ=255)

                                      tttCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATc                                                 >  2:5541861/1‑93 (MQ=255)

                                        tCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa                                               >  1:50582/1‑93 (MQ=255)

                                        tCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa                                               <  2:439581/93‑1 (MQ=255)

                                        tCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCa                                               <  2:4540936/93‑1 (MQ=255)

                                         cTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCAt                                              <  1:2825295/93‑1 (MQ=255)

                                         cTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCAt                                              >  2:5382263/1‑93 (MQ=255)

                                          tGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATc                                             <  1:4643362/93‑1 (MQ=255)

                                           gTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCa                                            >  1:4678648/1‑93 (MQ=255)

                                           gTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCa                                            >  1:4495190/1‑93 (MQ=255)

                                               cGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAg                                        <  1:3507532/93‑1 (MQ=255)

                                                   aaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGc                                    >  2:3779547/1‑93 (MQ=255)

                                                   aaGAGCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGc                                    >  1:3926968/1‑93 (MQ=255)

                                                    agagCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGcc                                   >  1:569496/1‑93 (MQ=255)

                                                      agCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCg                                 <  1:1445383/93‑1 (MQ=255)

                                                       gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                >  2:2151691/1‑93 (MQ=255)

                                                       gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                <  1:5999549/93‑1 (MQ=255)

                                                       gCATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCgg                                >  2:2882208/1‑93 (MQ=255)

                                                        cATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGc                               <  1:3169829/93‑1 (MQ=255)

                                                        cATCCGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGc                               <  2:34259/93‑1 (MQ=255)

                                                           ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            >  1:4858364/1‑93 (MQ=255)

                                                           ccGAAAACCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTa                            <  2:4761663/93‑1 (MQ=255)

                                                              aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         <  2:4015649/93‑1 (MQ=255)

                                                              aaaaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATa                         >  2:5442339/1‑93 (MQ=255)

                                                                aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       <  2:4034489/93‑1 (MQ=255)

                                                                aaCCGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACg                       >  1:2625170/1‑93 (MQ=255)

                                                                  ccGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcc                     >  2:2715908/1‑93 (MQ=255)

                                                                  ccGATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGcc                     <  1:2032227/93‑1 (MQ=255)

                                                                    gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   >  2:3005235/1‑93 (MQ=255)

                                                                    gATAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccac                   <  1:3970310/93‑1 (MQ=255)

                                                                     aTAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacc                  >  2:6068241/1‑93 (MQ=255)

                                                                      tAAACGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGccacca                 <  1:1469089/93‑1 (MQ=255)

                                                                       aaaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcat                >  1:2191975/1‑93 (MQ=255)

                                                                       aaaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcat                <  1:6194226/93‑1 (MQ=255)

                                                                        aaCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatc               >  2:5127335/1‑93 (MQ=255)

                                                                         aCGCGACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatca              <  1:6290755/93‑1 (MQ=255)

                                                                          cgcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatcat             <  1:165458/93‑1 (MQ=255)

                                                                          cgcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatcat             <  1:426486/93‑1 (MQ=255)

                                                                          cgcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatcat             <  1:4762281/93‑1 (MQ=255)

                                                                           gcgACGATTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACcatcata            >  2:1875387/1‑92 (MQ=255)

                                                                            cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           >  2:4884416/1‑93 (MQ=255)

                                                                            cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           >  1:1101130/1‑93 (MQ=255)

                                                                            cgacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatga           >  2:3587803/1‑93 (MQ=255)

                                                                             gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:3395398/93‑1 (MQ=255)

                                                                             gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          <  2:1203608/93‑1 (MQ=255)

                                                                             gacgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCatgat          >  1:329843/1‑93 (MQ=255)

                                                                               cgaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATgg        <  1:1754112/93‑1 (MQ=255)

                                                                                gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:3047812/1‑93 (MQ=255)

                                                                                gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:630371/1‑93 (MQ=255)

                                                                                gaTTAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGc       >  2:806424/1‑93 (MQ=255)

                                                                                   tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    >  2:3307894/1‑93 (MQ=255)

                                                                                   tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  2:2053306/93‑1 (MQ=255)

                                                                                   tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:983173/93‑1 (MQ=255)

                                                                                   tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:6292418/93‑1 (MQ=255)

                                                                                   tAAGCTGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCaaa    <  1:5944550/93‑1 (MQ=255)

TCTTCCGACTTGCTTTTCGCCTCCGGTGAATCAATCGCTTTCTGTTGCGGAAAGAGCATCCGAAAACCGATAAACGCGACGATTAAGCCGCCTGCAATTCGCAGACCGGGAATCGAAATGCCAAATGTATCCATCACCAGTTGCCCGGCGTAATACGCCACCATCATGATGGCAAATA  >  NC_000913/1616604‑1616781

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