BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,274	T→C	L537P (CTC→CCC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,274	0	T	C	98.2%	706.5 / ‑1.0	223	L537P (CTC→CCC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (2/2); new base C (110/109); total (112/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 230 total aligned reads displayed.

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCg                                                                                         <  1:2326941/93‑1 (MQ=255)

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCg                                                                                         >  1:485277/1‑93 (MQ=255)

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCg                                                                                         >  1:1062132/1‑93 (MQ=255)

gAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCg                                                                                         >  1:2601799/1‑93 (MQ=255)

 aGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCCCCGTTACTCAGCAACCCgg                                                                                        <  2:4098870/93‑1 (MQ=255)

  gTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCggg                                                                                       >  2:4567472/1‑93 (MQ=255)

    tATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGgcg                                                                                     >  2:4502989/1‑93 (MQ=255)

     aTTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgcgc                                                                                    <  1:3288283/93‑1 (MQ=255)

     aTTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACACGGccgc                                                                                    >  1:181953/1‑93 (MQ=255)

      ttCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   >  2:3626191/1‑93 (MQ=255)

       tCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGcac                                                                                  <  1:4577309/93‑1 (MQ=255)

       tCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGcac                                                                                  >  1:2942048/1‑93 (MQ=255)

          ggTGAAGCCGAAATGCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               <  1:1197171/93‑1 (MQ=255)

           gTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAtta                                                                              >  2:3850037/1‑91 (MQ=255)

            ttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCaa                                                                             >  2:2933291/1‑93 (MQ=255)

            ttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCaa                                                                             <  1:4239372/93‑1 (MQ=255)

            ttAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCaa                                                                             >  1:414699/1‑93 (MQ=255)

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGCGCGCACATCAacac                                                                          >  1:4002550/1‑93 (MQ=255)

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAacac                                                                          <  1:1539239/93‑1 (MQ=255)

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAacac                                                                          <  2:2242690/93‑1 (MQ=255)

               aGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAacac                                                                          >  2:1536629/1‑93 (MQ=255)

                  cGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGt                                                                       <  1:41931/93‑1 (MQ=255)

                  cGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGt                                                                       >  2:1321683/1‑93 (MQ=255)

                   gAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTa                                                                      >  2:2576258/1‑93 (MQ=255)

                   gAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTa                                                                      <  2:593145/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACGCTGTAtt                                                                    <  2:3279386/93‑1 (MQ=255)

                     aaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAtt                                                                    <  2:858862/93‑1 (MQ=255)

                      aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   <  2:2492247/93‑1 (MQ=255)

                      aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   >  1:3208056/1‑93 (MQ=255)

                      aTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAttt                                                                   >  1:4405241/1‑93 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  >  2:441527/1‑93 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  >  2:2716503/1‑93 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  >  2:761045/1‑93 (MQ=255)

                        tCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGa                                                                 <  1:483518/93‑1 (MQ=255)

                         cACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAg                                                                <  2:2611238/93‑1 (MQ=255)

                           ccGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTg                                                              <  1:4709413/93‑1 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             >  1:4027508/1‑93 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             >  2:1464505/1‑93 (MQ=255)

                            cGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTgg                                                             >  1:1709347/1‑93 (MQ=255)

                              tcTATTGCAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           <  1:1507369/92‑1 (MQ=255)

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           >  2:460654/1‑93 (MQ=255)

                                tATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGc                                                         >  1:1389165/1‑93 (MQ=255)

                                 tttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        <  2:2766778/92‑1 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        <  1:3505408/93‑1 (MQ=255)

                                   tACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCAACGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCAt                                                      <  1:1418037/93‑1 (MQ=255)

                                     cAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCa                                                    >  2:4643347/1‑93 (MQ=255)

                                      aaCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaa                                                   >  2:3083828/1‑93 (MQ=255)

                                      aaCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaa                                                   >  1:2581796/1‑93 (MQ=255)

                                       aCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaaa                                                  >  2:4055170/1‑93 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                >  2:1848152/1‑93 (MQ=255)

                                         cTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACt                                                >  1:2116384/1‑93 (MQ=255)

                                          tCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACtt                                               >  1:4144620/1‑93 (MQ=255)

                                           cGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTa                                              >  2:1140938/1‑93 (MQ=255)

                                           cGCTCATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTa                                              <  2:4166439/93‑1 (MQ=255)

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            >  1:1023614/1‑93 (MQ=255)

                                             cTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAtt                                            <  2:1268158/93‑1 (MQ=255)

                                               tataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                          <  2:411431/93‑1 (MQ=255)

                                               tataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                          <  2:3464309/93‑1 (MQ=255)

                                               tataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgg                                          >  1:2811953/1‑93 (MQ=255)

                                                ataACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTggg                                         <  2:4138460/93‑1 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        >  1:750189/1‑93 (MQ=255)

                                                 taACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTgggg                                        >  1:1791052/1‑93 (MQ=255)

                                                    cAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAAc                                     <  2:4020234/93‑1 (MQ=255)

                                                     aaaCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACg                                    <  1:4089008/93‑1 (MQ=255)

                                                       aCGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAg                                  <  1:3705696/93‑1 (MQ=255)

                                                         gcgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTa                                <  2:3216480/93‑1 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  1:4116972/93‑1 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               >  2:2445253/1‑93 (MQ=255)

                                                          cgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTaa                               <  1:109337/93‑1 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              >  1:1001399/1‑93 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              >  2:128208/1‑93 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              >  2:3351234/1‑93 (MQ=255)

                                                             gAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGt                            <  1:2172114/93‑1 (MQ=255)

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        >  2:572974/1‑93 (MQ=255)

                                                                 ttACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        >  1:3851009/1‑93 (MQ=255)

                                                                    cAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAg                     <  1:3096731/93‑1 (MQ=255)

                                                                     aaCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAga                    >  2:1345045/1‑92 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGCAGCCATCAAACTTATTGGGGAATGAGCAATGTCGATGAgta                   >  1:4471564/1‑91 (MQ=255)

                                                                       cTCACCGATACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCaa                  <  1:2254846/93‑1 (MQ=255)

                                                                        tCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAt                 >  1:1186014/1‑93 (MQ=255)

                                                                         cACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAtt                <  1:2655256/93‑1 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              >  1:454415/1‑93 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              <  1:1272156/93‑1 (MQ=255)

                                                                           ccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTcc              >  2:2177042/1‑93 (MQ=255)

                                                                            cGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCa             >  1:4452967/1‑93 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGCGTAATGTCGATGAGCAATTCCAg            >  2:1024978/1‑93 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            >  1:1772681/1‑93 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            <  2:51126/93‑1 (MQ=255)

                                                                             gTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAg            <  2:2477411/93‑1 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          >  2:937076/1‑93 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          >  1:482071/1‑93 (MQ=255)

                                                                                aCTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCt         <  1:4692707/93‑1 (MQ=255)

                                                                                aCTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCt         <  1:3531891/93‑1 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        <  1:4239508/93‑1 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        <  1:1209923/93‑1 (MQ=255)

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       >  1:1642495/1‑93 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACGTCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      <  2:2038381/93‑1 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      <  1:275767/93‑1 (MQ=255)

                                                                                    aGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt     <  1:65474/93‑1 (MQ=255)

                                                                                       aaCCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTaa  <  2:2978847/93‑1 (MQ=255)

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

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