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	99.6%	899.7 / NA	254	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (142/111); minor base A (0/1); total (142/112)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGa                                                                                         >  1:785790/1‑93 (MQ=255)

gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGa                                                                                         <  2:182893/93‑1 (MQ=255)

gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGa                                                                                         >  1:2438081/1‑93 (MQ=255)

gAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGa                                                                                         >  1:2102219/1‑93 (MQ=255)

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

  atatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaaa                                                                                       <  1:1021842/93‑1 (MQ=255)

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

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

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

   tatTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGacac                                                                                      >  1:3207514/1‑93 (MQ=255)

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

    atTGAGTTAGTGGGTTACCGTTATGCGACCAGCGGATATTTTAAGTTCGCTGATACAACATACAGGCGAATGAATGGCTACAACATCGAAaca                                                                                     <  2:3189794/93‑1 (MQ=255)

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

        aGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAgg                                                                                 <  2:2760457/93‑1 (MQ=255)

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

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

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               >  2:5423112/1‑93 (MQ=255)

          ttAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGAc                                                                               <  1:3396171/93‑1 (MQ=255)

           tAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                                              >  1:2644408/1‑93 (MQ=255)

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

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

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

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

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

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

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                                                 >  1:5357869/1‑93 (MQ=255)

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

                        ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                                                 >  1:1999579/1‑93 (MQ=255)

                         tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                                                <  1:5589962/93‑1 (MQ=255)

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

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

                            tCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAg                                                             >  2:5717259/1‑93 (MQ=255)

                            tCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGATATTCAGGTTAAg                                                             >  2:5011536/1‑93 (MQ=255)

                            tCGACCAGCGGATATTTTAATTTCCCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAg                                                             >  1:5639590/1‑93 (MQ=255)

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

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

                              gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                                           >  2:3634077/1‑93 (MQ=255)

                                       atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAc                                                  >  1:1933372/1‑93 (MQ=255)

                                         atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                                                <  2:4816114/93‑1 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              >  1:119943/1‑93 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              >  1:2466652/1‑93 (MQ=255)

                                               aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                          >  1:4241386/1‑93 (MQ=255)

                                               aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                                          <  2:1918606/93‑1 (MQ=255)

                                                tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTa                                         <  1:3127970/93‑1 (MQ=255)

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

                                                 ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                        <  2:769110/93‑1 (MQ=255)

                                                 ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                        >  1:4906619/1‑93 (MQ=255)

                                                 ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                                        <  1:4350427/93‑1 (MQ=255)

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

                                                   cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                                      <  1:1745786/93‑1 (MQ=255)

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

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

                                                     cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAcc                                    >  2:2096139/1‑93 (MQ=255)

                                                      tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                                   >  2:3351938/1‑93 (MQ=255)

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

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

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  2:1555086/93‑1 (MQ=255)

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

                                                       gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                                  <  1:422280/93‑1 (MQ=255)

                                                         tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc                                <  1:2056555/93‑1 (MQ=255)

                                                          acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCt                               >  2:2492631/1‑93 (MQ=255)

                                                           caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt                              <  2:4957765/93‑1 (MQ=255)

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

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

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             <  2:1305629/93‑1 (MQ=255)

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             >  1:3301720/1‑93 (MQ=255)

                                                            aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta                             >  1:2809560/1‑93 (MQ=255)

                                                             acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat                            >  2:1500665/1‑93 (MQ=255)

                                                              caTACAGACGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata                           <  2:756037/93‑1 (MQ=255)

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

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

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

                                                                   aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGCTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa                      >  1:5560748/1‑93 (MQ=255)

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

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

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

                                                                        aatgaatgGCTACAACATCGAAACGCAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCgg                 <  1:2751160/93‑1 (MQ=255)

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

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

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

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

                                                                           gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGaa              >  2:2341255/1‑93 (MQ=255)

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

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

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

                                                                              tgGCTACAACATAGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           <  2:4791411/93‑1 (MQ=255)

                                                                               ggCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACTACCTCGCTTATAACAAACGCGGGAAATTa          >  2:1122773/1‑93 (MQ=255)

                                                                                gCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc         >  1:1451736/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  1:2223753/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  2:3095928/1‑93 (MQ=255)

                                                                                  tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACaa       >  1:5762898/1‑93 (MQ=255)

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

                                                                                     aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc    >  2:3836430/1‑93 (MQ=255)

                                                                                     aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc    >  1:173357/1‑93 (MQ=255)

                                                                                     aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc    <  2:6295871/93‑1 (MQ=255)

                                                                                     aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTc    <  1:5740949/93‑1 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  2:287698/1‑93 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  2:263899/1‑93 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  1:6336255/1‑93 (MQ=255)

                                                                                      acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCa   >  1:5121365/1‑93 (MQ=255)

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

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