BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,311,831	0	C	T	20.9%	451.9 / 115.4	221	intergenic (+23/+17)	tonB/yciA	membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand): ref base C (85/89); new base T (24/22); total (110/111)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 227 total aligned reads displayed.

TGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913_3_bme_pgi/1311733‑1311930

tGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                                                                    >  1:1853003/1‑100 (MQ=255)

  aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGACATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   >  1:2151843/1‑99 (MQ=255)

  aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   <  1:2716900/99‑1 (MQ=255)

  aGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   <  2:2059645/99‑1 (MQ=255)

  aGCCGGGTAAGCCAGGCAGTGGGATTGCGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                  <  2:2746396/100‑1 (MQ=255)

   gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCg                                                                                                   >  2:1642416/1‑98 (MQ=255)

     cGGGTAAGCCAGGCAGTGGGAATGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAACGCCTCCGACCGGAgg                                                                                               >  2:2669880/1‑100 (MQ=255)

      gggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGc                                                                                              <  1:465497/100‑1 (MQ=255)

       ngTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                                                                             >  1:3111021/2‑100 (MQ=255)

       ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                                                                             <  1:534291/100‑1 (MQ=255)

       ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCt                                                                                             >  1:1787064/1‑100 (MQ=255)

         tAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCttt                                                                                           <  1:900824/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTACAATTAACGGCACCACCGAAATTCAGTAAGCAGAACGTCAACAGCCTCCGccc                                                                                                    >  2:3650848/1‑89 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  2:1460900/99‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  1:3056645/99‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                          <  1:198863/99‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:2291221/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:589423/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:1043189/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:4762099/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  1:1169052/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                          <  1:3431940/99‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  2:2724265/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  1:2901212/100‑1 (MQ=255)

           aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                         <  1:570569/100‑1 (MQ=255)

            gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                                                                        <  1:121490/100‑1 (MQ=255)

            gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                                                                        <  1:3284372/100‑1 (MQ=255)

            gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                        >  1:1558323/1‑100 (MQ=255)

            gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                        >  1:3439376/1‑100 (MQ=255)

              cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                       >  1:2393307/1‑99 (MQ=255)

              cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATACAGTAAGCAGAAAGTCAAAAGCCTCCGCCCGGAGGCTTTTGAc                                                                                       >  2:5316692/1‑99 (MQ=255)

              cAGGCAGTGGGACTGTGGTCAATAACCTGTTTAAAACTAACGGCACCACCGACATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAcc                                                                                                    >  2:4343748/1‑86 (MQ=255)

               aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                         <  2:740936/96‑1 (MQ=255)

                  cAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTa                                                                                  >  2:2145578/1‑100 (MQ=255)

                   aGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTAc                                                                                 <  2:1085934/100‑1 (MQ=255)

                   aGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTACGACCGGAGGCTTCCGACTATTa                                                                                  >  2:5414915/1‑99 (MQ=255)

                    gTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                >  2:1757206/1‑100 (MQ=255)

                     tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                >  2:4805348/1‑99 (MQ=255)

                     tGGGATTGTGGGGAATGTCCTGTGTAAAAGTAACGGCACCACCGAAATGCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                               <  2:2247014/100‑1 (MQ=255)

                            gtggtgAATATCGTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCATAAGCCTCCGACCTGAGGCTTTTGACTATTACTCAACAgg                                                                         <  1:5088575/99‑1 (MQ=255)

                            gtggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTACGCAGAAAGCTAAAAGCCTCCGACCGGAGGCTACTGCCTCCTAATCAACAgc                                                                         >  2:3858410/1‑98 (MQ=255)

                             tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                                                                        >  2:4896191/1‑99 (MQ=255)

                             tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                       <  1:1298062/100‑1 (MQ=255)

                                tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAgg                                                                    >  1:3618520/1‑100 (MQ=255)

                                tgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTACTACTCaa                                                                             >  2:11827/1‑91 (MQ=255)

                                      tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                              <  1:3300649/100‑1 (MQ=255)

                                      tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                              <  2:4359080/100‑1 (MQ=255)

                                      tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAg                                                              >  1:4569416/1‑100 (MQ=255)

                                      tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACCCAACAGGTAAGGCGCGAg                                                              >  2:3613910/1‑100 (MQ=255)

                                          gTTTAAAATTAACGGCACCACCGAACTGCAGTAAGCAGAAAGCCAAAAGCCACCGACCGGAGGCCTTTGACCATTACCCaa                                                                             >  2:5098818/1‑81 (MQ=255)

                                          gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                          >  1:5199920/1‑100 (MQ=255)

                                           tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGtttt                                                         >  1:4132064/1‑100 (MQ=255)

                                             tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       >  1:1707225/1‑100 (MQ=255)

                                              aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTcc                                                       <  1:240448/99‑1 (MQ=255)

                                              aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCATTTGACTATTACTCAACAGCCAAGGCGCGAGGTTTTCCt                                                      >  2:1242290/1‑100 (MQ=255)

                                               aaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                                                                        <  2:2885434/65‑1 (MQ=255)

                                               aaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                     >  1:4305887/1‑100 (MQ=255)

                                               aaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                                      <  2:5170923/99‑1 (MQ=255)

                                                aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                     >  1:2710325/1‑99 (MQ=255)

                                                  ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                  <  1:4026231/100‑1 (MQ=255)

                                                    aaCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAGGGCGCGAGGTTTTCCGTCAGGa                                                >  1:250235/1‑100 (MQ=255)

                                                     aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGa                                                <  2:4656816/99‑1 (MQ=255)

                                                     aCGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                               <  1:1488524/100‑1 (MQ=255)

                                                      cGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                               <  2:2395983/99‑1 (MQ=255)

                                                      cGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAt                                               >  1:2762192/1‑99 (MQ=255)

                                                        gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCa                                             <  2:3475773/99‑1 (MQ=255)

                                                         caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAc                                           >  1:3562051/1‑100 (MQ=255)

                                                           ccaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                          <  2:1047312/99‑1 (MQ=255)

                                                             accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgcg                                       <  1:1144058/100‑1 (MQ=255)

                                                             accGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                        >  1:1434788/1‑99 (MQ=255)

                                                              ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                        >  2:160447/1‑98 (MQ=255)

                                                              ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                                      >  1:2864700/1‑100 (MQ=255)

                                                              ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                                      >  1:4295446/1‑100 (MQ=255)

                                                                      cAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACtt                               <  2:3711404/99‑1 (MQ=255)

                                                                       aGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTa                              <  1:3715750/99‑1 (MQ=255)

                                                                        gTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa                            >  1:1647632/1‑100 (MQ=255)

                                                                         tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                           <  2:4564601/100‑1 (MQ=255)

                                                                         tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                           <  2:4233427/100‑1 (MQ=255)

                                                                           aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataa                         <  2:2503532/100‑1 (MQ=255)

                                                                             gtaataGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                       >  2:25447/6‑100 (MQ=37)

                                                                                 aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATCCTTAAATAATGCtt                    >  2:4028432/1‑99 (MQ=255)

                                                                                 aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                   <  1:1539911/100‑1 (MQ=255)

                                                                                  aGTCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                  >  1:939734/1‑100 (MQ=255)

                                                                                   gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg                 <  1:4921791/100‑1 (MQ=255)

                                                                                   ctCAAAAGCCTCCGAACGGAGGCTTTTAACTATTACTCAACAGGTAAGCCGCGCGGTTTTCCTTCAGGATCAACCGCGACATAcct                               >  2:3015565/2‑84 (MQ=255)

                                                                                    tCAAAAGCCTCCGGTCGGAGGCTTTTGACTTTTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGCCATACTTAAATCATGCCTCTGt                >  2:5055090/1‑100 (MQ=37)

                                                                                    tCAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt                >  1:4305340/1‑100 (MQ=255)

                                                                                     naaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  1:4556466/2‑100 (MQ=255)

                                                                                     cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  1:716738/1‑100 (MQ=255)

                                                                                     cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc               >  2:5463315/1‑100 (MQ=255)

                                                                                      aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:4849504/1‑100 (MQ=255)

                                                                                      aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:1131555/1‑100 (MQ=255)

                                                                                       aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:1303130/1‑99 (MQ=255)

                                                                                       aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  1:4233343/1‑100 (MQ=255)

                                                                                       aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  1:3827743/1‑100 (MQ=255)

                                                                                       aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGGTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  1:3947780/1‑100 (MQ=39)

                                                                                       aaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGc             >  2:2617972/1‑100 (MQ=255)

                                                                                            cTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg              >  2:4922808/1‑94 (MQ=255)

                                                                                               cGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAg      <  2:4541231/99‑1 (MQ=255)

                                                                                               cGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGc     <  2:918170/100‑1 (MQ=255)

TGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTT  >  NC_000913_3_bme_pgi/1311733‑1311930

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