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	99.7%	1150.2 / ‑4.9	332	N283S (AAC→AGC)	kefB	KefB
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (146/185); total (146/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.58e-01

Only 100 of 340 total aligned reads displayed.

GCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCA  >  NC_000913/3477307‑3477486

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

gcCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTg                                                                                         <  2:3147244/93‑1 (MQ=255)

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

  cagcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAg                                                                                       >  1:1273392/1‑93 (MQ=255)

   agcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGt                                                                                      >  2:8098326/1‑93 (MQ=255)

    gcagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAggg                                                                                     >  1:6430415/1‑91 (MQ=255)

     cagATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGa                                                                                    >  2:6679582/1‑93 (MQ=255)

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

       aaTACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACa                                                                                  <  2:6352018/92‑1 (MQ=255)

        aTACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACAt                                                                                 <  2:751094/93‑1 (MQ=255)

        aTACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACAt                                                                                 >  1:3020814/1‑93 (MQ=255)

          aCAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGc                                                                               >  2:6993902/1‑93 (MQ=255)

           cAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGcc                                                                              <  1:946912/93‑1 (MQ=255)

           cAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGcc                                                                              >  2:2815824/1‑93 (MQ=255)

           cAGCACGAGAATTTACACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGcc                                                                              >  1:8262217/1‑93 (MQ=255)

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

             gCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGa                                                                            <  2:4510089/93‑1 (MQ=255)

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

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

                 gagaATTTTCACCGCCCCCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACaga                                                                        <  1:2579844/93‑1 (MQ=255)

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

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

                 gagaATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACaga                                                                        >  1:7679589/1‑93 (MQ=255)

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

                   gaATTTTCACCGCCACCAGCACAACCACGCGTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAgaga                                                                      <  1:1146726/93‑1 (MQ=255)

                    aaTTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGAt                                                                     <  1:4033742/93‑1 (MQ=255)

                     aTTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATa                                                                    <  2:6805548/93‑1 (MQ=255)

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

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

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

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

                           cccGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGTTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAAc                                                              <  2:6143161/92‑1 (MQ=255)

                            ccGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACa                                                             <  1:6234969/93‑1 (MQ=255)

                             cGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACaa                                                            >  2:2300827/1‑93 (MQ=255)

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

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

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

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

                                    aGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagc                                                     >  2:8739414/1‑93 (MQ=255)

                                     gCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCTagca                                                    <  2:5600334/93‑1 (MQ=255)

                                     gCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagca                                                    >  2:4192142/1‑93 (MQ=255)

                                       acaACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGagcagc                                                  >  2:5718620/1‑93 (MQ=255)

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

                                         aaCCACGCTTATCTCTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCaa                                                <  1:5736733/93‑1 (MQ=255)

                                         aaCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCaa                                                <  1:723126/93‑1 (MQ=255)

                                         aaCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCaa                                                <  1:5254348/93‑1 (MQ=255)

                                          aCCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAg                                               >  1:7066251/1‑93 (MQ=255)

                                           ccACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGc                                              >  2:3133112/1‑93 (MQ=255)

                                           ccACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGc                                              >  2:2729012/1‑93 (MQ=255)

                                            cACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGcc                                             <  1:6591356/93‑1 (MQ=255)

                                              cGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCtt                                           >  2:7045916/1‑93 (MQ=255)

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

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

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

                                               gCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCttt                                          <  2:8242091/93‑1 (MQ=255)

                                               gCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCttt                                          <  1:1398891/93‑1 (MQ=255)

                                               gCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCttt                                          >  2:6642708/1‑93 (MQ=255)

                                                cTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTg                                         <  2:443690/93‑1 (MQ=255)

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

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

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

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

                                                      aCTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGa                                   <  2:3458082/93‑1 (MQ=255)

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

                                                       cTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGAt                                  >  1:1874867/1‑93 (MQ=255)

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

                                                       cTACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGAt                                  >  2:3245457/1‑93 (MQ=255)

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

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

                                                        tgcccACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATc                                 <  1:1717301/91‑1 (MQ=255)

                                                        tACCCACAACAGATGGGTATAAAGCGCCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATc                                 <  1:8608111/93‑1 (MQ=255)

                                                        tACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATc                                 <  1:2565712/93‑1 (MQ=255)

                                                        tACCCACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATc                                 <  1:2752797/93‑1 (MQ=255)

                                                          cccACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGa                               >  1:5879083/1‑93 (MQ=255)

                                                           ccACAACAGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGAt                              >  1:5389811/1‑93 (MQ=255)

                                                              caacaGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGTAAGCCTTTGAAGGGATCGATAGc                           >  2:4753921/1‑93 (MQ=255)

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

                                                               aacaGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGcc                          >  2:1190760/1‑93 (MQ=255)

                                                                 caGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGt                        <  1:2221997/93‑1 (MQ=255)

                                                                 caGATGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGCGATCGATAGCCGt                        <  2:1608122/93‑1 (MQ=255)

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

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

                                                                  aGATGAGTATAAAGCACCACGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGtt                       <  1:3299230/93‑1 (MQ=255)

                                                                     tGGGTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTcc                    <  1:4415347/93‑1 (MQ=255)

                                                                      gggTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCa                   >  1:3841693/1‑93 (MQ=255)

                                                                       ggTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAg                  <  2:39276/93‑1 (MQ=255)

                                                                        gTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt                 >  1:2041486/1‑93 (MQ=255)

                                                                        gTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt                 <  1:1056079/93‑1 (MQ=255)

                                                                        gTATAAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGt                 >  1:461852/1‑93 (MQ=255)

                                                                           taAAGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCa              <  2:5277148/93‑1 (MQ=255)

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

                                                                             aaGCACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTTATg            >  2:4889819/1‑93 (MQ=255)

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

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

                                                                                cACCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCg         <  2:7658745/93‑1 (MQ=255)

                                                                                 aCCCCGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGa        >  1:813990/1‑93 (MQ=255)

                                                                                    ccGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTTCAGTTCACGGCGatat     >  2:1825264/1‑93 (MQ=255)

                                                                                    ccGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGatat     <  1:7557563/93‑1 (MQ=255)

                                                                                     cGAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATAtt    <  1:677243/93‑1 (MQ=255)

                                                                                      gAGGCTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTc   <  2:1047118/93‑1 (MQ=255)

GCCAGCAGATACAGCACGAGAATTTTCACCGCCACCAGCACAACCACGCTTATCACTACCCACAACAGATGGGTATAAAGCACCCCGAGGTTGAGTGACATGCCGACAGAGATAAAGAACAAACCGAGCAGCAAGCCTTTGAAGGGATCGATAGCCGTTTCCAGTTCATGGCGATATTCA  >  NC_000913/3477307‑3477486

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