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	99.8%	1821.3 / ‑5.6	563	L537P (CTC→CCC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (247/315); total (248/315)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 574 total aligned reads displayed.

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

gAGTTATTCAGGTTAAGCCTAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTCCTCAGCAACCCg                                                                                         <  2:5204015/93‑1 (MQ=255)

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

 aGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCgg                                                                                        <  1:1572375/93‑1 (MQ=255)

  gTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCggg                                                                                       <  1:1909193/93‑1 (MQ=255)

   ttATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGgc                                                                                      >  1:8470117/1‑93 (MQ=255)

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

      ttCAGGTTAAGCCGAAATTCACCGACTAGTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGca                                                                                   <  1:217950/93‑1 (MQ=255)

         aGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACAt                                                                                <  1:888253/93‑1 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               >  2:5367234/1‑93 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               >  1:9080831/1‑93 (MQ=255)

          ggTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATc                                                                               >  2:2934636/1‑93 (MQ=255)

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

             tAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAac                                                                            >  1:6324116/1‑93 (MQ=255)

             tAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAac                                                                            <  2:2734699/93‑1 (MQ=255)

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

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

                gCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACt                                                                         >  2:262448/1‑93 (MQ=255)

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

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

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

                    aaaTTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTAt                                                                     >  2:1758674/1‑93 (MQ=255)

                       ttCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTg                                                                  <  1:6595240/93‑1 (MQ=255)

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

                        tCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGa                                                                 >  2:6766895/1‑93 (MQ=255)

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

                          cccGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCCCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               <  1:1783337/92‑1 (MQ=255)

                          aCCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGt                                                               <  2:5451771/93‑1 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            <  1:2018724/93‑1 (MQ=255)

                             gACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGt                                                            <  2:4674300/93‑1 (MQ=255)

                              aCTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTa                                                           <  2:7905549/93‑1 (MQ=255)

                               cTATTACACCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAg                                                          <  2:3894192/93‑1 (MQ=255)

                               cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAg                                                          >  1:7926310/1‑93 (MQ=255)

                               cTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAg                                                          <  1:7509258/93‑1 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  1:9568493/1‑93 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  2:3080837/1‑93 (MQ=255)

                                 aTTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGcc                                                        >  2:824155/1‑93 (MQ=255)

                                  ttACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCa                                                       <  2:6163267/93‑1 (MQ=255)

                                    aCAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATc                                                     >  2:5057439/1‑93 (MQ=255)

                                     cAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCa                                                    <  1:3826869/93‑1 (MQ=255)

                                      aaCCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCaa                                                   <  2:1456640/93‑1 (MQ=255)

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

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

                                           cGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTa                                              <  2:1381528/93‑1 (MQ=255)

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

                                            gCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAt                                             <  2:3980993/93‑1 (MQ=255)

                                            gCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTAt                                             <  2:5562599/93‑1 (MQ=255)

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

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

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

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

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

                                                  aaCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGa                                       >  1:4797120/1‑93 (MQ=255)

                                                   aCAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGaa                                      <  1:5961585/93‑1 (MQ=255)

                                                    cAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAAc                                     >  1:3642244/1‑93 (MQ=255)

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

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

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

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

                                                        cgcgGGAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGt                                 >  2:3472796/1‑93 (MQ=255)

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

                                                           gggAAATTACAACTCGCCGTTACACAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              <  1:8562145/93‑1 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCACACTTATTGGGGAACGAGTCAt                              >  1:5444806/1‑93 (MQ=255)

                                                           gggAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAAt                              <  2:2563840/93‑1 (MQ=255)

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

                                                             gAAATTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGt                            >  1:2594532/1‑93 (MQ=255)

                                                              aaaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           <  1:9283786/93‑1 (MQ=255)

                                                              aaaTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTc                           <  2:3716244/93‑1 (MQ=255)

                                                                aTTACACCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  2:92995/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  2:6402085/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         >  2:1144504/1‑93 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  2:5747186/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  1:3195588/93‑1 (MQ=255)

                                                                aTTACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGa                         <  2:4597376/93‑1 (MQ=255)

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

                                                                 ttACAACTCAACGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGAt                        >  2:883566/1‑93 (MQ=255)

                                                                  tACAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATg                       <  1:9897289/93‑1 (MQ=255)

                                                                   aCAACTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGa                      <  1:2438830/93‑1 (MQ=255)

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

                                                                     aaCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                    <  2:9481534/93‑1 (MQ=255)

                                                                     aaCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGc                    <  1:5557637/93‑1 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                   >  2:4387301/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                   >  1:530703/1‑93 (MQ=255)

                                                                      aCTCACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCa                   >  2:1647864/1‑93 (MQ=255)

                                                                        cccccGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAt                 <  1:8114959/90‑1 (MQ=255)

                                                                         cACCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAAtt                >  1:7084770/1‑93 (MQ=255)

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

                                                                          aCCGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTc               >  2:9051066/1‑93 (MQ=255)

                                                                            cGTTACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCa             <  2:7247046/93‑1 (MQ=255)

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

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

                                                                              ttACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAgg           >  2:4861612/1‑93 (MQ=255)

                                                                               tACTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGc          <  2:6152086/93‑1 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  2:7504063/1‑93 (MQ=255)

                                                                                 cTCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTg        >  2:8565647/1‑93 (MQ=255)

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

                                                                                  tCAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTgg       <  1:281756/93‑1 (MQ=255)

                                                                                   cAGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGa      >  2:8422887/1‑93 (MQ=255)

                                                                                    aGCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAt     >  1:121812/1‑93 (MQ=255)

                                                                                     gCAACCCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGAtt    <  2:3363305/93‑1 (MQ=255)

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

GAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACGCGGGAAATTACAACTCACCGTTACTCAGCAACTCGGGCGCACATCAACACTGTATTTGAGTGGTAGCCATCAAACTTATTGGGGAACGAGTAATGTCGATGAGCAATTCCAGGCTGGATTAA  >  NC_000913/4544184‑4544363

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