BRESEQ :: Evidence

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

Only 100 of 917 total aligned reads displayed.

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398

ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                                    >  2:689590‑M2/85‑99 (MQ=255)

gCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATg                                                                                                    <  1:1255636/99‑1 (MQ=255)

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

      gggggATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATa                                                                                              <  2:579138/99‑1 (MQ=255)

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

        gggATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGc                                                                                            >  1:1293506/1‑99 (MQ=255)

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

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

         gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTg                                                                                          >  2:2181514‑M2/76‑99 (MQ=255)

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

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

           gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGtt                                                                                        >  2:1661484‑M2/74‑99 (MQ=255)

           gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGtt                                                                                        >  1:2553778‑M2/74‑99 (MQ=255)

            tgtgCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATGGCTGtt                                                                                        <  2:256796/99‑1 (MQ=255)

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

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

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

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

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

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

                    gtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCtgtgtg                                                                               >  2:1556234‑M2/65‑99 (MQ=255)

                     tgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGa                                                                              >  1:1793399‑M2/64‑99 (MQ=255)

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

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

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

                         gATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAAt                                                                           <  2:1486010/99‑1 (MQ=255)

                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTg                                                                         <  1:2795430‑M2/41‑1 (MQ=255)

                          aTTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAAtt                                                                          <  2:2523022/99‑1 (MQ=255)

                           ttgtgaaggaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTTAAATTGt                                                                        <  2:815637‑M2/42‑1 (MQ=255)

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

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

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

                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     >  1:1001491‑M2/55‑99 (MQ=255)

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

                                ttGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                                    <  1:2056725/99‑1 (MQ=255)

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

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

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

                                    aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               >  1:1366315‑M2/49‑99 (MQ=255)

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

                                    aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               >  1:1495749‑M2/49‑99 (MQ=255)

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

                                     aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                              >  1:983261‑M2/48‑99 (MQ=255)

                                      aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                             >  1:1195845‑M2/47‑99 (MQ=255)

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

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

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

                                           cAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                                         >  1:2104545/1‑99 (MQ=255)

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

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

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

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

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

                                                ttttCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    <  1:2652107/99‑1 (MQ=255)

                                                   tCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                 >  2:268441/1‑99 (MQ=255)

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

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

                                                       gtcaagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            >  1:146298‑M2/30‑99 (MQ=255)

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

                                                          aagttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                                         <  2:1588961‑M2/73‑1 (MQ=255)

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

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

                                                            gttgtcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       <  1:1606600‑M2/75‑1 (MQ=255)

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

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

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

                                                                tcataattggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                                   >  2:2252286‑M2/21‑99 (MQ=255)

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

                                                                   taattggtaccgagctcGAATTCGTAATCGT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                <  1:530085‑M2/82‑1 (MQ=255)

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

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

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

                                                                       tggtaccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                            >  2:1397801‑M2/14‑99 (MQ=255)

                                                                          taccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         >  1:1414981‑M2/11‑99 (MQ=255)

                                                                          taccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         >  2:104173‑M2/11‑99 (MQ=255)

                                                                          gacgGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                          >  2:2859477/1‑99 (MQ=255)

                                                                            ccgagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                       >  2:1863302‑M2/9‑99 (MQ=255)

                                                                              gagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa                     >  2:1738063‑M2/7‑99 (MQ=255)

                                                                              gagctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa                     >  2:75349‑M2/7‑99 (MQ=255)

                                                                               ccAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa                     >  1:1924405/1‑99 (MQ=255)

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

                                                                                 ctcGAATTCGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc                  >  1:2616291‑M2/4‑99 (MQ=255)

                                                                                   tGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt                 <  1:2134931/99‑1 (MQ=255)

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

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

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

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

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

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

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

                                                                                      attcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             >  2:343149/4‑99 (MQ=255)

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

                                                                                       ttcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  2:393412/97‑1 (MQ=255)

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

                                                                                       ttcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  1:2182503/97‑1 (MQ=255)

                                                                                       ttcGTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  1:1133101/3‑99 (MQ=255)

                                                                                          gTAATCAT‑GTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         >  1:1971901/1‑99 (MQ=255)

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

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

                                                                                               cATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAAt     <  2:2429933/99‑1 (MQ=255)

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398

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