BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,291,253	T→C	R71R (CGT→CGC)	yraM →	putative glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,291,253	0	T	C	81.8%	603.4 / 137.0	286	R71R (CGT→CGC)	yraM	putative glycosylase
Reads supporting (aligned to +/- strand): ref base T (26/26); new base C (113/121); total (139/147)
Fisher's exact test for biased strand distribution p-value = 8.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 287 total aligned reads displayed.

CACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAA  >  NC_000913/3291163‑3291344

 aCGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCa                                                                                          <  1:5981021/93‑1 (MQ=255)

 aCGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCa                                                                                          <  1:4521269/93‑1 (MQ=255)

  cGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCAc                                                                                         >  1:285858/1‑93 (MQ=255)

    gcgcAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTg                                                                                       >  2:2487578/1‑93 (MQ=255)

     cgcAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTgg                                                                                      <  2:5040208/93‑1 (MQ=255)

      gcAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGt                                                                                     <  2:3097658/93‑1 (MQ=255)

       cAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTg                                                                                    <  1:464602/93‑1 (MQ=255)

          gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGaaa                                                                                 <  1:5065524/93‑1 (MQ=255)

          gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGaaa                                                                                 >  1:3591641/1‑93 (MQ=255)

          gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGaaa                                                                                 <  2:740003/93‑1 (MQ=255)

          gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGaaa                                                                                 <  1:2974122/93‑1 (MQ=255)

          gCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGaaa                                                                                 <  2:6968624/93‑1 (MQ=255)

           cTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAaag                                                                                >  2:5858545/1‑93 (MQ=255)

            tGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAaaga                                                                               >  1:2206079/1‑93 (MQ=255)

             gATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAaagaa                                                                              >  2:5811907/1‑93 (MQ=255)

             gATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAaagaa                                                                              <  1:2010510/93‑1 (MQ=255)

               ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAgg                                                                            <  1:1549546/93‑1 (MQ=255)

               ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAgg                                                                            >  2:1459839/1‑93 (MQ=255)

               ttCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAgg                                                                            >  2:1851023/1‑93 (MQ=255)

                 cTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTa                                                                          <  2:5264240/93‑1 (MQ=255)

                  tGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTaa                                                                         <  2:210582/93‑1 (MQ=255)

                  tGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTaa                                                                         >  2:5381725/1‑93 (MQ=255)

                  tGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTaa                                                                         >  2:127826/1‑93 (MQ=255)

                    ccTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTaaaa                                                                       <  1:6776004/93‑1 (MQ=255)

                     cTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAAc                                                                      <  1:5133633/93‑1 (MQ=255)

                     cTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAAc                                                                      <  2:705946/93‑1 (MQ=255)

                       tttATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTGCTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCg                                                                    <  2:425068/93‑1 (MQ=255)

                       tttATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCg                                                                    >  2:500438/1‑93 (MQ=255)

                       tttATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCg                                                                    <  2:5398295/93‑1 (MQ=255)

                          aTCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGc                                                                 <  1:673324/93‑1 (MQ=255)

                           tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCa                                                                >  2:164527/1‑93 (MQ=255)

                           tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCa                                                                >  2:5724248/1‑93 (MQ=255)

                           tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCa                                                                <  2:6480325/93‑1 (MQ=255)

                           tCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCa                                                                <  1:1819789/93‑1 (MQ=255)

                             ttCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAgg                                                              >  1:3462153/1‑93 (MQ=255)

                              tCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGc                                                             <  1:6720971/93‑1 (MQ=255)

                               cagcagATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCg                                                            >  2:5503141/1‑93 (MQ=255)

                                agcagATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCgg                                                           >  1:7190975/1‑93 (MQ=255)

                                 gcagATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGt                                                          >  1:5425995/1‑93 (MQ=255)

                                   agATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTg                                                        <  2:3587378/93‑1 (MQ=255)

                                    gATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGa                                                       >  2:2848683/1‑93 (MQ=255)

                                    gATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGa                                                       >  1:477823/1‑93 (MQ=255)

                                     aTGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAg                                                      <  2:732808/93‑1 (MQ=255)

                                       gcagcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgtt                                                    >  1:568635/1‑93 (MQ=255)

                                        cagcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttg                                                   <  1:4390910/93‑1 (MQ=255)

                                         agcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttgt                                                  >  2:2537635/1‑93 (MQ=255)

                                         agcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttgt                                                  >  2:4437170/1‑93 (MQ=255)

                                          gcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttgtt                                                 >  2:2185784/1‑93 (MQ=255)

                                          gcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttgtt                                                 >  1:2231810/1‑93 (MQ=255)

                                          gcaAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAgttgtt                                                 <  2:3694403/93‑1 (MQ=255)

                                           caAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGttt                                                >  1:3907444/1‑93 (MQ=255)

                                           caAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGttt                                                <  1:4716491/93‑1 (MQ=255)

                                            aaaGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTa                                               >  1:256478/1‑93 (MQ=255)

                                            aaaGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTa                                               >  2:2625461/1‑93 (MQ=255)

                                             aaGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTaa                                              <  1:2673025/93‑1 (MQ=255)

                                             aaGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTaa                                              >  1:6940380/1‑93 (MQ=255)

                                              aGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAAc                                             <  2:340804/93‑1 (MQ=255)

                                              aGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAAc                                             <  1:3275864/93‑1 (MQ=255)

                                               gCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAAcc                                            >  2:1929560/1‑93 (MQ=255)

                                               gCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAAcc                                            >  1:7427870/1‑93 (MQ=255)

                                               gCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAAcc                                            >  2:6874960/1‑93 (MQ=255)

                                                ctctGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCa                                           >  2:2740750/1‑93 (MQ=255)

                                                    gatgatACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTa                                       >  2:2312491/1‑93 (MQ=255)

                                                    gatgatACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTa                                       >  1:416402/1‑93 (MQ=255)

                                                     atgatACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTAc                                      <  2:3342117/93‑1 (MQ=255)

                                                     atgatACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTAc                                      <  2:6682593/93‑1 (MQ=255)

                                                     atgatACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTAc                                      <  2:1049761/93‑1 (MQ=255)

                                                        atACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGc                                   >  2:6352533/1‑93 (MQ=255)

                                                         tACCAGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCa                                  <  1:5894885/93‑1 (MQ=255)

                                                          aCCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCaa                                 <  2:7256449/93‑1 (MQ=255)

                                                             aGGATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGaa                              >  1:1512611/1‑93 (MQ=255)

                                                              ggATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAAc                             >  1:101475/1‑93 (MQ=255)

                                                               gATCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACt                            >  2:4750701/1‑93 (MQ=255)

                                                                aTCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTg                           >  1:7075348/1‑93 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          <  2:1715796/93‑1 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          <  1:994872/93‑1 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          >  1:894062/1‑93 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          <  2:5782729/93‑1 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          >  1:6549740/1‑93 (MQ=255)

                                                                 tCAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGa                          <  1:6080942/93‑1 (MQ=255)

                                                                  cAACTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGaa                         <  1:1673385/93‑1 (MQ=255)

                                                                   aaCTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAAc                        >  1:4294204/1‑93 (MQ=255)

                                                                     cTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGa                      <  2:3952875/93‑1 (MQ=255)

                                                                     cTGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCACAAGAACTGAACGa                      <  2:7488312/93‑1 (MQ=255)

                                                                      tGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGAt                     <  2:4508218/93‑1 (MQ=255)

                                                                      tGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGAt                     >  2:3131907/1‑93 (MQ=255)

                                                                      tGGCAATTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGAt                     >  1:403482/1‑93 (MQ=255)

                                                                          aaTTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTc                 >  2:3490370/1‑93 (MQ=255)

                                                                           aTTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCa                <  1:2371422/93‑1 (MQ=255)

                                                                           aTTACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCa                >  2:847801/1‑93 (MQ=255)

                                                                             tACTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGc              <  2:3236439/93‑1 (MQ=255)

                                                                               cTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGt            <  1:2354897/93‑1 (MQ=255)

                                                                               cTCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGt            <  1:1423483/93‑1 (MQ=255)

                                                                                tCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTc           >  2:6625289/1‑93 (MQ=255)

                                                                                tCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTc           >  1:339028/1‑93 (MQ=255)

                                                                                tCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTc           <  2:1488015/93‑1 (MQ=255)

                                                                                tCGCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGATGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTc           >  2:6263981/1‑93 (MQ=255)

                                                                                  gCCATTCGCGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTcgc         <  1:4876897/93‑1 (MQ=255)

                                                                                         gTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGaaaa  <  2:854614/93‑1 (MQ=255)

                                                                                         gTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGaaaa  >  1:1311588/1‑93 (MQ=255)

CACGGCGCAGGCTGATTCTGCCTTTTATCTTCAGCAGATGCAGCAAAGCTCTGATGATACCAGGATCAACTGGCAATTACTCGCCATTCGTGCACTGGTGAAAGAAGGTAAAACCGGGCAGGCGGTTGAGTTGTTTAACCAACTACCGCAAGAACTGAACGATGCTCAGCGTCGCGAGAAAA  >  NC_000913/3291163‑3291344

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