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	96.8%	572.2 / 3.7	186	Y525C (TAC→TGC)	yjbH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/4); new base G (91/89); total (93/93)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 189 total aligned reads displayed.

GCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATG  >  NC_000913/4236697‑4236875

gcgcGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACgg                                                                                        <  2:1551024/93‑1 (MQ=255)

  gcgAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCggtg                                                                                      <  1:2155779/93‑1 (MQ=255)

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

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

        aTGTGCAGAACGATGTGTATGTGGATAACCTTCAAGCCAACTACTCCCAGCATCTGGGCAACGGCTTCTACCTTCAGGTCTGCGGTGGTTaac                                                                                <  2:1657056/93‑3 (MQ=255)

        aTGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtc                                                                                <  1:1087111/93‑1 (MQ=255)

        aTGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtc                                                                                <  1:1493485/93‑1 (MQ=255)

        aTGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAtc                                                                                >  1:2557328/1‑93 (MQ=255)

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

              aGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGaac                                                                          >  1:746833/1‑92 (MQ=255)

              aGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGaaa                                                                          <  1:2572081/93‑1 (MQ=255)

               gAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAAc                                                                         <  2:1952367/93‑1 (MQ=255)

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

                 aCGATGTCTATGTGAATAACCTGAAAGCAAACTACTTCCAGAATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAACTCGAAACCa                                                                       <  1:4196468/93‑1 (MQ=255)

                  cGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCAt                                                                      >  2:937357/1‑93 (MQ=255)

                  cGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCAt                                                                      >  2:4426153/1‑93 (MQ=255)

                   gATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATg                                                                     >  2:3462718/1‑93 (MQ=255)

                   gATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATg                                                                     >  1:3879771/1‑93 (MQ=255)

                   gATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATg                                                                     >  1:327962/1‑93 (MQ=255)

                   gATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGCCTGCGGTGGTTATCTCGAAACCATg                                                                     <  1:333381/93‑1 (MQ=255)

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

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

                       tCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTg                                                                 <  2:3566433/93‑1 (MQ=255)

                           tgtgAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCgg                                                             >  1:3833901/1‑93 (MQ=255)

                            gtgAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCggg                                                            >  1:2507427/1‑92 (MQ=255)

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

                              taaTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGc                                                          <  1:1173888/92‑1 (MQ=255)

                              gAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGc                                                          <  1:856172/93‑1 (MQ=255)

                              gAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGc                                                          <  1:4534548/93‑1 (MQ=255)

                                aTACCCTGCAAGCCAACTCCTTCCAGCCGCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGTCGGTGCgg                                                        <  1:4636025/93‑1 (MQ=255)

                                aTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCgg                                                        <  2:4743140/93‑1 (MQ=255)

                                aTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCgg                                                        <  1:4636024/93‑1 (MQ=255)

                                aTAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCgg                                                        >  2:2132339/1‑93 (MQ=255)

                                 tAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGgcggggcggg                                                       >  2:1870168/1‑93 (MQ=255)

                                 tAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggg                                                       <  1:3467758/93‑1 (MQ=255)

                                 tAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggg                                                       <  1:2335044/93‑1 (MQ=255)

                                 aaaCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggg                                                       <  1:3167456/92‑1 (MQ=255)

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

                                   aCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggggg                                                     <  1:3524089/93‑1 (MQ=255)

                                   aCCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCggggg                                                     <  1:4055984/93‑1 (MQ=255)

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

                                    ccTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGc                                                    <  1:3440566/93‑1 (MQ=255)

                                    ccTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGc                                                    <  2:1969707/93‑1 (MQ=255)

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

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

                                      tGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAg                                                  >  1:4206546/1‑93 (MQ=255)

                                      tGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAg                                                  >  2:1077990/1‑93 (MQ=255)

                                        cAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGaa                                                >  2:2031527/1‑93 (MQ=255)

                                        cAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGaa                                                <  1:2014302/93‑1 (MQ=255)

                                         aaGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAg                                               <  1:2640293/93‑1 (MQ=255)

                                         aaGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAg                                               >  2:4507636/1‑93 (MQ=255)

                                         aaGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAg                                               <  1:1428123/93‑1 (MQ=255)

                                           gCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgtg                                             <  1:1577759/93‑1 (MQ=255)

                                            ccACCTACTTCGAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTAGCTCGAAACCATGTTTGGCAGTGCGGGGGCAGAAgtct                                            <  2:4407248/93‑3 (MQ=255)

                                            ccAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAgtgt                                            <  1:4615739/93‑1 (MQ=255)

                                             cAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgtt                                           <  1:4714953/93‑1 (MQ=255)

                                             cAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgtt                                           <  1:138940/93‑1 (MQ=255)

                                             cAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTAAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgtt                                           >  2:1212487/1‑93 (MQ=255)

                                              aaCTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgttg                                          >  1:2287875/1‑93 (MQ=255)

                                              aaCTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGtgttg                                          >  1:3762442/1‑93 (MQ=255)

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

                                                 tactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTAt                                       >  1:1006817/1‑93 (MQ=255)

                                                 tactTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTAt                                       >  1:4317346/1‑93 (MQ=255)

                                                   ctTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATcc                                     >  2:3453110/1‑92 (MQ=255)

                                                    ttCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGc                                    <  2:4177403/93‑1 (MQ=255)

                                                    ttCCAGCATCTGGGCAACGGCTTATACGGTCAGGTCTGCAGCGGTTATCTCGAAACCATGCTTGGCGGTGGGGGGGCAGAAGCGTTGTATCGc                                    >  2:1611148/1‑93 (MQ=255)

                                                     tCCAGCATCTGGGCCACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGcc                                   <  1:860273/93‑1 (MQ=255)

                                                      ccAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGccc                                  >  2:1371525/1‑93 (MQ=255)

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

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

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

                                                               gggCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAg                         >  2:4647362/1‑93 (MQ=255)

                                                                ggCAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGc                        <  1:1291833/93‑1 (MQ=255)

                                                                  cAACGGCTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAGCCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCaa                      <  2:112621/93‑1 (MQ=255)

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

                                                                    aCGGCTTCTACGGTCAGGTCTGCGGTGGTTAACTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGCATCGCCCGCTGGACAGCAACt                    >  2:562881/1‑93 (MQ=255)

                                                                     cgtcTTCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTg                   <  1:2087660/90‑1 (MQ=255)

                                                                          tCTACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGt              >  1:4738624/1‑93 (MQ=255)

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

                                                                            tACGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGttt            >  2:4019531/1‑93 (MQ=255)

                                                                             aCGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTGTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTg           <  2:261890/93‑1 (MQ=255)

                                                                             aCGGTCAGGTCTGAGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTg           <  1:4235063/93‑1 (MQ=255)

                                                                              cGGTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGTGCGTTTgg          >  2:2147317/1‑93 (MQ=255)

                                                                                gTCAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGTGCGCTTGGTc        >  2:3179099/1‑93 (MQ=255)

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

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

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

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

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

                                                                                  cAGGTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTg      <  2:1698424/93‑1 (MQ=255)

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

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

                                                                                    ggTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTggt    >  2:2367659/1‑92 (MQ=255)

                                                                                    ggTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGa    <  2:4568131/93‑1 (MQ=255)

                                                                                    ggTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGa    >  1:2723591/1‑93 (MQ=255)

                                                                                    ggTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGCCTGGa    >  2:1346822/1‑93 (MQ=255)

                                                                                     gTCTGCGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGAt   >  1:783853/1‑93 (MQ=255)

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

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

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

GCGCGAGTATGTGCAGAACGATGTCTATGTGAATAACCTGCAAGCCAACTACTTCCAGCATCTGGGCAACGGCTTCTACGGTCAGGTCTACGGTGGTTATCTCGAAACCATGTTTGGCGGTGCGGGGGCAGAAGTGTTGTATCGCCCGCTGGACAGCAACTGGGCGTTTGGTCTGGATG  >  NC_000913/4236697‑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: `‑`