BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,522,123	A→G	R152R (CGT→CGC)	b4298 ←	putative lyase/synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,522,123	0	A	G	84.3%	699.4 / 120.0	286	R152R (CGT→CGC)	b4298	putative lyase/synthase
Reads supporting (aligned to +/- strand): ref base A (24/21); new base G (131/110); total (155/131)
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.62e-01

Only 100 of 293 total aligned reads displayed.

CCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCT  >  NC_000913/4522032‑4522214

ccGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGAc                                                                                            <  1:5942119/93‑1 (MQ=255)

 cGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGTGCCAGGcg                                                                                           <  1:6478900/93‑1 (MQ=255)

 cGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGcg                                                                                           >  1:3344889/1‑93 (MQ=255)

 cGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGcg                                                                                           <  1:5599787/93‑1 (MQ=255)

 cGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGcg                                                                                           >  1:2145072/1‑93 (MQ=255)

   aCTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGcgcg                                                                                         >  2:4538231/1‑93 (MQ=255)

     ttAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTc                                                                                       >  2:768095/1‑93 (MQ=255)

      tAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCa                                                                                      >  2:6247575/1‑93 (MQ=255)

       aaCTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCAc                                                                                     <  2:4028070/93‑1 (MQ=255)

       aaCTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCAc                                                                                     <  2:2220234/93‑1 (MQ=255)

          tgtgTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGt                                                                                  <  2:6460483/93‑1 (MQ=255)

          tgtgTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGt                                                                                  >  1:220542/1‑93 (MQ=255)

          tgtgTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGt                                                                                  >  1:1480319/1‑93 (MQ=255)

           gtgtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGtt                                                                                 <  1:1129473/93‑1 (MQ=255)

           gtgtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGtt                                                                                 <  1:627276/93‑1 (MQ=255)

            tgtTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGttt                                                                                >  2:2585975/1‑93 (MQ=255)

              ttGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTcc                                                                              <  2:108519/93‑1 (MQ=255)

                 atcatcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCggg                                                                           <  1:5027057/93‑1 (MQ=255)

                  tcatcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCgggg                                                                          >  1:3531569/1‑93 (MQ=255)

                     tcGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTa                                                                       >  2:3030243/1‑93 (MQ=255)

                      cGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTaa                                                                      <  1:5151174/93‑1 (MQ=255)

                       gTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAg                                                                     >  2:2328120/1‑93 (MQ=255)

                       gTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAg                                                                     >  2:850569/1‑93 (MQ=255)

                        tACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGt                                                                    >  1:2646827/1‑93 (MQ=255)

                        tACGCAAGTGACCAACGCTGTCGATGGTGGCTTTGATGCCAACGATATTCTCGTTTTGCAGAGACAGGCGCGTCACGGTTTCCGGGGTTAAGt                                                                    >  2:3767144/1‑93 (MQ=255)

                         aCGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTc                                                                   <  1:657200/93‑1 (MQ=255)

                          cGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTcc                                                                  <  1:2969122/93‑1 (MQ=255)

                           gCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCt                                                                 <  2:1444556/93‑1 (MQ=255)

                           gCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCt                                                                 <  2:2127312/93‑1 (MQ=255)

                            cAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTg                                                                <  1:1391273/93‑1 (MQ=255)

                              aGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGAc                                                              >  1:2072/1‑93 (MQ=255)

                               gTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGAcc                                                             <  2:5170624/93‑1 (MQ=255)

                               gTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGAcc                                                             <  2:5621490/93‑1 (MQ=255)

                                   ccAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTc                                                         <  1:955478/93‑1 (MQ=255)

                                   ccAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTc                                                         >  1:2219696/1‑93 (MQ=255)

                                   ccAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTc                                                         <  2:2305347/93‑1 (MQ=255)

                                   ccAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTc                                                         >  1:3142244/1‑93 (MQ=255)

                                   caaaCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTc                                                         >  1:6139137/3‑93 (MQ=255)

                                     aaCGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAg                                                       >  2:3487135/1‑93 (MQ=255)

                                         cTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATcc                                                   >  1:177/1‑93 (MQ=255)

                                         cTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATcc                                                   >  1:3200884/1‑93 (MQ=255)

                                            tCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGa                                                >  1:5395954/1‑93 (MQ=255)

                                            tCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGa                                                >  2:2037105/1‑93 (MQ=255)

                                              gATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGaaa                                              >  2:5973640/1‑93 (MQ=255)

                                               aTGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTGCTGACCCGTCAGATCCGGAAAg                                             >  1:2365870/1‑93 (MQ=255)

                                               aTGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAg                                             >  2:4770853/1‑93 (MQ=255)

                                               aTGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGATTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAg                                             >  2:3795813/1‑93 (MQ=255)

                                                tGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGt                                            >  2:1184974/1‑93 (MQ=255)

                                                tGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGt                                            <  1:1066489/93‑1 (MQ=255)

                                                  gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTg                                          >  2:5541627/1‑93 (MQ=255)

                                                  gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTg                                          >  1:3655612/1‑93 (MQ=255)

                                                  gtgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTg                                          <  2:3600643/93‑1 (MQ=255)

                                                   tgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGt                                         >  2:530770/1‑93 (MQ=255)

                                                   tgtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGt                                         <  1:5083036/93‑1 (MQ=255)

                                                    gtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTa                                        >  2:91880/1‑93 (MQ=255)

                                                    gtCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTa                                        <  1:5515247/93‑1 (MQ=255)

                                                     tCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTAc                                       >  2:5776208/1‑93 (MQ=255)

                                                      cTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACa                                      <  1:2037105/93‑1 (MQ=255)

                                                       tttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAg                                     >  2:1705511/1‑93 (MQ=255)

                                                       tttGATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAg                                     <  2:4366600/93‑1 (MQ=255)

                                                          gATGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGAt                                  >  2:5915921/1‑93 (MQ=255)

                                                           aTGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATc                                 >  1:5300408/1‑93 (MQ=255)

                                                            tGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCa                                >  1:4025037/1‑93 (MQ=255)

                                                            tGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCa                                >  1:134095/1‑93 (MQ=255)

                                                            tGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCa                                <  2:3164036/93‑1 (MQ=255)

                                                            tGCCAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCa                                <  1:4483397/93‑1 (MQ=255)

                                                              ccAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCAcc                              >  2:1658143/1‑93 (MQ=255)

                                                               cAACGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCg                             >  2:5013098/1‑93 (MQ=255)

                                                                aaCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCgg                            >  1:1701375/1‑93 (MQ=255)

                                                                aaCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCgg                            >  2:5352513/1‑93 (MQ=255)

                                                                aaCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCgg                            <  1:4245447/93‑1 (MQ=255)

                                                                 aCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGt                           >  1:4413533/1‑93 (MQ=255)

                                                                 aCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGt                           <  1:6114018/93‑1 (MQ=255)

                                                                 aCGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGt                           <  2:1975269/93‑1 (MQ=255)

                                                                  cGATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTa                          <  1:6376121/93‑1 (MQ=255)

                                                                   gATATTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAg                         <  2:1215187/93‑1 (MQ=255)

                                                                     tatTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAggg                       <  2:4949805/93‑1 (MQ=255)

                                                                     tatTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAggg                       <  1:1022975/93‑1 (MQ=255)

                                                                     tatTCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAggg                       <  1:4284773/93‑1 (MQ=255)

                                                                       ttCTCGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTg                     >  1:3975451/1‑93 (MQ=255)

                                                                          tcGTTTTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACg                  <  2:464560/93‑1 (MQ=255)

                                                                            gTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCt                >  2:5462811/1‑93 (MQ=255)

                                                                            gTTGTGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCt                <  1:2500215/93‑1 (MQ=255)

                                                                             ttttGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTa               >  2:1551667/1‑93 (MQ=255)

                                                                             ttttGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTa               >  2:5226847/1‑93 (MQ=255)

                                                                              tttGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGGGACGCTAc              >  2:4585479/1‑93 (MQ=255)

                                                                                tGCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACgg            <  2:3308138/93‑1 (MQ=255)

                                                                                tGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACgg            <  1:4891286/93‑1 (MQ=255)

                                                                                 gCAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACggg           <  1:5875579/93‑1 (MQ=255)

                                                                                 gCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACggg           <  1:4617451/93‑1 (MQ=255)

                                                                                  cAGAGCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGc          >  2:4159304/1‑93 (MQ=255)

                                                                                   agagCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCg         >  2:2515872/1‑93 (MQ=255)

                                                                                   agagCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCg         >  2:4846794/1‑93 (MQ=255)

                                                                                    gagCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGa        >  2:3361857/1‑93 (MQ=255)

                                                                                      gCCAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATc      <  1:277847/93‑1 (MQ=255)

                                                                                       ccAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATct     <  2:2840731/93‑1 (MQ=255)

                                                                                       ccAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATct     <  1:2483018/93‑1 (MQ=255)

                                                                                       ccAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATct     <  1:722228/93‑1 (MQ=255)

                                                                                        cAGGCGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctg    >  2:3697867/1‑93 (MQ=255)

                                                                                          gACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATctgct  <  1:3286335/93‑1 (MQ=255)

CCGACTTAACTGTGTTGATCATCGTACGCAAGTGACCAACGCTGTCGATGGTGTCTTTGATGCCAACGATATTCTCGTTTTGCAGAGCCAGACGCGTCACGGTTTCCGGGGTTAAGTCCTGACCCGTCAGATCCGGAAAGTTGTACAGGATCACCGGTAGGGTGACGCTACGGGCGATCTGCT  >  NC_000913/4522032‑4522214

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