BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	366,290	C→T	92.0%	D6N (GAT→AAT)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,290	0	C	T	92.0%	3044.3 / 170.2	874	D6N (GAT→AAT)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand): ref base C (34/36); new base T (381/423); total (415/459)
Fisher's exact test for biased strand distribution p-value = 9.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 878 total aligned reads displayed.

ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366192‑366388

aCGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATc                                                                                                    <  2:1131897/99‑1 (MQ=255)

  tgcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCg                                                                                                  >  1:880252‑M2/93‑99 (MQ=255)

     accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                               <  1:226105‑M2/10‑1 (MQ=255)

     accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                               <  1:1078566‑M2/10‑1 (MQ=255)

     accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                               >  1:69820‑M2/90‑99 (MQ=255)

      gCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAAt                                                                                              >  2:747897/1‑99 (MQ=255)

        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  1:397145‑M2/13‑1 (MQ=255)

        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  2:1001524‑M2/13‑1 (MQ=255)

        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  1:856747‑M2/13‑1 (MQ=255)

        ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  2:1317748‑M2/13‑1 (MQ=255)

         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           <  1:558969‑M2/14‑1 (MQ=255)

         gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                           >  2:356695‑M2/86‑99 (MQ=255)

           gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         <  2:260356‑M2/16‑2 (MQ=255)

           gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         >  1:471644‑M2/84‑98 (MQ=255)

           gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         <  2:449268‑M2/16‑2 (MQ=255)

             cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                       <  2:1280753‑M2/18‑4 (MQ=255)

             cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                       <  1:532459‑M2/18‑4 (MQ=255)

             cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                       <  1:184540‑M2/18‑4 (MQ=255)

              agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                      <  2:718493‑M2/19‑5 (MQ=255)

               gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATa                                                                                    <  2:641867‑M2/20‑1 (MQ=255)

               gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATa                                                                                    >  1:1142744‑M2/80‑99 (MQ=255)

               gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATa                                                                                    >  2:1168251‑M2/80‑99 (MQ=255)

                cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAg                                                                                   <  1:1342927‑M2/21‑1 (MQ=255)

                cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAg                                                                                   <  1:1374477‑M2/21‑1 (MQ=255)

                 ccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGc                                                                                  <  2:88199‑M2/22‑1 (MQ=255)

                  cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCt                                                                                 <  2:1230375‑M2/23‑1 (MQ=255)

                   gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTg                                                                                >  1:254812‑M2/76‑99 (MQ=255)

                    gATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTg                                                                                <  2:787985/99‑1 (MQ=255)

                    agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                               >  1:36758‑M2/75‑99 (MQ=255)

                    agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                               <  1:1270705‑M2/25‑1 (MQ=255)

                    agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                               <  1:286941‑M2/25‑1 (MQ=255)

                     gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGtt                                                                              <  2:1343547‑M2/26‑1 (MQ=255)

                      cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGttt                                                                             >  1:1056284‑M2/73‑99 (MQ=255)

                       gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTc                                                                            <  2:811452‑M2/28‑1 (MQ=255)

                        agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTcc                                                                           <  2:569388‑M2/29‑1 (MQ=255)

                          ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtg                                                                         <  2:46524‑M2/31‑1 (MQ=255)

                          ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtg                                                                         >  1:603613‑M2/69‑99 (MQ=255)

                          ctgcAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCt                                                                          >  1:794452/1‑99 (MQ=255)

                            ccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtg                                                                       <  2:163770‑M2/33‑1 (MQ=255)

                            ccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtg                                                                       <  1:93722‑M2/33‑1 (MQ=255)

                             cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtgt                                                                      <  2:394870‑M2/34‑1 (MQ=255)

                              gtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtgtg                                                                     <  1:204544‑M2/35‑1 (MQ=255)

                               tgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGa                                                                    >  2:695149‑M2/64‑99 (MQ=255)

                                ggCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGa                                                                    >  2:1379525/1‑99 (MQ=255)

                                gccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGaa                                                                   <  2:1403361‑M2/37‑1 (MQ=255)

                                 ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGaaa                                                                  <  1:432964‑M2/38‑1 (MQ=255)

                                 ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGaaa                                                                  <  1:57917‑M2/38‑1 (MQ=255)

                                  cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAAt                                                                 <  1:1406088‑M2/39‑1 (MQ=255)

                                  cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAAt                                                                 <  1:547893‑M2/39‑1 (MQ=255)

                                  cGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGaaa                                                                  <  1:1257769/99‑1 (MQ=255)

                                   ggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAAtt                                                                >  1:117889‑M2/60‑99 (MQ=255)

                                      tgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGtt                                                             <  1:731510‑M2/43‑1 (MQ=255)

                                         gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                          <  1:1295565‑M2/46‑1 (MQ=255)

                                         gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                          >  2:230950‑M2/54‑99 (MQ=255)

                                          aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                         <  1:1213214‑M2/47‑1 (MQ=255)

                                          aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                         <  1:50541‑M2/47‑1 (MQ=255)

                                            gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                       >  1:975681‑M2/51‑99 (MQ=255)

                                             aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                      >  1:1015746‑M2/50‑99 (MQ=255)

                                             aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                      >  1:916160‑M2/50‑99 (MQ=255)

                                              aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                     >  2:609143‑M2/49‑99 (MQ=255)

                                               aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                    <  2:375294‑M2/52‑1 (MQ=255)

                                                aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                   <  2:856203‑M2/53‑1 (MQ=255)

                                                 agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                  >  1:71835‑M2/46‑99 (MQ=255)

                                                  gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                 <  1:131797‑M2/55‑1 (MQ=255)

                                                  gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                 >  2:32467‑M2/45‑99 (MQ=255)

                                                    gaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                               <  2:1259206‑M2/57‑1 (MQ=255)

                                                    gaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                               <  2:283166‑M2/57‑1 (MQ=255)

                                                     aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                              <  2:946016‑M2/58‑1 (MQ=255)

                                                     aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTGTCCTGTGTGAAATTGTTATCCGCTCTCAATTc                                              >  2:420117‑M2/42‑99 (MQ=255)

                                                      atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                             <  1:1360102‑M2/59‑1 (MQ=255)

                                                        gatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcac                                           >  1:819867‑M2/39‑99 (MQ=255)

                                                        gatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcac                                           >  2:1465056‑M2/39‑99 (MQ=255)

                                                         atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                          <  2:58876‑M2/62‑1 (MQ=255)

                                                         atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                          >  1:1030862‑M2/38‑99 (MQ=255)

                                                         atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                          <  2:1201416‑M2/62‑1 (MQ=255)

                                                            tagccgtcaagttgtcataattggtaccgcgctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                       <  2:979037‑M2/65‑1 (MQ=255)

                                                              gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                     >  1:672940‑M2/33‑99 (MQ=255)

                                                                  tcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:1229185‑M2/29‑99 (MQ=255)

                                                                  tcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 <  2:1249150‑M2/71‑1 (MQ=255)

                                                                   tCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:1122175/1‑99 (MQ=255)

                                                                   caagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGa                                <  1:1436543‑M2/72‑1 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              >  2:191972‑M2/26‑99 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              >  2:1116689‑M2/26‑99 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              <  2:1083006‑M2/74‑1 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              >  2:51384‑M2/26‑99 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              <  1:762010‑M2/74‑1 (MQ=255)

                                                                     agttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              >  1:1454523‑M2/26‑99 (MQ=255)

                                                                       ttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                            >  1:1106049‑M2/24‑99 (MQ=255)

                                                                         gtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGa                          >  1:160570‑M2/22‑99 (MQ=255)

                                                                              aattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                     >  2:652622‑M2/17‑99 (MQ=255)

                                                                               attggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    >  2:1032343‑M2/16‑99 (MQ=255)

                                                                               attggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    >  2:825576‑M2/16‑99 (MQ=255)

                                                                               attggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    >  2:924784‑M2/16‑99 (MQ=255)

                                                                                   gtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                >  2:1052923‑M2/12‑99 (MQ=255)

                                                                                   gtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                <  2:549650‑M2/88‑1 (MQ=255)

                                                                                          gctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc         <  2:483281‑M2/95‑1 (MQ=255)

                                                                                          gctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc         >  2:1441414‑M2/5‑99 (MQ=255)

                                                                                           ctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc        <  1:321921‑M2/96‑1 (MQ=255)

                                                                                             cGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg      <  1:629644‑M2/98‑1 (MQ=255)

                                                                                             cGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg      <  2:250971‑M2/98‑1 (MQ=255)

ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366192‑366388

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