BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,293,354	A→G	H202R (CAT→CGT)	tdk →	thymidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,293,354	0	A	G	98.4%	575.9 / ‑0.9	182	H202R (CAT→CGT)	tdk	thymidine kinase
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (89/90); total (92/90)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 187 total aligned reads displayed.

ATTGGTGGTAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCATCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCGGGCTA  >  NC_000913/1293263‑1293443

 ttGGTGGTAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTc                                                                                         >  2:4707420/1‑93 (MQ=255)

 ttGGTGGTAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTc                                                                                         <  1:1951516/93‑1 (MQ=255)

   ggtggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGc                                                                                       <  1:4479758/93‑1 (MQ=255)

    gtggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGcc                                                                                      >  1:2570881/1‑93 (MQ=255)

     tggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCa                                                                                     <  2:4470270/93‑1 (MQ=255)

     tggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCa                                                                                     <  2:1227346/93‑1 (MQ=255)

     tggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCa                                                                                     <  1:2931125/93‑1 (MQ=255)

      ggtAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCAc                                                                                    <  2:3604179/93‑1 (MQ=255)

        tAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGa                                                                                  >  1:1538489/1‑93 (MQ=255)

         aaTGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGAt                                                                                 <  1:4643604/93‑1 (MQ=255)

         aaTGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTCACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGAt                                                                                 >  1:2357758/1‑93 (MQ=255)

          aTGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGAtt                                                                                >  2:3652638/1‑93 (MQ=255)

           tGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATta                                                                               <  2:1392009/93‑1 (MQ=255)

           agAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATta                                                                               <  1:3861578/92‑1 (MQ=255)

            gAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATtaa                                                                              <  1:1012217/93‑1 (MQ=255)

              aCGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCATGCAACGCGTCGACACGATtaata                                                                            >  2:1781350/1‑93 (MQ=255)

              aCGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATtaata                                                                            >  2:1474175/1‑93 (MQ=255)

              aCGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATtaata                                                                            <  1:3290045/93‑1 (MQ=255)

                gATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAg                                                                          >  1:2682259/1‑93 (MQ=255)

                gATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAg                                                                          >  1:3575781/1‑93 (MQ=255)

                 aTACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGa                                                                         >  1:3664538/1‑93 (MQ=255)

                 aTACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCATCGCCACGATTAATAAGa                                                                         >  1:4013817/1‑93 (MQ=255)

                   aCGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAAt                                                                       <  1:1573584/93‑1 (MQ=255)

                   aCGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAAt                                                                       <  1:3545610/93‑1 (MQ=255)

                     gTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAAttt                                                                     <  1:2212458/93‑1 (MQ=255)

                     gTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAAttt                                                                     <  2:608077/93‑1 (MQ=255)

                     gTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAAttt                                                                     >  1:3240124/1‑93 (MQ=255)

                      tttCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTc                                                                    <  2:1729627/93‑1 (MQ=255)

                        tCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCtt                                                                  <  1:2028292/93‑1 (MQ=255)

                        tCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCtt                                                                  <  1:942669/93‑1 (MQ=255)

                          tGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAGCGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTa                                                                >  2:2810168/1‑93 (MQ=255)

                          tGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTa                                                                <  2:2293289/93‑1 (MQ=255)

                          tGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTa                                                                <  1:4745369/93‑1 (MQ=255)

                            tATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACt                                                              >  1:3405398/1‑93 (MQ=255)

                            tATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACt                                                              <  2:4616743/93‑1 (MQ=255)

                            tATGCCGTAAACACTATAAAGAGCCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACt                                                              <  2:2605525/93‑1 (MQ=255)

                             aTGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTg                                                             >  2:289203/1‑93 (MQ=255)

                              tGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGa                                                            >  2:3068735/1‑93 (MQ=255)

                               gCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGAc                                                           >  2:2019576/1‑93 (MQ=255)

                               gCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGAc                                                           <  2:3023750/93‑1 (MQ=255)

                                 cGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAg                                                         >  1:738533/1‑93 (MQ=255)

                                    aaaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGt                                                      <  2:679486/93‑1 (MQ=255)

                                    aaaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGt                                                      <  2:4259187/93‑1 (MQ=255)

                                     aaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTg                                                     <  2:1620279/93‑1 (MQ=255)

                                     aaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTg                                                     >  1:541217/1‑93 (MQ=255)

                                     aaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTg                                                     <  2:2790044/93‑1 (MQ=255)

                                     aaCACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTg                                                     <  1:2991691/93‑1 (MQ=255)

                                      acacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGa                                                    <  1:4446043/93‑1 (MQ=255)

                                      acacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGa                                                    <  2:4066650/93‑1 (MQ=255)

                                      acacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGa                                                    >  1:3781734/1‑93 (MQ=255)

                                      acacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGa                                                    <  2:1031032/93‑1 (MQ=255)

                                       cacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAg                                                   >  1:1103919/1‑93 (MQ=255)

                                       cacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAg                                                   >  2:3792437/1‑93 (MQ=255)

                                       cacTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAg                                                   <  1:984959/93‑1 (MQ=255)

                                        acTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGc                                                  <  2:2399903/93‑1 (MQ=255)

                                          tataAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAg                                                <  2:3088715/93‑1 (MQ=255)

                                          tataAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAg                                                <  1:3280941/93‑1 (MQ=255)

                                           ataAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAgg                                               >  1:1520611/1‑93 (MQ=255)

                                           ataAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAgg                                               >  1:3658026/1‑93 (MQ=255)

                                            taAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAggg                                              >  2:885442/1‑93 (MQ=255)

                                            taAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAggg                                              >  2:3948683/1‑93 (MQ=255)

                                              aaGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCa                                            >  1:800880/1‑93 (MQ=255)

                                               agagGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCAc                                           <  2:2136655/93‑1 (MQ=255)

                                                gagGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACt                                          <  2:1271688/93‑1 (MQ=255)

                                                gagGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACt                                          >  2:3415433/1‑93 (MQ=255)

                                                  ggCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACttt                                        <  1:2174126/93‑1 (MQ=255)

                                                   gCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACtttt                                       >  1:4106498/1‑93 (MQ=255)

                                                        aCAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCt                                  >  2:4190394/1‑93 (MQ=255)

                                                         cAAGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTg                                 >  1:2845604/1‑93 (MQ=255)

                                                          aaGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGt                                >  1:1404423/1‑93 (MQ=255)

                                                           aGTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTc                               >  1:2989256/1‑93 (MQ=255)

                                                            gTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCa                              <  2:520594/93‑1 (MQ=255)

                                                            gTCGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCa                              <  2:3902942/93‑1 (MQ=255)

                                                              cGACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGt                            <  2:3338962/93‑1 (MQ=255)

                                                               gACTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGtt                           >  1:2642547/1‑93 (MQ=255)

                                                                aCTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTc                          <  1:8366/93‑1 (MQ=255)

                                                                aCTCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTc                          <  1:3953722/93‑1 (MQ=255)

                                                                  tCATTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGt                        >  1:688973/1‑93 (MQ=255)

                                                                    aTTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGttt                      >  1:1012146/1‑93 (MQ=255)

                                                                    aTTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGttt                      >  2:505969/1‑93 (MQ=255)

                                                                    aTTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGttt                      >  1:4020449/1‑93 (MQ=255)

                                                                    aTTAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGttt                      >  2:1173539/1‑93 (MQ=255)

                                                                     ttAACGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGtttt                     >  1:1261185/1‑93 (MQ=255)

                                                                       aaCGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTAc                   >  1:3921088/1‑93 (MQ=255)

                                                                        aCGGCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACg                  >  1:2424774/1‑93 (MQ=255)

                                                                          ggCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCa                <  2:1121217/93‑1 (MQ=255)

                                                                          ggCTATTCAGGAAAGCCGTCGACACTGTTAATAAGAAGTTCTTTACGGACAGGGTGAGCAGGGCACTGTTAGCCTGTCAGTTCGTTTTACGCa                <  1:4672691/93‑1 (MQ=255)

                                                                           gCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCAc               >  1:3333893/1‑93 (MQ=255)

                                                                           gCTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCAc               <  1:168229/93‑1 (MQ=255)

                                                                            cTATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCAct              >  2:1937610/1‑93 (MQ=255)

                                                                             tATTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCActt             >  1:1032048/1‑93 (MQ=255)

                                                                              aTTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCActtc            >  1:2761585/1‑93 (MQ=255)

                                                                              aTTCAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCActtc            >  1:1433465/1‑93 (MQ=255)

                                                                                 cAGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACttcttc         <  2:3869473/93‑1 (MQ=255)

                                                                                  aGGAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTcc        >  1:4541568/1‑93 (MQ=255)

                                                                                   ggAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCg       >  2:3108390/1‑93 (MQ=255)

                                                                                    gAAAGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCgg      >  2:3526450/1‑93 (MQ=255)

                                                                                     aaaGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCggg     <  2:2598766/93‑1 (MQ=255)

                                                                                       aGGCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCGGGCt   >  2:4366870/1‑93 (MQ=255)

                                                                                        ggCGTCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCGGGCta  <  1:4435152/93‑1 (MQ=255)

ATTGGTGGTAATGAACGATACGTTTCTGTATGCCGTAAACACTATAAAGAGGCGTTACAAGTCGACTCATTAACGGCTATTCAGGAAAGGCATCGCCACGATTAATAAGAATTTCTTTACTGACAGGGTGAGCAGGGCACTTTTATCCTGTCAGTTCGTTTTACGCACTTCTTCCGGGCTA  >  NC_000913/1293263‑1293443

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