BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,934,800	T→C	L42P (CTG→CCG)	yebK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,934,800	0	T	C	99.6%	963.4 / NA	281	L42P (CTG→CCG)	yebK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (130/150); total (130/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

AACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCAC  >  NC_000913/1934710‑1934889

aaCATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCgg                                                                                         <  2:4649602/93‑1 (MQ=255)

aaCATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCgg                                                                                         <  1:1616089/93‑1 (MQ=255)

aaCATTTGAGCAAATCAGAGCGCAAAGTTGACGAGGTCATTCTGGCTTCGCCCGATAACTCGATCCATTCGAGTATTGCTGCTATGGCACCgg                                                                                         >  2:5806280/1‑93 (MQ=255)

 aCATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGa                                                                                        >  1:2561838/1‑93 (MQ=255)

  cATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGaa                                                                                       >  2:561910/1‑93 (MQ=255)

  cATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGaa                                                                                       <  2:7706112/93‑1 (MQ=255)

        aGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAAt                                                                                 <  2:5431187/93‑1 (MQ=255)

         gCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATg                                                                                >  2:8542681/1‑93 (MQ=255)

          cAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGt                                                                               <  1:2679354/93‑1 (MQ=255)

          cAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGt                                                                               <  1:556972/93‑1 (MQ=255)

           aaaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGACTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGtt                                                                              <  1:4939353/93‑1 (MQ=255)

            aaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTa                                                                             <  1:225664/93‑1 (MQ=255)

            aaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTa                                                                             <  2:1947152/93‑1 (MQ=255)

            aaTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTa                                                                             <  2:2031692/93‑1 (MQ=255)

             aTCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAg                                                                            <  1:2336544/93‑1 (MQ=255)

              tCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGATAgg                                                                           >  1:7411114/1‑92 (MQ=255)

               cAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCg                                                                          <  2:6975975/93‑1 (MQ=255)

                 aagCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGaa                                                                        <  2:1813823/92‑1 (MQ=255)

                     gcAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGa                                                                    >  2:2352085/1‑93 (MQ=255)

                     gcAAAGTGGGCGAGGTCATTCTGACTACGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGATGCCAATGTTAGCGAACCGa                                                                    <  1:8536001/93‑1 (MQ=255)

                       aaaGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACg                                                                  >  2:2258511/1‑93 (MQ=255)

                       aaaGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACg                                                                  >  1:509405/1‑93 (MQ=255)

                        aaGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGACCCGACgg                                                                 >  1:7189135/1‑93 (MQ=255)

                        aaGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACgg                                                                 >  1:5795464/1‑93 (MQ=255)

                         aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGt                                                                <  1:3386872/93‑1 (MQ=255)

                         aGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGt                                                                <  2:471009/93‑1 (MQ=255)

                          gTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTg                                                               >  2:7307951/1‑93 (MQ=255)

                          gTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTg                                                               >  2:5340475/1‑93 (MQ=255)

                          gTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTg                                                               >  2:2655426/1‑93 (MQ=255)

                           ttGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGa                                                              >  2:8320629/1‑93 (MQ=255)

                              ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATc                                                           >  1:3018261/1‑93 (MQ=255)

                              ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATc                                                           >  1:8073759/1‑93 (MQ=255)

                              ccGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATc                                                           <  1:4237225/93‑1 (MQ=255)

                               cGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCg                                                          >  2:7080169/1‑93 (MQ=255)

                               cGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCg                                                          <  2:4410308/93‑1 (MQ=255)

                                gAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGt                                                         <  1:3979743/93‑1 (MQ=255)

                                 aGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGtt                                                        >  1:1350009/1‑93 (MQ=255)

                                  ggTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGACCCGACGGTGAATCGttt                                                       <  1:76734/93‑1 (MQ=255)

                                   gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTc                                                      >  1:1420233/1‑93 (MQ=255)

                                   gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTc                                                      <  2:7925221/93‑1 (MQ=255)

                                   gTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTc                                                      <  1:6181881/93‑1 (MQ=255)

                                      aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGt                                                   <  2:1622155/93‑1 (MQ=255)

                                      aTTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGt                                                   <  1:8494265/93‑1 (MQ=255)

                                       ttCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTc                                                  >  1:5939839/1‑93 (MQ=255)

                                        tCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCg                                                 >  1:2128964/1‑93 (MQ=255)

                                         cTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCCCCGGAAGTCAATGCTAGCGAACCGACGGTGAATCGTTTCTGTTgc                                                >  1:6808639/1‑93 (MQ=255)

                                         cTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgc                                                <  2:3523005/93‑1 (MQ=255)

                                         cTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgc                                                >  1:503664/1‑93 (MQ=255)

                                          tGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgca                                               >  1:8337829/1‑93 (MQ=255)

                                          tGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgca                                               >  1:8655043/1‑93 (MQ=255)

                                           ggCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcag                                              >  2:6101134/1‑93 (MQ=255)

                                           ggCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcag                                              <  1:27340/93‑1 (MQ=255)

                                             cTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCgcagca                                            <  2:1272365/93‑1 (MQ=255)

                                              ttCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCAt                                           <  1:3815089/93‑1 (MQ=255)

                                               tCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATg                                          >  1:1960497/1‑93 (MQ=255)

                                                cGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATgg                                         >  2:6702202/1‑93 (MQ=255)

                                                 gCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGa                                        <  1:4359611/93‑1 (MQ=255)

                                                     gATAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcg                                    <  1:3346168/93‑1 (MQ=255)

                                                      aTAACGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgc                                   <  1:5806280/93‑1 (MQ=255)

                                                        aaCGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgc                                 <  2:5454714/93‑1 (MQ=255)

                                                        aaCGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgc                                 <  1:7920763/93‑1 (MQ=255)

                                                         aCGCGATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACAcgcgcg                                <  1:4527087/93‑1 (MQ=255)

                                                           gcgATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGt                              <  1:4701313/93‑1 (MQ=255)

                                                            cgATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtt                             <  1:6105330/93‑1 (MQ=255)

                                                             gATCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGttt                            >  1:6234656/1‑93 (MQ=255)

                                                              aTCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtttt                           >  1:5335409/1‑93 (MQ=255)

                                                              aTCCATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGtttt                           >  2:3155381/1‑93 (MQ=255)

                                                                ccATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                         <  2:6563265/93‑1 (MQ=255)

                                                                ccATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTcc                         <  2:4849082/93‑1 (MQ=255)

                                                                ccATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGGTTTcc                         >  2:880327/1‑93 (MQ=255)

                                                                 cATTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCt                        <  1:4454016/93‑1 (MQ=255)

                                                                  aTTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                       <  2:8448784/93‑1 (MQ=255)

                                                                  aTTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                       >  1:2892109/1‑93 (MQ=255)

                                                                  aTTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                       >  2:647838/1‑93 (MQ=255)

                                                                  aTTCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTg                       >  1:7136551/1‑93 (MQ=255)

                                                                   ttCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACTCGCGCGGTTTTCCTGa                      >  1:5510856/1‑93 (MQ=255)

                                                                   ttCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                      >  2:8743449/1‑93 (MQ=255)

                                                                   ttCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGa                      >  1:633680/1‑93 (MQ=255)

                                                                    tCGAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAt                     <  2:8152079/93‑1 (MQ=255)

                                                                      gAGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAttt                   <  2:2849363/93‑1 (MQ=255)

                                                                       aGTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGAtttt                  <  1:6651789/93‑1 (MQ=255)

                                                                        gTATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTa                 >  1:647934/1‑93 (MQ=255)

                                                                         tATTGCTGCTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTaa                >  2:7438802/1‑93 (MQ=255)

                                                                            tgctgcTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACt             >  1:7136047/1‑93 (MQ=255)

                                                                            tgctgcTATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACt             <  2:3948623/93‑1 (MQ=255)

                                                                              ctgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTc           <  1:3726344/93‑1 (MQ=255)

                                                                               tgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCa          <  2:3326123/93‑1 (MQ=255)

                                                                               tgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCa          <  2:820370/93‑1 (MQ=255)

                                                                               tgctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCa          >  1:6381443/1‑93 (MQ=255)

                                                                                gctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCAt         <  1:4686468/93‑1 (MQ=255)

                                                                                 ctATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATc        <  2:4145244/93‑1 (MQ=255)

                                                                                  tATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCt       <  2:1025658/93‑1 (MQ=255)

                                                                                  tATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCt       <  1:3442949/93‑1 (MQ=255)

                                                                                  tATGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCt       >  2:5781022/1‑93 (MQ=255)

                                                                                   aTGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTg      >  1:1068238/1‑93 (MQ=255)

                                                                                   aTGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTg      <  1:7235228/93‑1 (MQ=255)

                                                                                   aTGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTg      <  1:1215311/93‑1 (MQ=255)

                                                                                   aTGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTg      >  1:6905323/1‑93 (MQ=255)

                                                                                    tGGCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTgg     >  1:7442182/1‑93 (MQ=255)

                                                                                     ggCACCGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGc    <  1:1194472/93‑1 (MQ=255)

AACATTTGAGCAAATCAGAGCGCAAAGTTGCCGAGGTCATTCTGGCTTCGCCCGATAACGCGATCCATTCGAGTATTGCTGCTATGGCACTGGAAGCCAATGTTAGCGAACCGACGGTGAATCGTTTCTGTCGCAGCATGGACACGCGCGGTTTTCCTGATTTTAAACTTCATCTGGCAC  >  NC_000913/1934710‑1934889

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