BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,597,313	(C)_7→8	pseudogene (675/3861 nt)	yneO ←	pseudogene, AidA homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,597,306	1	.	C	96.2%	371.1 / 6.4	107	pseudogene (682/3861 nt)	yneO	pseudogene, AidA homolog
Reads supporting (aligned to +/- strand): ref base . (1/3); new base C (49/53); total (50/57)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400

agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa                                                                                                   >  1:803888/1‑100 (MQ=255)

 gagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑c                                                                                                >  2:2360131/1‑100 (MQ=255)

  agTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaaccaa‑‑cc                                                                                               >  2:1781105/1‑100 (MQ=255)

  agTCTGCGAACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCAACCGATAATCGTATTTTCAGCCGTaaccaa‑‑cc                                                                                               >  2:471107/1‑100 (MQ=255)

   gTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑c                                                                                             <  2:176694/100‑1 (MQ=255)

     cTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑ccc                                                                                           >  1:1338551/1‑100 (MQ=255)

      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccc                                                                                          <  2:803894/100‑1 (MQ=255)

      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccc                                                                                          >  1:291914/1‑100 (MQ=255)

      tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccc                                                                                          >  2:2323432/1‑100 (MQ=255)

          tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccccctc                                                                                      >  1:1277067/1‑97 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccccctcg                                                                                     <  1:1195010/100‑5 (MQ=255)

           aCCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑cccccctcg                                                                                     >  2:2538403/1‑96 (MQ=255)

             cGTGATGACATAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATTGTCTTTTCAGCCGTAACCAAC‑CC‑CCCCCTcgcg                                                                                    <  1:1903219/100‑1 (MQ=255)

             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTcgcg                                                                                    <  1:2770198/100‑1 (MQ=255)

             cGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTcgcg                                                                                    <  2:1338569/100‑1 (MQ=255)

              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGa                                                                                   <  2:1700346/100‑1 (MQ=255)

              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGa                                                                                   >  2:1467506/1‑100 (MQ=255)

              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGa                                                                                   <  2:2632523/100‑1 (MQ=255)

              gTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑cc‑ccccccctcgcg                                                                                  <  2:2764044/100‑8 (MQ=255)

                 atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATa                                                                                >  1:944566/1‑100 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCGCATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAg                                                                               <  1:93847/100‑1 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAg                                                                               >  2:2022811/1‑100 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAg                                                                               <  2:2439374/100‑1 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCGCGCGAATAg                                                                               <  2:1237639/100‑1 (MQ=255)

                  tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CC‑CCCCCTCGCGAATAGt                                                                              >  2:967236/1‑100 (MQ=255)

                   gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCC‑CCCCCTCGCGAATAg                                                                               <  1:781879/100‑1 (MQ=255)

                   gaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CC‑CCCCCTCGCGAATAGTc                                                                             <  1:403480/100‑1 (MQ=255)

                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCg                                                                            <  1:1754877/100‑1 (MQ=255)

                     cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCg                                                                            >  2:2587658/1‑100 (MQ=255)

                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTcccc                                                                     <  1:2531107/100‑1 (MQ=255)

                            ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTcccc                                                                     >  2:2711725/1‑100 (MQ=255)

                             tGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTcccc                                                                     >  2:1062788/1‑99 (MQ=255)

                              gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAg                                                                   >  2:754072/1‑100 (MQ=255)

                              gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATGGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAg                                                                   <  1:1946454/100‑1 (MQ=255)

                                 ccTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTc                                                                <  2:1711221/100‑1 (MQ=255)

                                   tGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCt                                                              <  2:1916201/100‑1 (MQ=255)

                                    gCACATTCAGGGTTCCGATACCGGTGACATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCCCCCCAGTTCCtt                                                             >  2:1767586/1‑100 (MQ=255)

                                      acaTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                           <  2:121744/100‑1 (MQ=255)

                                        aTTCAGGGTTCCGATCCCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                           <  1:2002858/98‑1 (MQ=255)

                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGACTAGTCGCCTCCCCCGTCCCt                                                              >  1:691205/1‑95 (MQ=255)

                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  2:646255/1‑100 (MQ=255)

                                        aTTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGAtt                                                         >  1:2576105/1‑100 (MQ=255)

                                          tCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAc                                                       <  1:113991/100‑1 (MQ=255)

                                           cAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTAcc                                                      >  1:2487364/1‑100 (MQ=255)

                                             gggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCaa                                                    <  1:2841303/100‑1 (MQ=255)

                                              ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                   >  2:90918/1‑100 (MQ=255)

                                              ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                   >  1:838074/1‑100 (MQ=255)

                                                ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:967227/100‑1 (MQ=255)

                                                 tCCGATAGCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGAGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTg                                                <  1:1561230/100‑1 (MQ=255)

                                                 acgattatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCCGTTCCTTGATTACCAATTTg                                                >  2:5803/9‑100 (MQ=255)

                                                  cgattatcGGTGGCAATGCCACCGATAATCGTATTTTTTTTTTTTTTTAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  1:703151/8‑100 (MQ=255)

                                                  ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  1:2846408/1‑100 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              >  2:1598468/1‑100 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              >  2:1463560/1‑100 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              <  1:1038650/100‑1 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaa                                              <  2:13464/100‑1 (MQ=255)

                                                   cGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                               >  2:1598342/1‑99 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             <  2:1500936/100‑1 (MQ=255)

                                                    gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  2:1622893/1‑100 (MQ=255)

                                                    attatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGaaa                                             >  2:1255856/6‑100 (MQ=255)

                                                     ttatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAt                                            >  1:1138074/5‑100 (MQ=255)

                                                      tatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAAtt                                           >  2:1674948/4‑100 (MQ=255)

                                                       atcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                          >  2:772339/3‑100 (MQ=255)

                                                        ccGGTGGCATTGCCACCAATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCa                                         >  1:5851/1‑100 (MQ=255)

                                                           gTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                      >  2:1600522/1‑100 (MQ=255)

                                                            tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTg                                     <  2:2581419/100‑1 (MQ=255)

                                                             ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                    >  2:2096590/1‑100 (MQ=255)

                                                             ggCATGGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CCTCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                    <  1:1231696/100‑1 (MQ=255)

                                                              gCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                   <  1:643251/100‑1 (MQ=255)

                                                               caatGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:772333/97‑1 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  2:604200/100‑1 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:1318694/100‑1 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                  <  1:2415748/100‑1 (MQ=255)

                                                               cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                   <  1:2030846/99‑1 (MQ=255)

                                                                    ccACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCa                             >  2:2740030/1‑100 (MQ=255)

                                                                     cACCGATAATCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCTCAGTTCCTTGCTTACCAATTTGAAATTCTATTGAGGAATctt                            <  2:1960922/100‑3 (MQ=255)

                                                                            aaTCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATATTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAttttt                     >  1:2426349/1‑100 (MQ=255)

                                                                             aTCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    <  1:2750406/100‑1 (MQ=255)

                                                                             aTCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttttt                    >  2:144080/1‑100 (MQ=255)

                                                                              tCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  1:1861923/100‑1 (MQ=255)

                                                                              tCGTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                   <  2:703157/100‑1 (MQ=255)

                                                                                gTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTTCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  1:991136/1‑100 (MQ=255)

                                                                                gTATTTTCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatt                 >  2:302397/1‑100 (MQ=255)

                                                                                   ttttCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  1:951951/100‑1 (MQ=255)

                                                                                   ttttCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:2846446/100‑1 (MQ=255)

                                                                                   ttttCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGc              <  2:5851/100‑1 (MQ=255)

                                                                                    tttCCGCCGGAACCAAACCC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAg               <  1:466064/99‑1 (MQ=255)

                                                                                    tttCAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGcc             >  2:1680058/1‑100 (MQ=255)

                                                                                       cAGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCCCt          >  1:2003601/1‑100 (MQ=255)

                                                                                        aGCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTc         >  1:248771/1‑100 (MQ=255)

                                                                                         gCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        <  1:1670501/100‑1 (MQ=255)

                                                                                         gCCGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCg        >  1:151800/1‑100 (MQ=255)

                                                                                          ccGTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGc       <  2:1785091/100‑1 (MQ=255)

                                                                                            gTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCa     >  1:682609/1‑100 (MQ=255)

                                                                                            gTAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCa     >  2:902310/1‑100 (MQ=255)

                                                                                             tAACCAAC‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAc    <  1:194097/100‑1 (MQ=255)

                                                                                              aaccaac‑cC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAcc   >  1:2130234/1‑100 (MQ=255)

                                                                                               accaac‑cC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  2:274670/100‑1 (MQ=255)

                                                                                               accaac‑cC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  <  2:274164/100‑1 (MQ=255)

                                                                                               accaac‑cC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCAccc  >  1:2846664/1‑100 (MQ=255)

AGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAA‑‑CC‑CCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCC  >  NC_000913/1597207‑1597400

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