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.7%	1081.0 / NA	316	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (169/146); minor base G (0/1); total (169/147)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 316 total aligned reads displayed.

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

    atTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                                     <  2:5796959/93‑1 (MQ=255)

     ttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                                    >  2:4937833/1‑93 (MQ=255)

     ttCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                                    <  2:6678379/93‑1 (MQ=255)

      tCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                                   <  2:7179085/93‑1 (MQ=255)

       cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  <  1:5519521/93‑1 (MQ=255)

       cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  <  2:7035161/93‑1 (MQ=255)

       cAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                                  <  1:6787899/93‑1 (MQ=255)

        aGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAgg                                                                                 >  1:2670950/1‑93 (MQ=255)

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

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

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

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

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

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

                 gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTa                                                                        >  2:5221913/1‑93 (MQ=255)

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

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

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

                   tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                                      >  2:289562/1‑93 (MQ=255)

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

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

                     ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCa                                                                    >  2:2041763/1‑93 (MQ=255)

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

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

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

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

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

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

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

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

                             cGACCAGCGGATAATTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATGGAAACACAGGACGGAGTTATTCAGGTTAAGc                                                            <  2:1207410/93‑1 (MQ=255)

                              gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                                           <  2:2083652/93‑1 (MQ=255)

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

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

                               aCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                                          <  1:1242882/93‑1 (MQ=255)

                                  aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                                       <  1:4934525/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                          ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                                               <  2:7365886/93‑1 (MQ=255)

                                           tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                                              <  2:1762234/93‑1 (MQ=255)

                                            ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                                             >  1:5470099/1‑93 (MQ=255)

                                             tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTa                                            >  2:6103474/1‑93 (MQ=255)

                                              aaTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                                           <  2:2937203/93‑1 (MQ=255)

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

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

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCTCCGACTATTACAAc                                     <  2:6453718/93‑1 (MQ=255)

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

                                                    gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAATTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                                     <  1:6939346/93‑1 (MQ=255)

                                                     cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAcc                                    <  1:2826589/93‑1 (MQ=255)

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

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

                                                        aTACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg                                 >  2:4971395/1‑93 (MQ=255)

                                                        aTACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg                                 >  2:4798943/1‑93 (MQ=255)

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

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

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

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

                                                               aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATaa                          <  1:99541/93‑1 (MQ=255)

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

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

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

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

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

                                                                       gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcgcg                  <  1:3234455/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                              tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           <  1:206107/93‑1 (MQ=255)

                                                                              tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAAtt           <  1:7036774/93‑1 (MQ=255)

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

                                                                                gCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTAc         <  2:2701298/93‑1 (MQ=255)

                                                                                 cTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACa        <  2:2318031/93‑1 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      <  1:1497213/93‑1 (MQ=255)

                                                                                   acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAAc      <  1:4090610/93‑1 (MQ=255)

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

                                                                                     aacaTCGAAACACAGGACGGGGGTATTTCGGTTAAGGCGGAATTTACCGACCATTACAACCTCGGCTATTACAAACGCGGGGAATTAAAACTc    >  1:2376797/1‑93 (MQ=255)

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

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

                                                                                       caTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAc  >  2:6361509/1‑93 (MQ=255)

GAATATTCAGTTAGTGGGTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAC  >  NC_000913/4544080‑4544259

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