BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,170	T→C	I502I (ATT→ATC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,170	0	T	C	100.0%	1435.5 / NA	410	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (231/179); total (231/179)

Only 100 of 414 total aligned reads displayed.

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

 aaTATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaa                                                                                        >  2:7867033/1‑93 (MQ=255)

  atatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaaa                                                                                       >  2:3546970/1‑93 (MQ=255)

   tatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAac                                                                                      >  1:7743249/1‑93 (MQ=255)

   tatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAac                                                                                      <  2:753394/93‑1 (MQ=255)

       cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  >  2:6611533/1‑93 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                <  2:7044492/93‑1 (MQ=255)

         gTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                                                <  2:442330/93‑1 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGTATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               <  2:3381109/93‑1 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              <  1:7894494/93‑1 (MQ=255)

            aGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACgg                                                                             >  2:2260515/1‑93 (MQ=255)

             gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGa                                                                            >  2:6433735/1‑93 (MQ=255)

             gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGa                                                                            <  1:7272189/93‑1 (MQ=255)

             gTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGa                                                                            >  1:5977149/1‑93 (MQ=255)

              tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAg                                                                           >  1:6563114/1‑93 (MQ=255)

              tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAg                                                                           <  1:5621581/93‑1 (MQ=255)

              tGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAg                                                                           <  2:6879342/93‑1 (MQ=255)

               gggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGt                                                                          <  2:4273426/93‑1 (MQ=255)

                ggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACTTCGAAACACAGGACGGAGtt                                                                         <  1:2530899/93‑1 (MQ=255)

                ggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                                         >  2:2697371/1‑93 (MQ=255)

                ggTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                                         <  1:3887441/93‑1 (MQ=255)

                  ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                                       >  1:3254556/1‑93 (MQ=255)

                  ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                                       <  2:130040/93‑1 (MQ=255)

                   tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                                      <  2:2841854/93‑1 (MQ=255)

                    aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                     <  1:3327443/93‑1 (MQ=255)

                    aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                     <  1:2948468/93‑1 (MQ=255)

                    aCCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                                     >  1:4325293/1‑93 (MQ=255)

                     ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCa                                                                    <  2:875652/93‑1 (MQ=255)

                       gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                                  >  2:7920019/1‑93 (MQ=255)

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                                                 <  1:8186684/93‑1 (MQ=255)

                         tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                                                >  1:5187088/1‑93 (MQ=255)

                          aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                                               >  1:4011829/1‑93 (MQ=255)

                          aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                                               <  2:8342280/93‑1 (MQ=255)

                           ttCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                                              <  1:5057618/93‑1 (MQ=255)

                           ttCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                                              >  1:3066655/1‑93 (MQ=255)

                           ttCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                                              >  1:2374430/1‑93 (MQ=255)

                             cGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGc                                                            <  1:2742338/93‑1 (MQ=255)

                              gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAgac                                                           <  1:681462/93‑3 (MQ=255)

                               aCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                          >  2:7139850/1‑93 (MQ=255)

                               aCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                          >  1:8122275/1‑93 (MQ=255)

                                ccAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGa                                                         <  1:2856645/93‑1 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        >  2:2438414/1‑93 (MQ=255)

                                 cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                                        >  1:474328/1‑93 (MQ=255)

                                  aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                                       >  2:6742142/1‑93 (MQ=255)

                                   gCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAt                                                      >  2:3017518/1‑93 (MQ=255)

                                    cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                                     <  2:200245/93‑1 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                    >  1:936171/1‑93 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                    >  1:3508777/1‑93 (MQ=255)

                                     ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                                    >  2:5123142/1‑93 (MQ=255)

                                      gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                                   >  2:1168508/1‑93 (MQ=255)

                                       atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAc                                                  <  1:8556999/93‑1 (MQ=255)

                                        tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                 >  2:2094784/1‑93 (MQ=255)

                                        tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                                                 >  1:1418792/1‑93 (MQ=255)

                                           tttAATTTCTCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACCCAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              <  2:6213664/93‑1 (MQ=255)

                                                tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTa                                         >  1:4891566/1‑93 (MQ=255)

                                                   cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                                      >  2:8420113/1‑93 (MQ=255)

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                                     <  2:3909297/93‑1 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   <  2:4735270/93‑1 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   <  2:8084934/93‑1 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCGGAANTTCACCGACTATTNCAACNtn                                  >  1:7528917/1‑92 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  >  2:583217/1‑93 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  >  2:5300831/1‑93 (MQ=255)

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  >  1:452282/1‑93 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                >  2:5772254/1‑93 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTAATTCAGGTTAAGCCGAAATTCCCCGACTATTAAAACCCCGc                                >  1:2186469/1‑93 (MQ=255)

                                                           caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt                              >  2:372619/1‑93 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGCAATTCACCGACTATTACAACCTCGCTtat                            >  1:5856844/1‑93 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat                            <  2:1881051/93‑1 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat                            <  1:4258022/93‑1 (MQ=255)

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          >  1:2054756/1‑93 (MQ=255)

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          >  2:7779267/1‑93 (MQ=255)

                                                                tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc                         <  2:1500364/93‑1 (MQ=255)

                                                                 aCAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACa                        <  2:5922639/93‑1 (MQ=255)

                                                                  cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaa                       <  2:7345659/93‑1 (MQ=255)

                                                                  cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaa                       >  2:4576071/1‑93 (MQ=255)

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      <  1:4987789/93‑1 (MQ=255)

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      >  1:1047421/1‑93 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     <  2:8577865/93‑1 (MQ=255)

                                                                    gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc                     >  1:7134310/1‑93 (MQ=255)

                                                                     tCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                    <  1:8806666/93‑1 (MQ=255)

                                                                     tCGAATGAATGGCTACAACATCAAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg                    <  2:3132506/93‑1 (MQ=255)

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  >  1:5259847/1‑93 (MQ=255)

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  >  1:6545653/1‑93 (MQ=255)

                                                                        aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                 >  1:3592932/1‑93 (MQ=255)

                                                                        aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                 <  1:5357546/93‑1 (MQ=255)

                                                                        aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                 >  2:4967067/1‑93 (MQ=255)

                                                                         atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCggg                >  2:5228486/1‑93 (MQ=255)

                                                                          tgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGa               >  1:1560409/1‑93 (MQ=255)

                                                                          tgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGa               >  2:7591385/1‑93 (MQ=255)

                                                                            aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaaa             <  1:1316561/93‑1 (MQ=255)

                                                                            aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaaa             <  1:7970290/93‑1 (MQ=255)

                                                                             atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt            >  1:6370973/1‑93 (MQ=255)

                                                                             atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt            >  1:7620106/1‑93 (MQ=255)

                                                                             atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAt            <  2:7688736/93‑1 (MQ=255)

                                                                              tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           >  2:6060684/1‑93 (MQ=255)

                                                                              tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           >  2:2389432/1‑93 (MQ=255)

                                                                               ggCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTa          >  2:1780775/1‑93 (MQ=255)

                                                                               ggCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTa          <  2:4059542/93‑1 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       <  1:180481/93‑1 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAATACaa       >  2:3622571/1‑93 (MQ=255)

                                                                                    caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACt     >  2:5635651/1‑93 (MQ=255)

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

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