BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,100,603	A→G	F90L (TTC→CTC)	yggH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,100,603	0	A	G	98.9%	554.6 / ‑2.7	178	F90L (TTC→CTC)	yggH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (87/87); total (89/89)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 186 total aligned reads displayed.

TGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAA  >  NC_000913/3100512‑3100692

tGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAAAACCCGGTGAATACACTTCAATGCCGAGGAAg                                                                                          >  1:477436/1‑93 (MQ=255)

   cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCCCCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTct                                                                                       <  2:2705155/93‑2 (MQ=255)

   cacaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTcc                                                                                       <  2:685657/93‑1 (MQ=255)

    acaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCt                                                                                      <  1:4739638/93‑1 (MQ=255)

    acaTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCt                                                                                      <  1:326404/93‑1 (MQ=255)

     caTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTg                                                                                     >  2:3701335/1‑93 (MQ=255)

     caTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTg                                                                                     <  1:100847/93‑1 (MQ=255)

      aTCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGc                                                                                    >  1:4055842/1‑93 (MQ=255)

        cACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTc                                                                                  <  1:2469159/93‑1 (MQ=255)

         aCGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCa                                                                                 >  1:3692139/1‑93 (MQ=255)

          cgcgCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAg                                                                                <  1:923733/93‑1 (MQ=255)

           gcgcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAgg                                                                               <  1:910561/93‑1 (MQ=255)

            cgcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAggg                                                                              <  2:2011345/93‑1 (MQ=255)

            cgcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAggg                                                                              <  2:421485/93‑1 (MQ=255)

             gcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGc                                                                             >  1:2659876/1‑93 (MQ=255)

             gcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGc                                                                             <  2:4628562/93‑1 (MQ=255)

             gcAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGc                                                                             >  2:540381/1‑93 (MQ=255)

              cAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCg                                                                            >  2:748608/1‑93 (MQ=255)

              cAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCg                                                                            <  1:980927/93‑1 (MQ=255)

                ggTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGAt                                                                          <  2:1822167/93‑1 (MQ=255)

                   tGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCtt                                                                       >  1:1850809/1‑93 (MQ=255)

                    gCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCttt                                                                      <  2:805840/93‑1 (MQ=255)

                     cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     >  1:1389158/1‑93 (MQ=255)

                     cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     <  1:19451/93‑1 (MQ=255)

                     cTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTa                                                                     <  2:2072730/93‑1 (MQ=255)

                      ttAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTac                                                                    >  1:1135061/1‑92 (MQ=255)

                         aaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGATTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCa                                                                 <  2:2811957/93‑1 (MQ=255)

                         aaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCa                                                                 >  2:738622/1‑93 (MQ=255)

                         aaCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCa                                                                 <  1:1877336/93‑1 (MQ=255)

                          aCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAt                                                                >  1:2852264/1‑93 (MQ=255)

                          aCCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAt                                                                <  1:4237383/93‑1 (MQ=255)

                           ccTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAtt                                                               <  1:2272548/93‑1 (MQ=255)

                           ccTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAtt                                                               <  1:1465119/93‑1 (MQ=255)

                           ccTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCAtt                                                               >  2:1136792/1‑93 (MQ=255)

                            cttcttCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTg                                                              >  1:1198899/1‑93 (MQ=255)

                               cttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGcca                                                           <  1:2305304/93‑1 (MQ=255)

                               cttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGcca                                                           >  2:4147255/1‑93 (MQ=255)

                                ttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccac                                                          >  1:2173885/1‑93 (MQ=255)

                                ttcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccac                                                          >  1:130476/1‑93 (MQ=255)

                                 tcATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGccacc                                                         >  1:3674208/1‑93 (MQ=255)

                                    tGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGc                                                      <  2:3166198/93‑1 (MQ=255)

                                     gcgcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCg                                                     <  2:4398599/93‑1 (MQ=255)

                                      cgcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGa                                                    <  1:103298/93‑1 (MQ=255)

                                      cgcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGa                                                    >  2:1873783/1‑93 (MQ=255)

                                      cgcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGa                                                    <  2:533711/93‑1 (MQ=255)

                                       gcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGAc                                                   >  1:4524888/1‑93 (MQ=255)

                                       gcAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGAc                                                   >  1:3233494/1‑93 (MQ=255)

                                        cAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACg                                                  >  2:3897816/1‑93 (MQ=255)

                                          gAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcc                                                <  2:4743960/93‑1 (MQ=255)

                                           aacccAGGGACGCACCAACGCCCGGTGAATGCCCTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGccc                                               <  1:3282382/90‑1 (MQ=255)

                                            agccagccACCCCCCAACCCCCGGTAAATGCATTTCAATCCCGAGGAGGTCCTGCTCAGGCCGATCTTTAGCCATTCCCACCAGCGACGcccc                                              <  2:1595998/93‑1 (MQ=255)

                                            aGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcccc                                              <  1:2983466/93‑1 (MQ=255)

                                            aGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcccc                                              <  2:2432028/93‑1 (MQ=255)

                                            aGCCAGGCACGCACCAACGCCCGATGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGcccc                                              <  1:1090494/93‑1 (MQ=255)

                                              ccAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCa                                            <  1:4291439/93‑1 (MQ=255)

                                              ccAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCa                                            >  1:1610450/1‑93 (MQ=255)

                                                aGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATg                                          >  1:2556740/1‑93 (MQ=255)

                                                 ggCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGc                                         >  1:3430878/1‑93 (MQ=255)

                                                 ggCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGc                                         >  2:208065/1‑93 (MQ=255)

                                                  gcacgcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCAACCAGCGACGCCCCCATGcc                                        >  1:1493729/1‑93 (MQ=255)

                                                   cacgcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCa                                       <  2:4534675/93‑1 (MQ=255)

                                                     cgcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCaaa                                     >  1:921/1‑93 (MQ=255)

                                                      gcacCAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCaaaa                                    <  1:462150/93‑1 (MQ=255)

                                                        cccAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAAcc                                  <  1:1426738/92‑1 (MQ=255)

                                                          cAACGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCaa                                >  2:1694124/1‑93 (MQ=255)

                                                           aaCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAt                               <  1:2115215/93‑1 (MQ=255)

                                                           aaCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAt                               >  2:1052038/1‑93 (MQ=255)

                                                            aCGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtc                              >  1:1729547/1‑93 (MQ=255)

                                                             cGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtct                             <  1:3643181/93‑1 (MQ=255)

                                                             cGCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtct                             <  2:3072557/93‑1 (MQ=255)

                                                              gCCCGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAAtctc                            >  1:2079343/1‑93 (MQ=255)

                                                                ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCTATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCaa                          <  2:4461266/93‑1 (MQ=255)

                                                                ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTaaa                          <  2:3673188/93‑1 (MQ=255)

                                                                ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCaa                          <  2:4686039/93‑1 (MQ=255)

                                                                ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCaa                          >  2:1233773/1‑93 (MQ=255)

                                                                ccGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCaa                          <  2:1573367/93‑1 (MQ=255)

                                                                 cGGTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAg                         >  1:1557425/1‑93 (MQ=255)

                                                                  ggTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCCATCTCAAGc                        >  2:3316506/1‑93 (MQ=255)

                                                                   gTGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCg                       >  1:2205129/1‑93 (MQ=255)

                                                                    tGAATGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGt                      <  2:247317/93‑1 (MQ=255)

                                                                      aaTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCa                    >  2:2117019/1‑93 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCGTTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   >  2:2809215/1‑93 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   >  2:2977616/1‑93 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   <  1:1689794/93‑1 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   <  2:56036/93‑1 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   <  2:3639032/93‑1 (MQ=255)

                                                                       aTGCACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCAc                   >  2:433470/1‑93 (MQ=255)

                                                                          cACTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCgg                >  1:4124551/1‑93 (MQ=255)

                                                                           ccTTCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGgc               <  2:865476/92‑1 (MQ=255)

                                                                              tCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCt            <  1:2287671/93‑1 (MQ=255)

                                                                              tCAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCt            >  1:2940003/1‑93 (MQ=255)

                                                                               cAATGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCtt           >  1:3073540/1‑93 (MQ=255)

                                                                                aaTGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTc          >  2:3484553/1‑93 (MQ=255)

                                                                                 aTGCCGAGGATGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCa         <  1:4576345/93‑1 (MQ=255)

                                                                                 aTGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCa         <  1:735110/93‑1 (MQ=255)

                                                                                 aTGCCGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACG‑CCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCAc        <  2:2617610/93‑1 (MQ=255)

                                                                                     cGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGc     >  2:2210688/1‑93 (MQ=255)

                                                                                     cGAGGAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCACCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGc     >  2:1321826/1‑93 (MQ=255)

                                                                                      gaggagGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGcc    >  1:858515/1‑93 (MQ=255)

                                                                                        ggAGGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCaa  >  2:3989385/1‑93 (MQ=255)

TGACACATCACGCGCAGGTTGCTTAAACCTTCTTCATGCGCAGAAGCCAGGCACGCACCAACGCCCGGTGAATGCACTTCAATGCCGAGGAAGTCCTGCTCAGGGCGATCTTTAGCCATTGCCACCAGCGACGCCCCCATGCCAAAACCAATCTCAAGCGTCACCGGCGCTTCACGGCCAA  >  NC_000913/3100512‑3100692

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