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	98.2%	392.2 / ‑4.7	111	K424E (AAA→GAA)	tyrS	tyrosyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (57/51); total (58/53)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

GTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGG  >  NC_000913/1715857‑1716048

gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCGTTTCCACTTAATGCATTATTCCCa                                                                                              >  2:144699/1‑100 (MQ=255)

gTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCa                                                                                              <  1:865832/100‑1 (MQ=255)

 tGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAg                                                                                             <  2:139236/100‑1 (MQ=255)

  gagaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGc                                                                                            <  2:593961/100‑1 (MQ=255)

   agaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTACAAGCa                                                                                           >  2:2238200/1‑100 (MQ=255)

    gaCTGGATAGCGAGAATATTTTTTTTTTTTTTTCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCaa                                                                                          >  2:608937/1‑100 (MQ=255)

    gaCTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCaa                                                                                          >  1:1036556/1‑100 (MQ=255)

       tGGATAGCGAGAATATTTTTCATCATGTTACCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCTCCTTTCCACTTAATGCATTATTCCCAGCAAATc                                                                                       <  2:1407773/100‑1 (MQ=255)

            ngCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACa                                                                                  >  1:670295/2‑100 (MQ=255)

             gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAg                                                                                 >  1:217112/1‑100 (MQ=255)

             gCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAg                                                                                 >  2:699204/1‑100 (MQ=255)

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

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

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

               gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGt                                                                                >  2:1727982/1‑99 (MQ=255)

               gagaATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTa                                                                               <  2:1704956/100‑1 (MQ=255)

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

                    tatTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCt                                                                          <  1:608938/100‑1 (MQ=255)

                     atTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCtt                                                                         <  2:195164/100‑1 (MQ=255)

                     atTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCtt                                                                         >  1:1833519/1‑100 (MQ=255)

                      tttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCttt                                                                        >  2:1736097/1‑100 (MQ=255)

                       ttttCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGATAGTAATTCtttt                                                                       <  1:1539085/100‑1 (MQ=255)

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

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

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

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

                           catcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTAcc                                                                   <  1:2126536/100‑1 (MQ=255)

                            atcatGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcg                                                                  >  2:188570/1‑100 (MQ=255)

                              catGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACcgcg                                                                >  2:869933/1‑100 (MQ=255)

                                tGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGAc                                                              <  2:110660/100‑1 (MQ=255)

                                tGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGAc                                                              <  2:572447/100‑1 (MQ=255)

                                 gTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTGCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACg                                                             >  2:1447040/1‑100 (MQ=255)

                                    cccTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAg                                                          >  1:988675/1‑100 (MQ=255)

                                      cTGTATGAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTCCCACTTAGTGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGt                                                         <  2:1261170/99‑1 (MQ=255)

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

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

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

                                           taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                   <  2:296578/100‑1 (MQ=255)

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

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

                                           taAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTa                                                   <  1:1433181/100‑1 (MQ=255)

                                             aaaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaa                                                 >  2:263099/1‑100 (MQ=255)

                                              aaaCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                <  2:1746068/100‑1 (MQ=255)

                                               ggaaaGGGGGAGTGAGAAATCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTaaaa                                                >  1:1740465/5‑99 (MQ=37)

                                                aCCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACg                                              >  2:1018400/1‑100 (MQ=255)

                                                 ccAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGa                                             <  1:847018/100‑1 (MQ=255)

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

                                                     ggggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaa                                         >  1:1861875/1‑100 (MQ=255)

                                                      gggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaaa                                        >  1:94668/1‑100 (MQ=255)

                                                      gggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCaaa                                        >  1:1319055/1‑100 (MQ=255)

                                                      gggAGTGATTACTCACT‑CCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                       <  2:1534822/100‑1 (MQ=255)

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

                                                       ggAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAAc                                       >  1:274051/1‑100 (MQ=255)

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

                                                        tcttccgatctcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGCAATTCTTTTTCCCGCGACGCAGTAAGGTTAAACGACCaaata                                      <  1:1576234/92‑3 (MQ=255)

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

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

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

                                                          gTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGa                                    <  1:188572/100‑1 (MQ=255)

                                                           tGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGAc                                   >  1:640336/1‑100 (MQ=255)

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

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

                                                             agaaatcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAATAAGGTAAAACAACCAAACAgccaa                                 >  1:1146085/6‑96 (MQ=255)

                                                             aTTTCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCCGACAGTACTTCTTTTTACCGCGACGCAGTAAGGTAAACCGCCCCAACAGACGa                                 >  1:1363808/1‑100 (MQ=255)

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

                                                              tttCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                <  1:1305213/100‑1 (MQ=255)

                                                              tttCTCACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAt                                >  1:446725/1‑100 (MQ=255)

                                                                 ctcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                             >  2:1378566/1‑100 (MQ=255)

                                                                 ctcACTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtctt                             >  2:1835318/1‑100 (MQ=255)

                                                                   cACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttct                           <  2:442492/100‑1 (MQ=255)

                                                                     cTCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGaagatcgga                         >  2:29131/1‑92 (MQ=255)

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

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

                                                                      tCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttct                        <  2:1863170/100‑1 (MQ=255)

                                                                      tCCCCCTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGAtcttcttct                        <  1:1447044/100‑1 (MQ=255)

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

                                                                       cccccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATcttcttctt                       <  1:1852668/100‑1 (MQ=255)

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

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

                                                                         cccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTa                     <  1:586931/100‑1 (MQ=255)

                                                                          ccTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaa                    <  2:1920626/100‑1 (MQ=255)

                                                                           cTTTCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTaaa                   <  1:1041290/100‑1 (MQ=255)

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

                                                                            tttCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaag                  <  2:2312327/100‑1 (MQ=255)

                                                                            tttCCACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaag                  <  2:1644312/100‑1 (MQ=255)

                                                                               ccACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag               <  2:1861869/100‑1 (MQ=255)

                                                                               ccACTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAaagaag               >  2:585680/1‑100 (MQ=255)

                                                                                cACTTAATGCATTATTCCCAGCAAATCAGAGAGTAATTCTTTTTACGGCGACGCAGAAAGGTAAAACGACCAAACAGACGATCTTGTTCTTTAAAGAAGt              >  1:2239985/1‑100 (MQ=255)

                                                                                 aCTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTa             >  1:1520975/1‑100 (MQ=255)

                                                                                 aCTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTa             >  2:1770444/1‑100 (MQ=255)

                                                                                  cTTAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTAt            >  2:2305201/1‑100 (MQ=255)

                                                                                    tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc          >  1:1886335/1‑100 (MQ=255)

                                                                                    tAATGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTc          >  1:980429/1‑100 (MQ=255)

                                                                                     aaTGCATTATTCCCAGCCCATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa         <  1:1018402/100‑1 (MQ=255)

                                                                                     aaTGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCa         <  2:1119707/100‑1 (MQ=255)

                                                                                       tGCATTATTCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAgg       >  2:35938/1‑100 (MQ=255)

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

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

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

                                                                                            tattCCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCgg  >  1:884240/1‑100 (MQ=255)

GTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCTGTATAAAAACCAGGGGGAGTGATTTCTCACTCCCCCTTTCCACTTAATGCATTATTTCCAGCAAATCAGACAGTAATTCTTTTTACCGCGACGCAGTAAGGTAAAACGACCAAACAGACGATCTTCTTCTTTAAAGAAGTATTCAGGATCGG  >  NC_000913/1715857‑1716048

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: `‑`