BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	139,187	T→C	58.6%	N680S (AAT→AGT)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	139,187	0	T	C	58.6%	114.8 / 210.2	186	N680S (AAT→AGT)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (39/38); new base C (45/64); total (84/102)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 188 total aligned reads displayed.

CCAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAC  >  NC_000913/139089‑139282

ccAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCAt                                                                                                 <  2:1220497/99‑1 (MQ=255)

 cAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTa                                                                                               <  2:368470/100‑1 (MQ=255)

  aGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAg                                                                                              <  2:562502/100‑1 (MQ=255)

   gCATCGGCGTACCCATATTGTACGGCCCCGGAACCGGCATCGCGAACCGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTa                                                                                               <  2:3387719/98‑1 (MQ=255)

   gCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGt                                                                                             >  1:2349810/1‑100 (MQ=255)

        ggCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCa                                                                                        >  1:1052676/1‑100 (MQ=255)

        ggCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCa                                                                                        <  2:1959683/100‑1 (MQ=255)

        ggCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCa                                                                                        <  2:3928371/100‑1 (MQ=255)

        ggCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCa                                                                                        >  2:537404/1‑100 (MQ=255)

          cATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGa                                                                                      <  1:2159379/100‑1 (MQ=255)

            tACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATc                                                                                    <  1:1958927/100‑1 (MQ=255)

             aCCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATcc                                                                                   <  1:762809/100‑1 (MQ=255)

              cccATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCa                                                                                  >  2:2299600/1‑100 (MQ=255)

              cccATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCa                                                                                  <  1:221731/100‑1 (MQ=255)

                cATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAgc                                                                                >  2:716716/1‑100 (MQ=255)

                    ttGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCg                                                                            >  1:82066/1‑100 (MQ=255)

                     tGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCga                                                                           <  1:3337654/100‑1 (MQ=255)

                      gAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCgag                                                                          >  1:162417/1‑100 (MQ=255)

                       aaCGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCgaga                                                                         <  2:1232433/100‑1 (MQ=255)

                        aCGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGat                                                                        >  2:149407/1‑100 (MQ=255)

                         cGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGata                                                                       <  2:470429/100‑1 (MQ=255)

                         cGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGata                                                                       <  2:1253404/100‑1 (MQ=255)

                         cGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGata                                                                       >  1:1523271/1‑100 (MQ=255)

                           gCACCGGTACAGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTGGTTTTCAGGTCCAGCGCCGAGAtata                                                                     >  2:1636821/1‑100 (MQ=255)

                           gCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGAtata                                                                     <  2:3570039/100‑1 (MQ=255)

                           gCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGAtata                                                                     <  2:1384515/100‑1 (MQ=255)

                           gCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGAtata                                                                     <  2:386428/100‑1 (MQ=255)

                           cttccgatctcGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGAtata                                                                     <  2:1694691/90‑1 (MQ=255)

                            cACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGAtata                                                                     <  2:344649/99‑1 (MQ=255)

                               cGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAAccc                                                                 <  2:1409094/100‑1 (MQ=255)

                               cGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAAccc                                                                 <  2:162416/100‑1 (MQ=255)

                                ggAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAAcccc                                                                >  1:1397086/1‑100 (MQ=255)

                                   aCCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGCTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATg                                                             <  2:2881666/100‑1 (MQ=255)

                                    ccGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgc                                                            >  1:3735057/1‑100 (MQ=255)

                                    ccGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgc                                                            <  2:3877438/100‑1 (MQ=255)

                                     cGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATaagatcggaag                                                           >  1:1694703/1‑91 (MQ=255)

                                     cGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgct                                                           >  1:3338645/1‑100 (MQ=255)

                                     cGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgct                                                           <  1:1248687/100‑1 (MQ=255)

                                      ggCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctg                                                          >  1:3523617/1‑100 (MQ=255)

                                      ggCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctg                                                          >  2:2133706/1‑100 (MQ=255)

                                       gCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctgg                                                         <  1:2914744/100‑1 (MQ=255)

                                       gCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctgg                                                         <  1:1828239/100‑1 (MQ=255)

                                        cATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctggc                                                        >  1:391193/1‑100 (MQ=255)

                                          tCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctggctg                                                      >  2:790898/1‑100 (MQ=255)

                                          tCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctggctg                                                      <  2:3036634/100‑1 (MQ=255)

                                          tCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATgctggctg                                                      <  1:1152269/100‑1 (MQ=255)

                                            gggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGtt                                                    >  1:3053142/1‑100 (MQ=255)

                                             ggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGttt                                                   <  1:1115377/100‑1 (MQ=255)

                                             ggAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGttt                                                   >  2:1181243/1‑100 (MQ=255)

                                              gAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTa                                                  >  1:3835069/1‑100 (MQ=255)

                                              gAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTa                                                  >  1:3518816/1‑100 (MQ=255)

                                               aaCGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTAc                                                 <  1:1939832/100‑1 (MQ=255)

                                               aaCGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTAc                                                 <  2:2941945/100‑1 (MQ=255)

                                                aCGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACa                                                <  2:74319/100‑1 (MQ=255)

                                                 cGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACAt                                               <  1:597940/100‑1 (MQ=255)

                                                  ggCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATg                                              <  2:3708333/100‑1 (MQ=255)

                                                    cATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGc                                            <  2:686806/100‑1 (MQ=255)

                                                    cATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGc                                            <  1:1932548/100‑1 (MQ=255)

                                                      tACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAg                                          <  2:3599162/100‑1 (MQ=255)

                                                      tACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCACCGCCGAGATATCACCCCATGCTGGCTGTTTACAt                                               >  2:2041040/1‑95 (MQ=255)

                                                       aCTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGa                                         >  1:175614/1‑100 (MQ=255)

                                                        cTGTCCTGCGGCGTACCAATCCGTTTCTTCCACACCACTTCATTAGTCTTCAGATCCAGCGCCGCGATATAACCCCATGCTGGCTGTTTACATGTCAGAc                                        >  2:940949/1‑100 (MQ=255)

                                                         tGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGAcc                                       >  2:2104401/1‑100 (MQ=255)

                                                          gTCCTGCGGCGTACCAATACGTTTATTCCACACCACTTCATTAGTTTTCAGCTCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGAcc                                       >  2:3542737/1‑99 (MQ=255)

                                                          gTCCTGCAGCGAACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCa                                      >  2:2159070/1‑100 (MQ=255)

                                                            ccTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCaaa                                    >  1:1294930/1‑100 (MQ=255)

                                                             cTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAAt                                   <  1:1259152/100‑1 (MQ=255)

                                                             cTGCGGCGTACCAATACGTTTCTTCCACACCACTTCACTAGTTCTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCGGTGTACATGGCAGACCaaa                                    >  2:622866/1‑99 (MQ=255)

                                                              tGCGGTGTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATg                                  <  1:3809283/100‑1 (MQ=255)

                                                              tGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATg                                  >  1:3413319/1‑100 (MQ=255)

                                                               gcggcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATgg                                 <  1:3774881/100‑1 (MQ=255)

                                                                cggcgTACCAATACGTTTCTTCAACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGt                                >  2:2394131/1‑100 (MQ=255)

                                                                 ggcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTg                               <  1:1008730/100‑1 (MQ=255)

                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  2:3733505/100‑1 (MQ=255)

                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              >  1:4061634/1‑100 (MQ=255)

                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  2:1313503/100‑1 (MQ=255)

                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  1:3119006/100‑1 (MQ=255)

                                                                  gcgTACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTga                              <  1:2513361/100‑1 (MQ=255)

                                                                     tACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTgcga                            >  2:2084713/1‑99 (MQ=255)

                                                                     tACCAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGagag                           <  1:3242510/100‑1 (MQ=255)

                                                                       ccAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGa                         <  2:3410417/100‑1 (MQ=255)

                                                                        cAATACGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGaa                        >  1:3146297/1‑100 (MQ=255)

                                                                           tACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACgg                     <  2:303976/100‑1 (MQ=255)

                                                                            aCGTTTCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACggg                    <  1:2977127/100‑1 (MQ=255)

                                                                               tttCTTCCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTg                 <  1:684456/100‑1 (MQ=255)

                                                                                ttcttcCACACCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGa                >  1:3270662/1‑100 (MQ=255)

                                                                                  cttcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGc              <  1:2032057/100‑1 (MQ=255)

                                                                                    tcCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGt            <  1:154097/100‑1 (MQ=255)

                                                                                     ccACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTg           <  1:3396542/100‑1 (MQ=255)

                                                                                       acacCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGac         <  2:1941933/100‑1 (MQ=255)

                                                                                       acacCACTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCCGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGac         >  2:334965/1‑100 (MQ=255)

                                                                                       acacCACTTCACTAGTGTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGac         <  2:2636378/100‑1 (MQ=255)

                                                                                          ccacTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACAcc      >  1:3299167/1‑100 (MQ=255)

                                                                                           cacTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCa     >  1:363259/1‑100 (MQ=255)

                                                                                            acTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCAt    >  1:2075070/1‑100 (MQ=255)

                                                                                             cTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATa   >  2:2265219/1‑100 (MQ=255)

                                                                                             cTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGCTGAGAGGAACGGGTTGAGCGTGACACCATa   >  2:1774900/1‑100 (MQ=255)

                                                                                             cTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATa   >  1:1800633/1‑100 (MQ=255)

                                                                                             cTTCACTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATa   >  2:89720/1‑100 (MQ=255)

                                                                                              ttCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGGGAGGAACGGGTTGAGCGTGACACCATAc  >  1:2591001/1‑100 (MQ=255)

CCAGCATCGGCATACCCATATTGAACGGCACCGGAACCGGCATCGGGAACGGCATACTGTCCTGCGGCGTACCAATACGTTTCTTCCACACCACTTCATTAGTTTTCAGATCCAGCGCCGAGATATAACCCCATGCTGGCTGTTTACATGGCAGACCAAATGGTGAGAGGAACGGGTTGAGCGTGACACCATAC  >  NC_000913/139089‑139282

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