BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	366,300	Δ1 bp	81.5%	coding (6/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	.	81.5%	1960.9 / inf	655	coding (6/3075 nt)	lacZ	beta‑galactosidase
Reads supporting (aligned to +/- strand): ref base G (56/65); new base . (280/254); total (336/319)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 685 total aligned reads displayed.

GAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366204‑366398

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

  aaGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCa                                                                                                >  1:1518254/1‑99 (MQ=255)

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

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

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

       gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTg                                                                                          <  1:251402‑M2/24‑1 (MQ=255)

        agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                                         >  2:1951227‑M2/75‑99 (MQ=255)

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

        agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                                         >  2:2308115‑M2/75‑99 (MQ=255)

        agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGt                                                                                         >  2:473329‑M2/75‑99 (MQ=255)

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

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

         gcgagttccgtgccggttatgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGtt                                                                                        <  1:1186209‑M2/26‑1 (MQ=255)

         aTGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGt                                                                                         >  2:1803933/1‑99 (MQ=255)

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

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

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

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

              ctgcAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCTCGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCt                                                                                    <  1:1896547/99‑1 (MQ=255)

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

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

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

                        gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  2:251209‑M2/59‑99 (MQ=255)

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

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

                         ttAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  2:558246/1‑99 (MQ=255)

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

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

                             gTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     >  1:2087740/1‑99 (MQ=255)

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

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

                              aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   >  1:875270‑M2/53‑99 (MQ=255)

                               agaaaaagtgaatgatgtaaccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                                  >  1:2339510‑M2/52‑99 (MQ=255)

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

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

                                   tAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               >  1:953092/1‑99 (MQ=255)

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

                                   aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                              >  2:762627‑M2/48‑99 (MQ=255)

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

                                     agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            <  1:834313‑M2/54‑1 (MQ=255)

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

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

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

                                        ccAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:2167295/1‑99 (MQ=255)

                                          aGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        <  2:1571922/99‑1 (MQ=255)

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

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

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

                                              tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   >  1:268439‑M2/37‑99 (MQ=255)

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

                                                tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                 >  2:2474795‑M2/35‑99 (MQ=255)

                                                tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                 <  1:745180‑M2/65‑1 (MQ=255)

                                                 tCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                 <  1:2297571/99‑1 (MQ=255)

                                                 agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                                >  1:2210794‑M2/34‑99 (MQ=255)

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

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

                                                    cAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCAAAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACAt                                              <  2:1170985/99‑1 (MQ=255)

                                                     gtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            <  2:1694369‑M2/70‑1 (MQ=255)

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

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

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

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

                                                           gacgTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       <  2:628776/99‑1 (MQ=255)

                                                               ttGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                                   >  1:2420015/1‑99 (MQ=255)

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

                                                               cataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                                  >  2:2364177‑M2/20‑99 (MQ=255)

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

                                                                  tAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                <  2:148248/99‑1 (MQ=255)

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

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

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

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

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

                                                                         accgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                        >  1:1321102‑M2/10‑99 (MQ=255)

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

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

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

                                                                               aGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  1:2260811/99‑1 (MQ=255)

                                                                               aGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  1:11234/99‑1 (MQ=255)

                                                                                tcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt                 >  2:2136491‑M2/3‑99 (MQ=255)

                                                                                tcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt                 >  1:1538022‑M2/3‑99 (MQ=255)

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

                                                                                  gAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg                <  2:1909351/99‑1 (MQ=255)

                                                                                    attcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             >  1:644785/4‑99 (MQ=255)

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

                                                                                    aTCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              <  1:1795253/99‑1 (MQ=255)

                                                                                    aTCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  2:1726185/1‑99 (MQ=255)

                                                                                     tCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             >  2:1935072/1‑99 (MQ=255)

                                                                                     tCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             >  2:268315/1‑99 (MQ=255)

                                                                                       cGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGc          >  2:1997907/1‑99 (MQ=255)

                                                                                        gTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         <  1:1591501/99‑1 (MQ=255)

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

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

                                                                                          aaTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTa       >  1:718855/1‑99 (MQ=255)

                                                                                          aaTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTa       >  2:2054313/1‑99 (MQ=255)

                                                                                           aTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  2:1536047/1‑99 (MQ=255)

                                                                                           aTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  1:1242519/1‑99 (MQ=255)

                                                                                           aTCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  2:1449502/1‑99 (MQ=255)

                                                                                               tGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAtga   >  1:522693/1‑99 (MQ=255)

                                                                                                ggTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAtgag  <  1:1999742/99‑1 (MQ=255)

GAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366204‑366398

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