BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	366,301	Δ1 bp	coding (5/3075 nt)	lacZ ←	beta‑galactosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,300	0	G	.	99.7%	95148.1 / 287.9	24907	coding (6/3075 nt)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand): ref base G (31/45); new base . (12623/12196); total (12658/12249)
Fisher's exact test for biased strand distribution p-value = 8.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 26255 total aligned reads displayed.

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398

ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                                    <  1:1633085‑M2/15‑1 (MQ=255)

 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   >  2:1849210‑M2/84‑98 (MQ=255)

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

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

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

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

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

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

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

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

            cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGttt                                                                                       <  2:1800862‑M2/27‑1 (MQ=255)

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

             gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTc                                                                                      >  1:757535‑M2/72‑99 (MQ=255)

              agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTcc                                                                                     >  2:1749446‑M2/71‑99 (MQ=255)

              agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTcc                                                                                     >  2:1386887‑M2/71‑99 (MQ=255)

               gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCt                                                                                    >  1:314378‑M2/70‑99 (MQ=255)

               gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCt                                                                                    <  1:2179430‑M2/30‑1 (MQ=255)

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

                 tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgt                                                                                  <  2:1030297‑M2/32‑1 (MQ=255)

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

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

                       ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGaaa                                                                            >  1:62263‑M2/62‑99 (MQ=255)

                        cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAAt                                                                           >  2:2442643‑M2/61‑99 (MQ=255)

                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  1:2449815‑M2/59‑99 (MQ=255)

                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  1:113896‑M2/59‑99 (MQ=255)

                           ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        >  1:486154‑M2/58‑99 (MQ=255)

                             gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                                      <  2:1475148‑M2/44‑1 (MQ=255)

                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     <  1:186136‑M2/45‑1 (MQ=255)

                                 agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                                  >  1:2137842‑M2/52‑99 (MQ=255)

                                  gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 <  2:2487172‑M2/49‑1 (MQ=255)

                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                >  2:819548‑M2/50‑99 (MQ=255)

                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                <  2:198146‑M2/50‑1 (MQ=255)

                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                >  2:490534‑M2/50‑99 (MQ=255)

                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                >  2:729419‑M2/50‑99 (MQ=255)

                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                <  2:636299‑M2/50‑1 (MQ=255)

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

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

                                        gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                           >  1:235949‑M2/45‑99 (MQ=255)

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

                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:575585‑M2/44‑99 (MQ=255)

                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:1904773‑M2/44‑99 (MQ=255)

                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:80262‑M2/44‑99 (MQ=255)

                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:1409976‑M2/44‑99 (MQ=255)

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

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

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

                                           aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        >  2:1526688‑M2/42‑99 (MQ=255)

                                            atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                                       >  2:1893145‑M2/41‑99 (MQ=255)

                                            atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                                       >  2:345781‑M2/41‑99 (MQ=255)

                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      >  1:444434‑M2/40‑99 (MQ=255)

                                              gatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcac                                                     <  1:356199‑M2/61‑1 (MQ=255)

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

                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    >  2:1070565‑M2/38‑99 (MQ=255)

                                               atgaagccgtcagcttgtcataattggtcccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTCTTCGCTCACAATTCcaca                                                    >  2:956934‑M2/38‑99 (MQ=255)

                                                tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   >  2:44229‑M2/37‑99 (MQ=255)

                                                   agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                                <  1:1590271‑M2/66‑1 (MQ=255)

                                                    gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:2041754‑M2/67‑1 (MQ=255)

                                                    gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:84524‑M2/67‑1 (MQ=255)

                                                    gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:26696‑M2/67‑1 (MQ=255)

                                                     ccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACAt                                              <  2:2268421‑M2/68‑1 (MQ=255)

                                                     ccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACAt                                              >  2:589210‑M2/32‑99 (MQ=255)

                                                      cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             >  1:1000052‑M2/31‑99 (MQ=255)

                                                        tcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  2:2398512‑M2/29‑99 (MQ=255)

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

                                                                tcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                                   <  2:865708‑M2/79‑1 (MQ=255)

                                                                 cataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                                  <  1:1682877‑M2/80‑1 (MQ=255)

                                                                   taattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                <  1:2493677‑M2/82‑1 (MQ=255)

                                                                   taattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                >  1:2222887‑M2/18‑99 (MQ=255)

                                                                   taattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                >  2:911835‑M2/18‑99 (MQ=255)

                                                                    aattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                               >  1:1094980‑M2/17‑99 (MQ=255)

                                                                    aattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                               <  1:1322111‑M2/83‑1 (MQ=255)

                                                                     attggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                              >  1:1771392‑M2/16‑99 (MQ=255)

                                                                      ttggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                             >  2:1937804‑M2/15‑99 (MQ=255)

                                                                      ttggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                             <  2:2120408‑M2/85‑1 (MQ=255)

                                                                      aaaCGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                              <  1:1227849/99‑1 (MQ=255)

                                                                        ggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                           <  1:1954335‑M2/87‑1 (MQ=255)

                                                                        ggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                           <  1:2315308‑M2/87‑1 (MQ=255)

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

                                                                          taccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:1334565‑M2/89‑1 (MQ=255)

                                                                            ccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                       >  1:1995757‑M2/9‑99 (MQ=255)

                                                                             cgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa                      >  1:1815666‑M2/8‑99 (MQ=255)

                                                                             cgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa                      >  2:1247645‑M2/8‑99 (MQ=255)

                                                                             cgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa                      <  2:1569294‑M2/92‑1 (MQ=255)

                                                                               agctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    >  1:939517‑M2/6‑99 (MQ=255)

                                                                               agctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    <  1:1454675‑M2/94‑1 (MQ=255)

                                                                               agctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    >  1:1727825‑M2/6‑99 (MQ=255)

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

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

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

                                                                                    gAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgg               >  1:2092271/1‑99 (MQ=255)

                                                                                     aaTTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  1:89365/1‑99 (MQ=255)

                                                                                     aaTTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  1:1275981/1‑99 (MQ=255)

                                                                                      attcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             <  2:1702724/96‑1 (MQ=255)

                                                                                       ttcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  2:1613375/3‑99 (MQ=255)

                                                                                           tAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:684261/1‑99 (MQ=255)

                                                                                           tAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  1:2402934/1‑99 (MQ=255)

                                                                                           tAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:1292916/1‑99 (MQ=255)

                                                                                           tAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  1:984933/1‑99 (MQ=255)

                                                                                             aTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      <  2:775322/99‑1 (MQ=255)

                                                                                             aTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      <  2:2028376/99‑1 (MQ=255)

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398

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