BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,236,786	A→G	Y525C (TAC→TGC)	yjbH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,236,786	0	A	G	99.7%	1248.4 / ‑5.4	387	Y525C (TAC→TGC)	yjbH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (173/212); total (174/213)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 405 total aligned reads displayed.

TGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATG  >  NC_000913/4236696‑4236875

tgCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCg                                                                                         >  1:7055429/1‑93 (MQ=255)

 gcgcGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCgg                                                                                        <  2:5356527/93‑1 (MQ=255)

 gcgcGAGTATGTGCAGAACGATGTCTAAGTTAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCgg                                                                                        <  1:8006654/93‑1 (MQ=255)

  cgcgAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCggt                                                                                       <  1:9779893/93‑1 (MQ=255)

   gcgAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCggtg                                                                                      >  1:455806/1‑93 (MQ=255)

     gAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCggtggt                                                                                    <  1:2464180/93‑1 (MQ=255)

     gAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACggtggt                                                                                    >  1:9269322/1‑93 (MQ=255)

      aGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGtt                                                                                   >  1:4340049/1‑93 (MQ=255)

       gTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTa                                                                                  <  1:9447465/93‑1 (MQ=255)

        tATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAt                                                                                 >  2:7793933/1‑93 (MQ=255)

          tgtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACTGCTTCTACGGTCAGGTCTGCGGTGGTTAtct                                                                               >  1:4920101/1‑93 (MQ=255)

          tgtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtct                                                                               <  2:2294240/93‑1 (MQ=255)

          tgtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtct                                                                               <  1:611340/93‑1 (MQ=255)

           gtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtctc                                                                              <  1:8476232/93‑1 (MQ=255)

           gtgCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtctc                                                                              <  2:2722834/93‑1 (MQ=255)

                  aCGATGTGTATGTGAATAACATGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCa                                                                       <  1:2603767/93‑1 (MQ=255)

                  aCGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCa                                                                       <  1:8997664/93‑1 (MQ=255)

                  aCGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCa                                                                       <  1:6417516/93‑1 (MQ=255)

                   cGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCAt                                                                      <  2:9809441/93‑1 (MQ=255)

                     aTGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGt                                                                    >  1:3147881/1‑93 (MQ=255)

                     aTGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGt                                                                    <  2:6779270/93‑1 (MQ=255)

                     aTGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGt                                                                    >  1:1566917/1‑93 (MQ=255)

                     aTGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGt                                                                    >  1:2119883/1‑93 (MQ=255)

                      tGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGtt                                                                   <  2:9505776/93‑1 (MQ=255)

                      tGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGtt                                                                   >  2:8162703/1‑93 (MQ=255)

                      tGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGtt                                                                   <  2:7566323/93‑1 (MQ=255)

                       gTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGttt                                                                  <  2:3820824/93‑1 (MQ=255)

                         cTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTgg                                                                <  2:2941921/93‑1 (MQ=255)

                           aTGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCg                                                              <  2:6402004/93‑1 (MQ=255)

                            tgtgAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCgg                                                             <  1:5048109/93‑1 (MQ=255)

                            tgtgAATAACCTGCAAGCCAACTACTGCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCgg                                                             <  1:3024148/93‑1 (MQ=255)

                             gtgAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGt                                                            <  2:7342688/93‑1 (MQ=255)

                             gtgAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGt                                                            <  1:6606166/93‑1 (MQ=255)

                               caaTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGc                                                          <  2:5171615/92‑1 (MQ=255)

                                aaTAACCTGCAAGCCAACTACTTCCAGCGTCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCg                                                         >  1:9098237/1‑93 (MQ=255)

                                aaTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCg                                                         <  2:3575396/93‑1 (MQ=255)

                                  tAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggg                                                       >  2:4474509/1‑93 (MQ=255)

                                   aaCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCgggg                                                      <  1:7758770/93‑1 (MQ=255)

                                   aaCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCgggg                                                      >  2:210569/1‑93 (MQ=255)

                                    aCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggggg                                                     >  1:7224441/1‑93 (MQ=255)

                                      cTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCa                                                   <  2:3453399/93‑1 (MQ=255)

                                           aGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgt                                              >  2:3871061/1‑93 (MQ=255)

                                           aGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgg                                              >  2:3151931/1‑92 (MQ=255)

                                           aGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgg                                              >  1:1020404/1‑92 (MQ=255)

                                            gCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgtg                                             >  1:7551563/1‑93 (MQ=255)

                                             ccAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTTGCGGTGCGGGGGCAGAAgtgt                                            >  2:5451730/1‑93 (MQ=255)

                                              cAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgtt                                           >  1:5684445/1‑93 (MQ=255)

                                                actactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgttgt                                         >  2:2452954/1‑93 (MQ=255)

                                                actactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgttgt                                         <  2:3301050/93‑1 (MQ=255)

                                                 ctactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTa                                        <  1:3294877/93‑1 (MQ=255)

                                                  tactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTAt                                       <  1:1471933/93‑1 (MQ=255)

                                                  tactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTAt                                       <  1:5329818/93‑1 (MQ=255)

                                                   actTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATc                                      <  2:6254669/93‑1 (MQ=255)

                                                     ttCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGc                                    >  1:1081080/1‑93 (MQ=255)

                                                     ttCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGc                                    >  1:5439134/1‑93 (MQ=255)

                                                       ccAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGccc                                  <  1:4025355/93‑1 (MQ=255)

                                                          gCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCt                               >  1:5895819/1‑93 (MQ=255)

                                                          gCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCt                               <  2:6176702/93‑1 (MQ=255)

                                                           gaTCTGGGCAACGGCGTCGACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGCGTTGTATCGCCCGCTg                              <  1:7174915/92‑1 (MQ=255)

                                                            aTCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTgg                             >  1:3865483/1‑93 (MQ=255)

                                                             tCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGa                            >  2:5810481/1‑93 (MQ=255)

                                                              cTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGAc                           <  1:1166343/93‑1 (MQ=255)

                                                               tGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACa                          >  1:5810695/1‑93 (MQ=255)

                                                               tGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACa                          <  2:9347284/93‑1 (MQ=255)

                                                                gggCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCCGAAGTGTTGTATCGCCCGCTGGACAg                         >  1:1913378/1‑93 (MQ=255)

                                                                 ggCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGc                        >  1:521601/1‑93 (MQ=255)

                                                                  gCAACGACATCTGCGGTCAGTTCTTCGGTGGTAATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCa                       <  1:1787353/93‑1 (MQ=255)

                                                                   cAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCaa                      <  2:1582048/93‑1 (MQ=255)

                                                                   cAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCaa                      <  2:8368720/93‑1 (MQ=255)

                                                                    aaaggCTTTTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAAc                     <  2:2003557/90‑1 (MQ=255)

                                                                    aaCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAAc                     <  1:922323/93‑1 (MQ=255)

                                                                    aaCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAAc                     <  1:7825035/93‑1 (MQ=255)

                                                                    aaCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAAc                     <  1:9574696/93‑1 (MQ=255)

                                                                     aCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACt                    >  2:412122/1‑93 (MQ=255)

                                                                     aCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACt                    >  2:4306848/1‑93 (MQ=255)

                                                                     aCGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACt                    >  1:5148078/1‑93 (MQ=255)

                                                                         cTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGc                <  1:1264624/93‑1 (MQ=255)

                                                                         cTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGc                >  2:8896577/1‑93 (MQ=255)

                                                                             tACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGttt            <  1:3679212/93‑1 (MQ=255)

                                                                             tACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGttt            <  2:10033713/93‑1 (MQ=255)

                                                                              tcgctCAGGTCTGCGGTTGTTATCTCGAAGCCATGTTTGGCGGTGCGGGGGAAGAAGTGTTGTATCGCCCGCTGGACAGCTACTGGGCGTTTg           <  2:1813450/89‑1 (MQ=255)

                                                                              aCGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTg           <  2:9075851/93‑1 (MQ=255)

                                                                              aCGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTg           <  2:1278114/93‑1 (MQ=255)

                                                                                ggTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGt         <  1:2736474/93‑1 (MQ=255)

                                                                                ggTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGt         <  2:6947041/93‑1 (MQ=255)

                                                                                ggTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGt         >  1:2023586/1‑93 (MQ=255)

                                                                                 gTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTc        <  2:6576211/93‑1 (MQ=255)

                                                                                 gTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTc        <  1:133107/93‑1 (MQ=255)

                                                                                 gTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTc        >  2:8013746/1‑93 (MQ=255)

                                                                                  tCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCt       >  2:2818638/1‑93 (MQ=255)

                                                                                   cAGGTCGTCGGTGGTTAGCTCGAAACCATGTTTGGCTGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTg      <  1:8371585/93‑1 (MQ=255)

                                                                                    aGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGGCAGCAACTGGGCGTTTGGTCTgg     >  2:5504951/1‑93 (MQ=255)

                                                                                    aGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTgg     <  1:8744541/93‑1 (MQ=255)

                                                                                    aGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTgg     <  1:1621713/93‑1 (MQ=255)

                                                                                    aGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTgg     >  2:8631254/1‑93 (MQ=255)

                                                                                    aGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTgg     <  2:3387227/93‑1 (MQ=255)

                                                                                      gTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGAt   <  1:9863318/93‑1 (MQ=255)

                                                                                      gTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGAt   <  1:2890914/93‑1 (MQ=255)

                                                                                       tCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATg  >  2:7041144/1‑93 (MQ=255)

                                                                                       tCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATg  <  2:4799196/93‑1 (MQ=255)

TGCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATG  >  NC_000913/4236696‑4236875

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