BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,953	T→C	K424E (AAA→GAA)	tyrS ←	tyrosyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,953	0	T	C	97.0%	337.9 / 0.0	101	K424E (AAA→GAA)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (49/48); total (51/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01

Only 100 of 102 total aligned reads displayed.

CGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCG  >  NC_000913/1715856‑1716047

cGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTccc                                                                                              <  1:2042302/100‑1 (MQ=255)

  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCTCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAg                                                                                            >  1:1444334/1‑100 (MQ=255)

  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCan                                                                                            <  1:1650010/100‑2 (MQ=255)

  tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAg                                                                                            >  1:1578603/1‑100 (MQ=255)

   gagaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGc                                                                                           >  2:1297485/1‑100 (MQ=255)

       ttccgatctgagaATATTTTTCATCATGTCCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAAt                                                                                       <  1:1892974/91‑1 (MQ=255)

       cTGGATAGCGAGAATATTTTTCATCATGTTCGCTGTATAAAAGCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAAt                                                                                       >  1:1379741/1‑100 (MQ=255)

        tGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATc                                                                                      <  1:987157/100‑1 (MQ=255)

            tAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGAc                                                                                  <  1:178937/100‑1 (MQ=255)

              gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAg                                                                                <  1:113460/100‑1 (MQ=255)

              gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCAATCCCCCTTTCCACTTAATGCATTATTCCAAGCAAATCAGACAg                                                                                >  1:1218459/1‑100 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCCCTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                               >  2:929512/1‑100 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                               >  2:1863485/1‑100 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                               >  1:519176/1‑100 (MQ=255)

               cGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                               <  1:399169/100‑1 (MQ=255)

                gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTa                                                                              >  1:1695626/1‑100 (MQ=255)

                gagaATATTTTTCATCATGTCCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAAtagatcggaa                                                                              >  2:1892956/1‑91 (MQ=255)

                  gaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAAt                                                                            >  2:1179266/1‑100 (MQ=255)

                  gaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAAt                                                                            <  2:1756564/100‑1 (MQ=255)

                     tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAAtt                                                                           >  2:1742353/1‑98 (MQ=255)

                     tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCt                                                                         >  1:1806094/1‑100 (MQ=255)

                        ttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCtttt                                                                      <  2:122827/100‑1 (MQ=255)

                        ttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCtttt                                                                      >  2:763666/1‑100 (MQ=255)

                         tttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCttttt                                                                     >  1:343808/1‑100 (MQ=255)

                         tttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCttttt                                                                     <  2:152506/100‑1 (MQ=255)

                           tcatcaTGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                   <  2:1898322/100‑1 (MQ=255)

                           tcatcaTGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAc                                                                   >  1:1096698/1‑100 (MQ=255)

                             atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAcct                                                                 >  2:1402257/1‑99 (MQ=255)

                             atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                 <  2:1319896/100‑1 (MQ=255)

                             atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                 <  1:1298929/100‑1 (MQ=255)

                             atcatGTTCCCCGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                 <  2:123071/100‑1 (MQ=255)

                              tcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTCATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                 >  1:1143586/1‑99 (MQ=255)

                                atGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGa                                                              >  1:1577724/1‑100 (MQ=255)

                                    tCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCa                                                          >  1:1713722/1‑100 (MQ=255)

                                       cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                       <  2:1578593/100‑1 (MQ=255)

                                       cTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTa                                                       <  2:1806078/100‑1 (MQ=255)

                                        tGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCCACGCAGTa                                                       >  2:1694845/1‑99 (MQ=255)

                                          tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCATATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                    <  2:1096693/100‑1 (MQ=255)

                                          tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                    <  2:1713709/100‑1 (MQ=255)

                                          tataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAgg                                                    >  2:1325235/1‑100 (MQ=255)

                                           ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                   >  2:296589/1‑100 (MQ=255)

                                           ataAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGt                                                   <  2:2044182/100‑1 (MQ=255)

                                            taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                  >  1:166002/1‑100 (MQ=255)

                                             aaaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaa                                                 >  2:1875663/1‑100 (MQ=255)

                                                aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAAc                                              <  1:981829/100‑1 (MQ=255)

                                                aaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAAc                                              <  1:1182530/100‑1 (MQ=255)

                                                 acaaGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACg                                             <  2:244922/97‑1 (MQ=255)

                                                 aCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACg                                             <  2:1043876/100‑1 (MQ=255)

                                                  ccAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACACAGTAAGGTAAAACGa                                            <  2:1405596/100‑1 (MQ=255)

                                                   cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                           <  1:1297492/100‑1 (MQ=255)

                                                   cAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAc                                           <  2:185030/100‑1 (MQ=255)

                                                     gggtgAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTCTTACCGCGACGCAGTAAGGTAAAACGACCa                                         <  2:1218453/100‑1 (MQ=255)

                                                     gggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                         >  1:1021193/1‑100 (MQ=255)

                                                     gggggAGTGATTTATCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCa                                         >  1:162219/1‑100 (MQ=255)

                                                      ggggAGTGNTTTCTCACTCCCCCTTTNCACTNAANGCATTATTCNCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAANGACCaa                                        >  1:1584641/1‑100 (MQ=255)

                                                        ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                      <  2:427175/100‑1 (MQ=255)

                                                        ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                      <  1:763669/100‑1 (MQ=255)

                                                           gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGa                                   >  1:1237744/1‑100 (MQ=255)

                                                           ctGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACatg                                   <  1:190354/99‑3 (MQ=255)

                                                            tgagaaatcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                  <  2:1565858/93‑1 (MQ=255)

                                                             gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACat                                    >  2:190353/1‑96 (MQ=255)

                                                             gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACg                                 >  2:1297951/1‑100 (MQ=255)

                                                             gATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACg                                 <  2:882709/100‑1 (MQ=255)

                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGa                                <  1:1863502/100‑1 (MQ=255)

                                                              aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGa                                >  2:188434/1‑100 (MQ=255)

                                                                ttCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtc                              <  1:235178/100‑1 (MQ=255)

                                                                   acACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttc                           >  1:260186/2‑100 (MQ=255)

                                                                    cACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCCGACAGTAATTCTTTTTACTGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                          >  2:1426565/1‑100 (MQ=255)

                                                                    cACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                          <  2:162218/100‑1 (MQ=255)

                                                                     aCTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttctt                         >  1:1028073/1‑100 (MQ=255)

                                                                      cTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                        >  1:1789722/1‑100 (MQ=255)

                                                                      cTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                        >  1:1712988/1‑100 (MQ=255)

                                                                      cTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttc                        <  1:73785/100‑1 (MQ=255)

                                                                        cccccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttct                       <  2:1695613/99‑1 (MQ=255)

                                                                        cccccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                      >  2:2008243/1‑100 (MQ=255)

                                                                        cccccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                      >  2:1458258/1‑100 (MQ=255)

                                                                         ccccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCttt                     <  1:1425265/100‑1 (MQ=255)

                                                                          cccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                    >  2:504026/1‑100 (MQ=255)

                                                                          cccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                    >  1:365981/1‑100 (MQ=255)

                                                                           ccTTTCCACTTAATGCATTATTCCCAGAAAATCAGACAGTAATTCTTTTGTCCGCGTCTCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                   <  1:134143/100‑1 (MQ=255)

                                                                            cTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                  >  1:64956/1‑100 (MQ=255)

                                                                             tttCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaag                 <  2:1444325/100‑1 (MQ=255)

                                                                             tttCCACTTAATGCATTATTCCCAGCACATCAGACAGTAATTCTTTTTACCGCGACTCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                   <  2:1662370/98‑1 (MQ=255)

                                                                              ttCCACTTAATGCATTATTCCCTGCAAATCACACAGTAATTCTTTTTGCCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaaga                <  2:1379732/100‑1 (MQ=255)

                                                                              ttCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaaga                <  2:1021188/100‑1 (MQ=255)

                                                                               tCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa               >  1:1336572/1‑100 (MQ=255)

                                                                               tCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaa               >  1:903681/1‑100 (MQ=255)

                                                                                ccACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag              >  1:1395692/1‑100 (MQ=255)

                                                                                ccACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGAACAAAAAAACGATCTTCTTCTTTAaaaaag              >  1:1848222/1‑100 (MQ=255)

                                                                                 cACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGt             <  2:343807/100‑1 (MQ=255)

                                                                                     tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc         <  1:296590/100‑1 (MQ=255)

                                                                                     tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc         <  2:1015970/100‑1 (MQ=255)

                                                                                       aTGCATTATTGCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg       <  2:1336563/100‑1 (MQ=255)

                                                                                       aTGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg       <  1:1613116/100‑1 (MQ=255)

                                                                                       aTGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAg       <  1:1612691/100‑1 (MQ=255)

                                                                                         gCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa     >  1:1175272/1‑100 (MQ=255)

                                                                                          cATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGa     >  1:576874/1‑99 (MQ=255)

                                                                                           attattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATc   >  1:548331/1‑100 (MQ=255)

                                                                                           attattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATc   <  2:1579122/100‑1 (MQ=255)

                                                                                            ttattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCg  <  1:1875680/100‑1 (MQ=255)

CGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCG  >  NC_000913/1715856‑1716047

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