BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	216,444	C→T	S187N (AGC→AAC)	yaeF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	216,444	0	C	T	99.5%	742.5 / NA	207	S187N (AGC→AAC)	yaeF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (92/114); total (93/114)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 211 total aligned reads displayed.

TTGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGCTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCG  >  NC_000913/216354‑216533

    cAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGc                                                                                     <  2:395466/93‑1 (MQ=255)

     aaaGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCa                                                                                    >  1:2680578/1‑93 (MQ=255)

      aaGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCAc                                                                                   >  1:3968104/1‑93 (MQ=255)

       aGGCATCCGTGTCAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACt                                                                                  <  1:2283995/93‑1 (MQ=255)

         gCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGt                                                                                <  1:2081508/93‑1 (MQ=255)

          cATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGtt                                                                               <  2:4266863/93‑1 (MQ=255)

            tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                             >  1:1049724/1‑93 (MQ=255)

            tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                             >  1:2734161/1‑93 (MQ=255)

            tCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTgg                                                                             <  2:3332329/93‑1 (MQ=255)

             ccGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACAGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGc                                                                            <  1:1929588/93‑1 (MQ=255)

                tGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaa                                                                         >  1:805198/1‑93 (MQ=255)

                tGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaa                                                                         >  2:2284235/1‑93 (MQ=255)

                 gACAAATTCCGAACAAAACCACGAGTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaa                                                                        <  2:4105682/93‑1 (MQ=255)

                 gACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaa                                                                        >  2:2149285/1‑93 (MQ=255)

                  aCAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaaa                                                                       >  2:2677703/1‑93 (MQ=255)

                  aCAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGaaaa                                                                       >  2:2921560/1‑93 (MQ=255)

                   cAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAAt                                                                      >  1:3787272/1‑93 (MQ=255)

                   cAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAAt                                                                      <  1:2960405/93‑1 (MQ=255)

                   cAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAAt                                                                      <  1:31724/93‑1 (MQ=255)

                     aaTTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATcc                                                                    >  2:3480191/1‑93 (MQ=255)

                       ttCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTc                                                                  >  2:3732263/1‑93 (MQ=255)

                       ttCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTc                                                                  <  1:3443344/93‑1 (MQ=255)

                         ccGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCgg                                                                <  1:115283/93‑1 (MQ=255)

                          cGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGga                                                               <  1:153371/93‑1 (MQ=255)

                           gAACACAAGCACGTCTTGTTGGCTCCTTCGCCCACACTGCTTATCTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGgag                                                              <  2:3597512/93‑1 (MQ=255)

                           gAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGgag                                                              >  1:1209050/1‑93 (MQ=255)

                            aaCAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGgaga                                                             >  2:3692165/1‑93 (MQ=255)

                             aaaaaaCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGaa                                                            >  2:4488788/3‑93 (MQ=255)

                             aCAAAACCACGCCTTTTTGCCTCCTTCGCCCACACTGCTTAACTGCGCTTTCCCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGaa                                                            <  2:383501/93‑1 (MQ=255)

                             aCAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGaa                                                            <  1:2137040/93‑1 (MQ=255)

                              cAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAAc                                                           <  1:1730186/93‑1 (MQ=255)

                                aaaaCACGACTTTTTGTCTCCTTTGGCCACACTGCTTAAATTCGGCTTTGCCAGGCCGGTGACGCACTGTTGGGCAAAATCCCCCGAGGACgg                                                         >  2:1235436/1‑93 (MQ=255)

                                  aCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGAt                                                       >  2:588626/1‑93 (MQ=255)

                                  aCCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGAt                                                       <  1:4104861/93‑1 (MQ=255)

                                     ccGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCa                                                    <  2:335656/92‑1 (MQ=255)

                                     aCGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCa                                                    >  1:3723082/1‑93 (MQ=255)

                                     aCGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCa                                                    >  1:3850096/1‑93 (MQ=255)

                                      cGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAg                                                   <  1:1056557/93‑1 (MQ=255)

                                       gACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGt                                                  <  1:3635503/93‑1 (MQ=255)

                                        aCTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTg                                                 <  1:3607153/93‑1 (MQ=255)

                                         cTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGa                                                <  2:4069409/93‑1 (MQ=255)

                                          tttttGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAg                                               <  1:2069753/93‑1 (MQ=255)

                                           ttttGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGc                                              >  2:631820/1‑93 (MQ=255)

                                            tttGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGca                                             <  2:3851616/93‑1 (MQ=255)

                                              tGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGcaca                                           <  1:224787/93‑1 (MQ=255)

                                                tctcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATc                                         <  1:2536198/93‑1 (MQ=255)

                                                tctcCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATc                                         <  1:1457608/93‑1 (MQ=255)

                                                   ccTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGACAATCCTCGGAGAACGGATTCAGTGAGCACATCTgg                                      >  1:2305379/1‑93 (MQ=255)

                                                   ccTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTgg                                      <  2:3724988/93‑1 (MQ=255)

                                                    cTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGc                                     >  1:47016/1‑93 (MQ=255)

                                                    cTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGc                                     <  1:3122761/93‑1 (MQ=255)

                                                     ttCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCg                                    >  2:198106/1‑93 (MQ=255)

                                                      tCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGa                                   <  2:1111585/93‑1 (MQ=255)

                                                       cGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAg                                  >  1:531474/1‑93 (MQ=255)

                                                        gCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGt                                 >  2:2267711/1‑93 (MQ=255)

                                                        gCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGt                                 <  2:1250767/93‑1 (MQ=255)

                                                        gCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGt                                 <  1:4234805/93‑1 (MQ=255)

                                                         cccACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTc                                >  2:2886400/1‑93 (MQ=255)

                                                          ccACACTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCa                               >  1:3155/1‑93 (MQ=255)

                                                           cacaCTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAc                              >  2:761386/1‑93 (MQ=255)

                                                            acacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCAcc                             <  1:3367204/93‑1 (MQ=255)

                                                             cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                            >  2:2400343/1‑93 (MQ=255)

                                                             cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                            >  1:4156964/1‑93 (MQ=255)

                                                             cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                            >  1:1855850/1‑93 (MQ=255)

                                                             cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                            >  2:1911911/1‑93 (MQ=255)

                                                             cacTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCa                            >  2:3017679/1‑93 (MQ=255)

                                                              acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                           >  2:3661252/1‑93 (MQ=255)

                                                              acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                           <  2:4048193/93‑1 (MQ=255)

                                                              acTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCAt                           <  2:3062788/93‑1 (MQ=255)

                                                               cTGCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATa                          >  2:4492112/1‑93 (MQ=255)

                                                                 gCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATaaa                        >  2:3223094/1‑93 (MQ=255)

                                                                 gCTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATaaa                        >  1:1970489/1‑93 (MQ=255)

                                                                  cTTAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAg                       <  1:3396402/93‑1 (MQ=255)

                                                                   ttAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAgg                      <  1:894095/93‑1 (MQ=255)

                                                                    tAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAggg                     >  2:3168696/1‑93 (MQ=255)

                                                                    tAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAggg                     <  2:1428696/93‑1 (MQ=255)

                                                                    tAACTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAggg                     <  1:1244694/93‑1 (MQ=255)

                                                                      aCTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGaa                   <  2:2786950/93‑1 (MQ=255)

                                                                      aCTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGaa                   <  2:3816726/93‑1 (MQ=255)

                                                                       cggcgcTTTCGCCAGGCCGTTGACGTACTGTCGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCTCCATAAAGGGAAt                  <  2:3714451/91‑1 (MQ=255)

                                                                       cTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAAt                  <  2:528359/93‑1 (MQ=255)

                                                                       cTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAAt                  >  1:564091/1‑93 (MQ=255)

                                                                       cTGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAAt                  <  1:1687075/93‑1 (MQ=255)

                                                                        tGCGCTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATa                 <  2:455531/93‑1 (MQ=255)

                                                                         gcgcTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATaa                <  2:2858028/93‑1 (MQ=255)

                                                                            cTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAAtt             <  1:1486917/93‑1 (MQ=255)

                                                                            cTTTCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAAtt             <  2:1359503/93‑1 (MQ=255)

                                                                            cTTTCGCCAGGCCGTTGACGCACTATTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAAtt             >  2:2241623/1‑93 (MQ=255)

                                                                             tttCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTc            <  2:4529372/93‑1 (MQ=255)

                                                                               tCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGa          <  2:1093909/93‑1 (MQ=255)

                                                                               tCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGa          <  2:502540/93‑1 (MQ=255)

                                                                               tCGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGa          <  1:4390543/93‑1 (MQ=255)

                                                                                cGCCAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGAc         <  1:304894/93‑1 (MQ=255)

                                                                                  ccAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACaa       >  2:2529865/1‑93 (MQ=255)

                                                                                   cAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAAt      >  1:3256712/1‑93 (MQ=255)

                                                                                   cAGGCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAAt      <  2:2822925/93‑1 (MQ=255)

                                                                                      gCCGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGcc   <  2:3996466/93‑1 (MQ=255)

                                                                                       ccGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCcg  <  1:679827/93‑1 (MQ=255)

                                                                                       ccGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCcg  >  1:4603196/1‑93 (MQ=255)

                                                                                       ccGTTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCcg  <  2:1499162/93‑1 (MQ=255)

TTGGCAAAGGCATCCGTGACAAATTCCGAACAAAACCACGACTTTTTGTCTCCTTCGCCCACACTGCTTAACTGCGCTTTCGCCAGGCCGCTGACGCACTGTTGGCGAAAATCCTCGGAGAACGGATTCAGTGAGCACATCTGGCGAGTCACCATAAAGGGAATAAATTCGACAATGCCG  >  NC_000913/216354‑216533

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