BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,477,397	T→C	N283S (AAC→AGC)	kefB ←	KefB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,477,397	0	T	C	83.6%	612.7 / 109.0	274	N283S (AAC→AGC)	kefB	KefB
Reads supporting (aligned to +/- strand): ref base T (18/27); new base C (101/128); total (119/155)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCAC  >  NC_000913/3477305‑3477487

gcgcCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGt                                                                                            <  1:5476847/93‑1 (MQ=255)

  gcCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTg                                                                                          <  2:6346334/93‑1 (MQ=255)

  gcCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTg                                                                                          >  1:5316964/1‑93 (MQ=255)

   ccAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGa                                                                                         >  2:844087/1‑93 (MQ=255)

   ccAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGa                                                                                         <  1:1509711/93‑1 (MQ=255)

    cagcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAg                                                                                        <  2:5533891/93‑1 (MQ=255)

    cagcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAg                                                                                        <  1:6411931/93‑1 (MQ=255)

      gcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTg                                                                                      >  2:5025074/1‑93 (MQ=255)

      gcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTg                                                                                      <  2:3957079/93‑1 (MQ=255)

       cagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGa                                                                                     <  2:2327888/93‑1 (MQ=255)

        agATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGAc                                                                                    >  2:34276/1‑93 (MQ=255)

        agATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGAc                                                                                    >  1:2638750/1‑93 (MQ=255)

        agATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGAc                                                                                    <  2:3343794/93‑1 (MQ=255)

        agATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGAc                                                                                    >  1:1463202/1‑93 (MQ=255)

         gATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACa                                                                                   >  2:4473478/1‑93 (MQ=255)

         gATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACa                                                                                   <  2:3772583/93‑1 (MQ=255)

         gATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACa                                                                                   >  1:2519976/1‑93 (MQ=255)

           tACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATg                                                                                 >  1:86301/1‑93 (MQ=255)

            aCAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGc                                                                                <  1:1158692/93‑1 (MQ=255)

            aCAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGc                                                                                <  1:2774328/93‑1 (MQ=255)

              aGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCg                                                                              >  1:2279403/1‑93 (MQ=255)

              aGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCg                                                                              >  2:5342249/1‑93 (MQ=255)

               gCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGa                                                                             >  2:1534385/1‑93 (MQ=255)

                cACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGAc                                                                            >  1:5539528/1‑93 (MQ=255)

                 aCGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACa                                                                           >  2:5982339/1‑93 (MQ=255)

                 aCGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACa                                                                           >  2:6877254/1‑93 (MQ=255)

                  cGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACag                                                                          >  2:1322736/1‑93 (MQ=255)

                  cGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACag                                                                          >  1:3050545/1‑93 (MQ=255)

                  cGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACag                                                                          >  1:4160740/1‑93 (MQ=255)

                   gagaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACaga                                                                         <  2:6383758/93‑1 (MQ=255)

                   gagaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACaga                                                                         <  2:2115384/93‑1 (MQ=255)

                   gagaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACaga                                                                         <  1:2251073/93‑1 (MQ=255)

                     gaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAgaga                                                                       <  2:666555/93‑1 (MQ=255)

                       aTTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATa                                                                     >  1:6655696/1‑93 (MQ=255)

                        ttttCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATaa                                                                    >  2:1889744/1‑93 (MQ=255)

                         tttCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATaaa                                                                   <  2:3823790/93‑1 (MQ=255)

                         tttCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATaaa                                                                   <  1:163186/93‑1 (MQ=255)

                           tCATCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGa                                                                 >  1:6974308/1‑93 (MQ=255)

                           tCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGa                                                                 <  1:3930603/93‑1 (MQ=255)

                            cACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGaa                                                                <  2:1687578/93‑1 (MQ=255)

                             aCCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAAc                                                               <  1:1313982/93‑1 (MQ=255)

                              ccGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACa                                                              >  2:6939276/1‑93 (MQ=255)

                                 ccaccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAAc                                                           <  2:2716841/93‑1 (MQ=255)

                                 ccaccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAAc                                                           >  1:4088718/1‑93 (MQ=255)

                                   accaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCg                                                         <  1:5119535/93‑1 (MQ=255)

                                   accaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCg                                                         <  1:809978/93‑1 (MQ=255)

                                    ccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGa                                                        <  2:1727795/93‑1 (MQ=255)

                                    ccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGa                                                        <  2:5940313/93‑1 (MQ=255)

                                    ccaGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGa                                                        >  2:1107101/1‑93 (MQ=255)

                                     caGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGCCATGCCGACAGAGAAAAAGAACAAACCGag                                                       <  1:7399687/93‑1 (MQ=255)

                                      aGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagc                                                      >  1:7034388/1‑93 (MQ=255)

                                       gCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagca                                                     <  1:3383330/93‑1 (MQ=255)

                                        cacaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagcag                                                    >  1:2097038/1‑93 (MQ=255)

                                          caACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAgcagca                                                  >  2:4042158/1‑93 (MQ=255)

                                            aCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAg                                                <  1:6391931/93‑1 (MQ=255)

                                               aCGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCt                                             >  1:2017376/1‑93 (MQ=255)

                                               aCGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCAAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCt                                             >  2:4943620/1‑93 (MQ=255)

                                                cGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCtt                                            <  1:6475124/93‑1 (MQ=255)

                                                cGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCtt                                            <  1:3202722/93‑1 (MQ=255)

                                                cGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACACGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCtt                                            <  1:3123869/93‑1 (MQ=255)

                                                 gCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCttt                                           >  2:7209822/1‑93 (MQ=255)

                                                  cTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTg                                          <  2:3520784/93‑1 (MQ=255)

                                                   ttATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGa                                         >  2:4197363/1‑93 (MQ=255)

                                                    tATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGaa                                        >  1:381396/1‑93 (MQ=255)

                                                    tATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGaa                                        >  1:6006436/1‑93 (MQ=255)

                                                    tATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGaa                                        <  2:6494871/93‑1 (MQ=255)

                                                      tCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAgg                                      <  2:7086201/93‑1 (MQ=255)

                                                       cACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAggg                                     >  2:4247615/1‑93 (MQ=255)

                                                       cACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAggg                                     <  1:4354127/93‑1 (MQ=255)

                                                        aCTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGa                                    <  2:2366550/93‑1 (MQ=255)

                                                        aCTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGa                                    >  1:5610722/1‑93 (MQ=255)

                                                         cTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGAt                                   <  1:4610181/93‑1 (MQ=255)

                                                           aCCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCg                                 <  2:2027989/93‑1 (MQ=255)

                                                           aCCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCg                                 <  1:6846630/93‑1 (MQ=255)

                                                             ccACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGAt                               <  2:5117780/93‑1 (MQ=255)

                                                              cacaACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATa                              <  1:2456991/93‑1 (MQ=255)

                                                              cacaACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATa                              <  2:4504359/93‑1 (MQ=255)

                                                                caacaGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGc                            >  1:7329494/1‑93 (MQ=255)

                                                                caacaGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGc                            <  2:6149906/93‑1 (MQ=255)

                                                                 aacaGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGcc                           <  2:4729229/93‑1 (MQ=255)

                                                                    aGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGtt                        <  2:512010/93‑1 (MQ=255)

                                                                    aGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGtt                        <  2:4841214/93‑1 (MQ=255)

                                                                     gATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGttt                       >  1:888626/1‑93 (MQ=255)

                                                                     gATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGttt                       >  2:1365274/1‑93 (MQ=255)

                                                                      aTGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTc                      <  2:3083140/93‑1 (MQ=255)

                                                                       tGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTcc                     >  2:2142641/1‑93 (MQ=255)

                                                                        gggTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCa                    >  1:2907470/1‑93 (MQ=255)

                                                                         ggTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAg                   <  1:2683518/93‑1 (MQ=255)

                                                                             taAAGGACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCa               >  1:4703695/1‑93 (MQ=255)

                                                                             taAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCa               >  1:4550019/1‑93 (MQ=255)

                                                                              aaaGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTG‑AGGGATCGATAGCCGTTTCCAGTTCATg             <  1:3462480/93‑1 (MQ=255)

                                                                              aaaGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCAt              <  1:1012321/93‑1 (MQ=255)

                                                                               aaGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATg             >  1:654822/1‑93 (MQ=255)

                                                                               aaGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATg             >  1:3550114/1‑93 (MQ=255)

                                                                                aGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATgg            >  2:2916091/1‑93 (MQ=255)

                                                                                 gCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGc           <  2:3294372/93‑1 (MQ=255)

                                                                                        gAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTc    >  2:3717619/1‑93 (MQ=255)

                                                                                         aGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCa   <  2:5734132/93‑1 (MQ=255)

                                                                                         aGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCa   <  2:2510188/93‑1 (MQ=255)

                                                                                          ggTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCAc  <  1:3866186/93‑1 (MQ=255)

GCGCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCAC  >  NC_000913/3477305‑3477487

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