BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,461,874	T→C	P104P (CCT→CCC)	b1399 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,461,874	0	T	C	98.8%	513.9 / ‑3.1	162	P104P (CCT→CCC)	b1399	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (84/76); total (84/78)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 173 total aligned reads displayed.

CATCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGG  >  NC_000913/1461784‑1461963

cATCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGc                                                                                         >  1:3718340/1‑93 (MQ=255)

 aTCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCa                                                                                        <  2:1966643/93‑1 (MQ=255)

 aTCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCa                                                                                        >  1:663711/1‑93 (MQ=255)

  tCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCAt                                                                                       >  2:3149058/1‑93 (MQ=255)

  tCGGGCGACGCAGCTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGCAAGAAAAATTAATCCCGCCGGCCAACCCGCag                                                                                       >  1:3516289/1‑92 (MQ=255)

   cGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATg                                                                                      <  2:3171226/93‑1 (MQ=255)

    gggCGACGCAGTTTCTATAGCCTCAGTGATAAATGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATgg                                                                                     >  1:1856398/1‑93 (MQ=255)

    gggCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATgg                                                                                     >  2:2488455/1‑93 (MQ=255)

     ggCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATggg                                                                                    >  1:3468562/1‑93 (MQ=255)

       cGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGAt                                                                                  <  1:861123/93‑1 (MQ=255)

        gACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATg                                                                                 >  2:4179513/1‑93 (MQ=255)

          cGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGt                                                                               >  2:4582143/1‑93 (MQ=255)

            cAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGGAAAATTTATCTCTCAGAGCAACCCGCATGGGATGGGaa                                                                             >  1:4658539/1‑93 (MQ=255)

            cAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGAAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTaa                                                                             >  2:3682467/1‑93 (MQ=255)

             aGTTTCTATAGCCTCAGTGATAAAGGATTGCGCCTGACGCGACCGGCAGACAGTAAAATTTATCGCTCAGATTCACCAGCATGGAATGGCaag                                                                            >  2:3061085/1‑92 (MQ=255)

              gTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAAt                                                                           >  2:239126/1‑93 (MQ=255)

               tttCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATg                                                                          >  2:1426681/1‑93 (MQ=255)

                ttCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATgg                                                                         >  1:4150933/1‑93 (MQ=255)

                  cTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATgggt                                                                       >  1:1263175/1‑91 (MQ=255)

                   tataGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTc                                                                      >  2:3081654/1‑93 (MQ=255)

                    ataGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTcc                                                                     >  2:1715658/1‑93 (MQ=255)

                     taGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAAGGGCTCCt                                                                    >  2:4032962/1‑93 (MQ=255)

                      aGCCTCAGTGATAAAGGCTTGCGCCTGGCGCGACGCGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTg                                                                   <  1:1261169/93‑1 (MQ=255)

                      aGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTg                                                                   <  1:4680288/93‑1 (MQ=255)

                       gCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGt                                                                  <  2:4675131/93‑1 (MQ=255)

                        ccTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGtt                                                                 >  1:4708721/1‑93 (MQ=255)

                        ccTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGtt                                                                 >  2:2456540/1‑93 (MQ=255)

                            aGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAACTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGctc                                                             <  2:817397/93‑1 (MQ=255)

                             gTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGctct                                                            <  1:1420930/93‑1 (MQ=255)

                             gTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGAAACCCGCATGGGATGGCAAATGGCTCCTGTTGctct                                                            >  1:1657088/1‑93 (MQ=255)

                              tGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCtctc                                                           <  1:3662856/93‑1 (MQ=255)

                                aTAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCgg                                                         <  2:4483586/93‑1 (MQ=255)

                                  aaaGGCTTGCGCGTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGaa                                                       <  2:571512/93‑1 (MQ=255)

                                  aaaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCTCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGaa                                                       <  1:2819249/93‑1 (MQ=255)

                                  aaaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGaa                                                       <  2:1025532/93‑1 (MQ=255)

                                   aaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAg                                                      >  2:3705193/1‑93 (MQ=255)

                                   aaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAg                                                      <  2:1751838/93‑1 (MQ=255)

                                   aaGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAg                                                      <  1:415359/93‑1 (MQ=255)

                                     ggCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGt                                                    >  1:2246685/1‑93 (MQ=255)

                                      gCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGtt                                                   >  1:358836/1‑93 (MQ=255)

                                       cTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGttt                                                  <  1:866723/93‑1 (MQ=255)

                                        ttGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTa                                                 >  1:2207049/1‑93 (MQ=255)

                                         tGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAg                                                <  2:3025585/93‑1 (MQ=255)

                                           cgcCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGAt                                              >  2:777042/1‑93 (MQ=255)

                                             ccTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTCTTGCTCTCGGAAGGTTTAGATaa                                            >  2:3391732/1‑93 (MQ=255)

                                              cTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATaaa                                           <  2:562341/93‑1 (MQ=255)

                                               tGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAt                                          >  2:2082724/1‑93 (MQ=255)

                                               tGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAt                                          <  2:351879/93‑1 (MQ=255)

                                               tGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAAt                                          >  1:2510842/1‑93 (MQ=255)

                                               tGACGCGACCGGTAGAACGCAAAATTTATCGGCCAGAGCAACCCGCATTGGATGGTAAAGTGCTCCTGTTGCTATCGGAAGGTTTAGATAAAt                                          <  2:351880/93‑1 (MQ=255)

                                                gACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATc                                         <  2:1277075/93‑1 (MQ=255)

                                                gACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATc                                         <  1:4326924/93‑1 (MQ=255)

                                                gACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATc                                         <  2:332393/93‑1 (MQ=255)

                                                 aCGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCa                                        <  1:2694259/93‑1 (MQ=255)

                                                 aCGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCa                                        <  2:481606/93‑1 (MQ=255)

                                                  cgcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCaa                                       <  2:2232611/93‑1 (MQ=255)

                                                  cgcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCaa                                       <  2:2087489/93‑1 (MQ=255)

                                                  cgcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCaa                                       <  1:2196295/93‑1 (MQ=255)

                                                   gcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAAc                                      <  1:1231564/93‑1 (MQ=255)

                                                   gcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAAc                                      <  1:2265429/93‑1 (MQ=255)

                                                   gcgACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAAc                                      >  1:179278/1‑93 (MQ=255)

                                                     gACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgc                                    >  2:3278223/1‑93 (MQ=255)

                                                      gcGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgct                                   <  1:2415940/92‑1 (MQ=255)

                                                      aCGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgct                                   >  2:2793290/1‑93 (MQ=255)

                                                       ggggAAGAAAAAAAAATTTACCGCGCAGGGCAACCCGCATGGGATGGTAACTGCCTCGTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctg                                  <  2:309896/92‑1 (MQ=255)

                                                       cGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACgctg                                  <  1:716321/93‑1 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             >  2:1135729/1‑93 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             <  1:2019423/93‑1 (MQ=255)

                                                            aGAAAGTAAAATTTATCGCGCACAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGa                             >  1:2262174/1‑93 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAGGGTTTAGATAAATCAACGCTGGCTGAt                            >  2:1895739/1‑93 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            <  2:421113/93‑1 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            <  1:1090097/93‑1 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            >  1:478677/1‑93 (MQ=255)

                                                             gAAAGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGAt                            >  1:4196118/1‑93 (MQ=255)

                                                               aaGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGt                          <  2:1332011/93‑1 (MQ=255)

                                                               aaGTAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGt                          >  2:2650536/1‑93 (MQ=255)

                                                                  tAAAATTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTcca                       >  1:327502/1‑91 (MQ=255)

                                                                   aaaaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaa                      <  1:3630676/93‑1 (MQ=255)

                                                                     aaTTTATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCaaaaa                    >  1:4656131/1‑93 (MQ=255)

                                                                        ttATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCa                 <  1:3326676/93‑1 (MQ=255)

                                                                         tATCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAg                >  2:399745/1‑93 (MQ=255)

                                                                          aTCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGt               >  2:4205252/1‑93 (MQ=255)

                                                                           tCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGtt              <  1:2744583/93‑1 (MQ=255)

                                                                           tCGCGCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGtt              <  2:3915618/93‑1 (MQ=255)

                                                                            cgcgCAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTg             <  1:1070870/93‑1 (MQ=255)

                                                                              cgcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGAt           <  2:3739431/93‑1 (MQ=255)

                                                                              cgcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGAt           >  2:4040341/1‑93 (MQ=255)

                                                                               gcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATc          >  2:3544978/1‑93 (MQ=255)

                                                                               gcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATc          >  2:3583484/1‑93 (MQ=255)

                                                                               gcAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATc          >  2:1508079/1‑93 (MQ=255)

                                                                                cAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCt         <  1:3190759/93‑1 (MQ=255)

                                                                                cAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCt         <  1:3751785/93‑1 (MQ=255)

                                                                                cAGAGCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCt         <  2:2857848/93‑1 (MQ=255)

                                                                                 agagCAACGCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTg        <  2:1198319/93‑1 (MQ=255)

                                                                                 agagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTg        <  2:4055910/93‑1 (MQ=255)

                                                                                 agagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTg        >  1:388091/1‑93 (MQ=255)

                                                                                  gagCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTgg       >  1:529474/1‑93 (MQ=255)

                                                                                    gCAACCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGGTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCa     >  2:1356388/1‑93 (MQ=255)

                                                                                      aaCCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAg   <  1:3152190/93‑1 (MQ=255)

                                                                                       aCCCGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAgg  >  1:2966029/1‑93 (MQ=255)

CATCGGGCGACGCAGTTTCTATAGCCTCAGTGATAAAGGCTTGCGCCTGACGCGACGGGCAGAAAGTAAAATTTATCGCGCAGAGCAACCTGCATGGGATGGTAAATGGCTCCTGTTGCTCTCGGAAGGTTTAGATAAATCAACGCTGGCTGATGTCAAAAAGCAGTTGATCTGGCAAGG  >  NC_000913/1461784‑1461963

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: `‑`