BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	723,954	G→A	G83G (GGC→GGT)	kdpC ←	KdpC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	723,954	0	G	A	99.2%	403.1 / ‑5.6	120	G83G (GGC→GGT)	kdpC	KdpC
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (49/70); total (49/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 124 total aligned reads displayed.

TGGCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTAA  >  NC_000913/723866‑724043

tGGCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAg                                                                                       >  1:765961/1‑93 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     <  1:1478579/93‑1 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     <  2:3086274/93‑1 (MQ=255)

  gCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGaa                                                                                     >  1:3799482/1‑93 (MQ=255)

    aTCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGc                                                                                   <  1:3472956/93‑1 (MQ=255)

      ccGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTCGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCACCCACCAGAAGCCt                                                                                 <  1:1009905/93‑1 (MQ=255)

      ccGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCt                                                                                 <  2:1497237/93‑1 (MQ=255)

      ccGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCt                                                                                 <  1:3319123/93‑1 (MQ=255)

      ccGGGGTAGCGGCCCGTAATGCAGCATCGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTGACTGACCGCCAGATTGCTCCCGCCAGAAGCCt                                                                                 <  1:1009906/93‑1 (MQ=255)

       cGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtg                                                                                <  2:4067768/93‑1 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               >  2:1460143/1‑93 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               >  1:2250881/1‑93 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               <  2:2618145/93‑1 (MQ=255)

        gggTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCtgt                                                                               <  2:2042650/93‑1 (MQ=255)

          gTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTgg                                                                             >  2:4099740/1‑93 (MQ=255)

          gTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTgg                                                                             >  2:4099739/1‑93 (MQ=255)

            tAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAt                                                                           <  2:1928658/93‑1 (MQ=255)

            tAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAt                                                                           >  2:1792067/1‑93 (MQ=255)

            tAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAt                                                                           <  1:1622962/93‑1 (MQ=255)

            tAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAt                                                                           <  1:1622961/93‑1 (MQ=255)

             aGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGAtt                                                                          <  1:3235764/93‑1 (MQ=255)

                ggCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATtata                                                                       <  1:3850252/93‑1 (MQ=255)

                 gCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAg                                                                      <  1:1913286/93‑1 (MQ=255)

                  cgcGTAATGCAGCAACGCGTGCGGCTAGTAGTTTATCCGGCTCAGTGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAgg                                                                     <  1:2363909/91‑1 (MQ=255)

                  cccGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAgg                                                                     >  2:1157512/1‑93 (MQ=255)

                    cGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGc                                                                   >  1:4280959/1‑93 (MQ=255)

                     gTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCa                                                                  <  2:1393743/93‑1 (MQ=255)

                       aaTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAtt                                                                <  1:2246888/93‑1 (MQ=255)

                       aaTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAtt                                                                >  2:2994517/1‑93 (MQ=255)

                       aaTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAtt                                                                <  1:2246887/93‑1 (MQ=255)

                       aaTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAtt                                                                <  1:1743460/93‑1 (MQ=255)

                        aTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAttt                                                               <  1:2852373/93‑1 (MQ=255)

                        aTGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCAttt                                                               <  1:3326162/93‑1 (MQ=255)

                         tGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTc                                                              <  2:589180/93‑1 (MQ=255)

                         tGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTc                                                              <  1:4176776/93‑1 (MQ=255)

                         tGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTc                                                              >  2:520898/1‑93 (MQ=255)

                         tGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTc                                                              <  1:1632603/93‑1 (MQ=255)

                           cagcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTg                                                            <  1:547391/93‑1 (MQ=255)

                             gcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGcc                                                          <  2:1276308/93‑1 (MQ=255)

                             gcaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGcc                                                          >  2:1044084/1‑93 (MQ=255)

                             acaACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGcc                                                          <  1:786605/92‑1 (MQ=255)

                              caACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCg                                                         >  1:3446031/1‑93 (MQ=255)

                               aaCGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGt                                                        >  1:3955244/1‑93 (MQ=255)

                                  gcgTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGc                                                     >  2:3497912/1‑93 (MQ=255)

                                  gcgTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTACCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGc                                                     >  2:1342319/1‑93 (MQ=255)

                                    gTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCg                                                   >  2:3540537/1‑93 (MQ=255)

                                     tGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGa                                                  >  2:3969194/1‑93 (MQ=255)

                                      gCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGAc                                                 <  2:2337860/93‑1 (MQ=255)

                                      gCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGAc                                                 <  2:1799572/93‑1 (MQ=255)

                                       cGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACg                                                <  1:3067279/93‑1 (MQ=255)

                                        ggCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACgg                                               >  2:4277720/1‑93 (MQ=255)

                                        ggCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACgg                                               <  1:2439388/93‑1 (MQ=255)

                                            aTTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGa                                           >  2:2729164/1‑93 (MQ=255)

                                               tgtctATCCAGCTCAGTGTTAGTGACCGCTAGATTGCTCCGACCAGAAGCCTGTGGATTATAGGGTATTTCTGCCGTTGCCGCCGGGCGACCa                                        <  2:3994247/89‑1 (MQ=255)

                                               aGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCa                                        <  2:3073221/93‑1 (MQ=255)

                                                  ttATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGa                                     <  2:1496534/93‑1 (MQ=255)

                                                   tATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGaa                                    <  1:2505319/93‑1 (MQ=255)

                                                     tCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGaaaa                                  >  1:1833498/1‑93 (MQ=255)

                                                     tCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCAGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGaaaa                                  <  1:3558469/93‑1 (MQ=255)

                                                      ccAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAAt                                 <  1:2773398/93‑1 (MQ=255)

                                                      ccAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAAt                                 <  2:4024001/93‑1 (MQ=255)

                                                       catcTCAGGGTTACTGACCGCCAGATTTCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGGCGGGCGACCATGAAAATa                                >  2:1004861/4‑93 (MQ=255)

                                                       cAGCTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATa                                <  1:4729342/93‑1 (MQ=255)

                                                          cTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGcc                             <  2:3869594/93‑1 (MQ=255)

                                                          cTCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGcc                             >  1:1505735/1‑93 (MQ=255)

                                                           tCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCTACGGGCGACCATGGAAATAGCCg                            <  1:828988/93‑1 (MQ=255)

                                                           tCAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCg                            >  2:1321310/1‑93 (MQ=255)

                                                            cAGGTTCACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           <  1:3252401/93‑1 (MQ=255)

                                                            cAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           >  2:1026502/1‑93 (MQ=255)

                                                            cAGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGt                           <  1:2206644/93‑1 (MQ=255)

                                                             aGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtt                          >  1:4619566/1‑93 (MQ=255)

                                                             aGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtt                          >  2:460637/1‑93 (MQ=255)

                                                             aGGGTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtt                          <  1:4570364/93‑1 (MQ=255)

                                                              gggTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGtgg                         >  1:3681847/1‑91 (MQ=255)

                                                               ggTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGTTTATGGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGc                        <  2:3617565/93‑1 (MQ=255)

                                                                gTTACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGcc                       >  2:4363410/1‑93 (MQ=255)

                                                                 ttACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCg                      >  1:3905026/1‑93 (MQ=255)

                                                                 ttACTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCg                      <  1:1119836/93‑1 (MQ=255)

                                                                   aCTGACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGt                    <  1:361116/93‑1 (MQ=255)

                                                                      gACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTaaa                 >  1:1646035/1‑93 (MQ=255)

                                                                      gACCGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCAGGTaaa                 >  1:2978408/1‑93 (MQ=255)

                                                                        ccGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAAt               <  2:4352537/93‑1 (MQ=255)

                                                                        ccGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAAt               >  2:72248/1‑93 (MQ=255)

                                                                         cGCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAAtt              >  1:1586084/1‑93 (MQ=255)

                                                                          gCCAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTc             <  1:2397647/93‑1 (MQ=255)

                                                                            cAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAACATTCtc           >  1:656370/1‑92 (MQ=255)

                                                                            cAGATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTg           <  2:134532/93‑1 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         <  1:4176855/93‑1 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  2:4333022/1‑93 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGcc         >  1:2227188/1‑93 (MQ=255)

                                                                              gATTGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCGGcc         >  1:2922537/1‑93 (MQ=255)

                                                                                ttGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCg       >  1:1005228/1‑93 (MQ=255)

                                                                                ttGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCg       >  1:1005229/1‑93 (MQ=255)

                                                                                 tGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGa      <  2:3192902/93‑1 (MQ=255)

                                                                                 tGCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGa      >  1:1233009/1‑93 (MQ=255)

                                                                                  gCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAt     <  2:3391800/93‑1 (MQ=255)

                                                                                  gCTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAt     <  2:3442559/93‑1 (MQ=255)

                                                                                   cTCCCACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGAtt    <  2:4252161/93‑1 (MQ=255)

                                                                                     cccACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTaa  <  1:4702210/93‑1 (MQ=255)

                                                                                     cccACCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTaa  <  2:2683652/93‑1 (MQ=255)

TGGCATCCGGGTTAGCGGCCCGTAATGCAGCAACGCGTGCGGCTATTAGTTTATCCAGCTCAGGGTTACTGACCGCCAGATTGCTCCCGCCAGAAGCCTGTGGATTATAGGGCATTTCTGCCGTTGCCGACGGGCGACCATGAAAATAGCCGTTGCCGGTAAAATTCTGCCCGATTAA  >  NC_000913/723866‑724043

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