BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,408,560	G→A	A144T (GCC→ACC)	yjfJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,408,560	0	G	A	99.0%	687.5 / NA	206	A144T (GCC→ACC)	yjfJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/2); major base A (109/93); minor base C (0/2); total (109/97)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

cTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAAcc                                                                                         <  2:2653961/93‑1 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       <  1:4084822/93‑1 (MQ=255)

  ggAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCAc                                                                                       >  1:2352759/1‑93 (MQ=255)

  ggAAGGCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCAc                                                                                       <  1:1317028/93‑1 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      <  2:597539/93‑1 (MQ=255)

   gAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACt                                                                                      <  1:3479194/93‑1 (MQ=255)

    aaGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTg                                                                                     >  1:2119402/1‑93 (MQ=255)

     aGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGa                                                                                    >  2:2232334/1‑93 (MQ=255)

     aGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGa                                                                                    <  1:939670/93‑1 (MQ=255)

     aGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGa                                                                                    >  1:4375582/1‑93 (MQ=255)

      gTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGaa                                                                                   <  1:3162847/93‑1 (MQ=255)

        ctcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGc                                                                                 <  2:1894006/93‑1 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                >  2:4246557/1‑93 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                <  1:2365829/93‑1 (MQ=255)

         tcCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGcc                                                                                >  2:3388847/1‑93 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGTCAt                                                                              >  1:4666589/1‑93 (MQ=255)

           ccGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAt                                                                              <  1:2809407/93‑1 (MQ=255)

             gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGc                                                                            <  1:4733966/93‑1 (MQ=255)

             gTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGc                                                                            <  1:4301453/93‑1 (MQ=255)

              tGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCa                                                                           <  1:2007249/93‑1 (MQ=255)

              tGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCa                                                                           <  1:2772992/93‑1 (MQ=255)

               gATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAg                                                                          >  1:3223670/1‑93 (MQ=255)

                aTGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGc                                                                         <  1:1784455/93‑1 (MQ=255)

                 tGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCg                                                                        <  2:3160487/93‑1 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       <  2:636229/93‑1 (MQ=255)

                  ggCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGt                                                                       <  2:2607910/93‑1 (MQ=255)

                   gCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTg                                                                      <  2:2969494/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    <  2:4648150/93‑1 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  2:2533298/1‑93 (MQ=255)

                     gTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGca                                                                    >  1:2446952/1‑93 (MQ=255)

                      tGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcac                                                                   >  1:3277432/1‑93 (MQ=255)

                      tGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcac                                                                   <  2:719888/93‑1 (MQ=255)

                       ggAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcaca                                                                  >  2:1204458/1‑93 (MQ=255)

                       ggAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGcaca                                                                  >  1:1160163/1‑93 (MQ=255)

                          aaaaGCAGTTACAGCGACAGGGAAGCGTATCCATCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               <  1:639140/93‑1 (MQ=255)

                          aaaaGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTAGTAAAAACCACTGAAGCCATGCAGCGTGCacaaca                                                               <  2:1717119/93‑1 (MQ=255)

                             aGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGc                                                            <  2:520692/93‑1 (MQ=255)

                              gCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCt                                                           >  2:1317218/1‑93 (MQ=255)

                                  ttACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaa                                                       >  1:3054583/1‑93 (MQ=255)

                                     cagcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaa                                                    <  2:136417/93‑1 (MQ=255)

                                      agcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaac                                                   <  1:442898/93‑1 (MQ=255)

                                      agcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTaacaac                                                   <  2:2331542/93‑1 (MQ=255)

                                       gcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  <  1:4456664/93‑1 (MQ=255)

                                       gcgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAAcc                                                  >  1:1383823/1‑93 (MQ=255)

                                        cgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACct                                                 <  2:4041794/93‑1 (MQ=255)

                                        cgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACct                                                 <  2:3603741/93‑1 (MQ=255)

                                        cgacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACct                                                 >  1:3557713/1‑93 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  1:1230580/93‑1 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctc                                                <  2:1072510/93‑1 (MQ=255)

                                         gacagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGCAACAACctc                                                >  2:1584235/1‑93 (MQ=255)

                                          acagcgCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACctct                                               <  2:832469/93‑1 (MQ=255)

                                           cagcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTa                                              <  2:1383565/93‑1 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            <  1:132349/93‑1 (MQ=255)

                                             gcgcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTAcc                                            >  1:188148/1‑93 (MQ=255)

                                               gcagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGt                                          <  1:1516675/93‑1 (MQ=255)

                                                cagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGtt                                         <  1:3046960/93‑1 (MQ=255)

                                                cagcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGtt                                         >  2:4512981/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  2:2094257/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  2:2759552/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  1:1962047/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       >  2:3958757/1‑93 (MQ=255)

                                                  gcgTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCAAGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  1:2584167/93‑1 (MQ=255)

                                                  gcgTATCGCTCAGTTAAAGCAGCAAATGAAAGTAGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTgg                                       <  1:4199345/93‑1 (MQ=255)

                                                         gCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCt                                >  2:1797120/1‑93 (MQ=255)

                                                          cTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCtt                               <  2:3846523/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  2:876030/1‑93 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              <  1:217779/93‑1 (MQ=255)

                                                           tCAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTc                              >  1:4500081/1‑93 (MQ=255)

                                                            cAGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcc                             <  2:2289059/93‑1 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            >  2:1465094/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            >  2:1642427/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            >  2:3597615/1‑93 (MQ=255)

                                                             aGTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCAAAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCa                            <  2:4396042/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTcccg                           >  2:350783/1‑91 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           >  2:1454435/1‑93 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           <  1:50986/93‑1 (MQ=255)

                                                              gTTTGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAg                           >  1:4601944/1‑93 (MQ=255)

                                                               tttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGc                          <  1:893473/93‑1 (MQ=255)

                                                                ttGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCg                         <  1:3786227/93‑1 (MQ=255)

                                                                 tGAGCAGCAAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGt                        <  2:1962268/93‑1 (MQ=255)

                                                                       gcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                  >  1:2080185/1‑93 (MQ=255)

                                                                       gcaAATGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGAc                  <  2:858328/93‑1 (MQ=255)

                                                                        caAATGGAAGTCGTTAAACCCCCTGAAGCCATGCAGCGTGCACAACAGTCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACa                 <  2:378597/93‑1 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                >  2:1162842/1‑93 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                <  1:3725190/93‑1 (MQ=255)

                                                                         aaaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAg                >  1:1080753/1‑93 (MQ=255)

                                                                          aaTGGAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGc               >  1:3374788/1‑93 (MQ=255)

                                                                             ggAAGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGc            <  2:1489938/93‑1 (MQ=255)

                                                                              gAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCa           <  2:3894614/93‑1 (MQ=255)

                                                                                aGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGa         <  2:2999708/93‑1 (MQ=255)

                                                                                aGTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGa         >  1:1022613/1‑93 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        <  2:3699411/93‑1 (MQ=255)

                                                                                 gTCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGaa        >  2:3447001/1‑93 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       >  1:3107702/1‑93 (MQ=255)

                                                                                  tCGTTAAAACCACTGAAGCCATCAAGCGTGCACAAAAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAAt       <  1:4059515/93‑1 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    <  2:3739834/93‑1 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    >  2:4136438/1‑93 (MQ=255)

                                                                                     ttAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCt    <  2:468445/93‑1 (MQ=255)

                                                                                      tAAAACCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCTACAGCGGCAGAATCCta   >  1:782412/1‑92 (MQ=255)

                                                                                       aaaaCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTa  <  1:1565416/93‑1 (MQ=255)

CTGGAAGTCTCCCGTGATGGCGTGGAAAAAGCAGTTACAGCGACAGCGCAGCGTATCGCTCAGTTTGAGCAGCAAATGGAAGTCGTTAAAGCCACTGAAGCCATGCAGCGTGCACAACAGGCTGTAACAACCTCTACCGTTGGCGCATCTTCCAGCGTTTCGACAGCGGCAGAATCCTTA  >  NC_000913/4408470‑4408649

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