BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	366,290	C→T	93.5%	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	93.5%	9061.3 / inf	2565	D6N (GAT→AAT)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand): ref base C (72/95); new base T (1198/1199); total (1271/1294)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 2569 total aligned reads displayed.

ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGG  >  NC_000913/366192‑366387

  tgcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCg                                                                                                 <  1:1050895‑M2/7‑1 (MQ=255)

   gcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGt                                                                                                >  2:1031683‑M2/92‑99 (MQ=255)

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

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

          ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataatcggtaccgagctcGAATTCGTAATCATg                                                                                         >  1:558086‑M2/85‑99 (MQ=255)

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

            ccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCatgtc                                                                                       >  1:1911398‑M2/83‑97 (MQ=255)

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

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

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

              agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                     >  2:605005‑M2/81‑95 (MQ=255)

              agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                     >  2:908839‑M2/81‑95 (MQ=255)

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

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

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

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

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

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

                      tgtgCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGtt                                                                             <  1:858238/99‑1 (MQ=255)

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

                         gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCt                                                                         >  2:2175243‑M2/70‑99 (MQ=255)

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

                          ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtg                                                                        >  2:1581535‑M2/69‑99 (MQ=255)

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

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

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

                          ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtg                                                                        >  2:878084‑M2/69‑99 (MQ=255)

                             cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtgt                                                                     >  2:126009‑M2/66‑99 (MQ=255)

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

                              gtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtgtg                                                                    >  1:2319026‑M2/65‑99 (MQ=255)

                               tgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGa                                                                   <  2:1739916‑M2/36‑1 (MQ=255)

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

                                  cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAAt                                                                >  1:1147122‑M2/61‑99 (MQ=255)

                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGt                                                             <  2:2266068‑M2/42‑1 (MQ=255)

                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGt                                                             <  2:2424322‑M2/42‑1 (MQ=255)

                                     ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGt                                                             <  1:1279005‑M2/42‑1 (MQ=255)

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

                                      tgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGtt                                                            >  1:293824‑M2/57‑99 (MQ=255)

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

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

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

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

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

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

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

                                            gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                      >  2:1383452‑M2/51‑99 (MQ=255)

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

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

                                                aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                  >  2:1292247‑M2/47‑99 (MQ=255)

                                                 agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                 >  2:367098‑M2/46‑99 (MQ=255)

                                                   tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                               <  1:366413‑M2/56‑1 (MQ=255)

                                                    gaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                              >  1:1740791‑M2/43‑99 (MQ=255)

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

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

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

                                                      atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                            >  1:1673137‑M2/41‑99 (MQ=255)

                                                       tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                           <  2:523243‑M2/60‑1 (MQ=255)

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

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

                                                             agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                     >  2:470931‑M2/34‑99 (MQ=255)

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

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

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

                                                                cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                  <  1:499559‑M2/69‑1 (MQ=255)

                                                                 aGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                  <  2:1340996/99‑1 (MQ=255)

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

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

                                                                    aagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                              >  2:64515‑M2/27‑99 (MQ=255)

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

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

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

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

                                                                      gttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                            >  2:1550272‑M2/25‑99 (MQ=255)

                                                                          tcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                        >  1:522‑M2/21‑99 (MQ=255)

                                                                            ataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGc                      >  2:1707439‑M2/19‑99 (MQ=255)

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

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

                                                                              tAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                     >  1:989126/1‑99 (MQ=255)

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

                                                                                ttggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                  >  1:1092107‑M2/15‑99 (MQ=255)

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

                                                                                   gtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt               >  1:75618‑M2/12‑99 (MQ=255)

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

                                                                                     accgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt             <  1:1190420‑M2/90‑1 (MQ=255)

                                                                                      ccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa            <  1:899932‑M2/91‑1 (MQ=255)

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

                                                                                         agctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAtaaagtgtnnnn         <  2:1477769‑M2/94‑5 (MQ=255)

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

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

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

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

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

                                                                                          gctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc        >  1:1467285‑M2/5‑99 (MQ=255)

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

                                                                                           ctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc       >  2:328148‑M2/4‑99 (MQ=255)

                                                                                             cGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg     >  2:571032‑M2/2‑99 (MQ=255)

                                                                                             cGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg     >  1:1787598‑M2/2‑99 (MQ=255)

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

                                                                                              gAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgg    <  1:2400381/99‑1 (MQ=255)

                                                                                               aaTTCGTAATCAT‑GTCATATCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg   <  2:414574/99‑1 (MQ=255)

ACGCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGG  >  NC_000913/366192‑366387

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