BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,506,318	(C)_5→6	coding (1134/1227 nt)	yhfW ←	putative mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,506,313	1	.	C	98.2%	754.1 / 4.1	222	coding (1139/1227 nt)	yhfW	putative mutase
Reads supporting (aligned to +/- strand): ref base . (2/2); new base C (117/101); total (119/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 259 total aligned reads displayed.

AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGG  >  NC_000913/3506221‑3506399

 aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑c                                                                                       <  2:4361453/93‑1 (MQ=255)

 aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑c                                                                                       >  1:6154283/1‑93 (MQ=255)

  aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑cc                                                                                      <  2:2415020/93‑1 (MQ=255)

  aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTCGCGCGCGGAAAAATTCACACACGGT‑AG‑cc                                                                                      >  1:1396489/1‑93 (MQ=255)

   cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑cac                                                                                     >  1:5211482/1‑91 (MQ=255)

    tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑cccc                                                                                    >  1:5714830/1‑93 (MQ=255)

    tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑cccc                                                                                    >  1:4153554/1‑93 (MQ=255)

     cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑ccccc                                                                                   <  1:1598550/93‑1 (MQ=255)

     cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑ccccc                                                                                   <  1:1776015/93‑1 (MQ=255)

         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑ccccccaca                                                                               <  1:3025986/93‑5 (MQ=255)

         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑ccccccaca                                                                               >  2:196103/1‑89 (MQ=255)

         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑ccccccaca                                                                               >  2:5525476/1‑89 (MQ=255)

            tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACGATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCa                                                                             >  2:427242/1‑93 (MQ=255)

             gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAg                                                                            >  2:2725670/1‑93 (MQ=255)

              tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGa                                                                           <  1:754922/93‑1 (MQ=255)

                tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGaaa                                                                         >  1:2778521/1‑93 (MQ=255)

                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAg                                                                        >  1:2073828/1‑93 (MQ=255)

                 ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAg                                                                        >  2:1342751/1‑93 (MQ=255)

                  cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGc                                                                       >  1:5875923/1‑93 (MQ=255)

                   tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCg                                                                      >  2:2898357/1‑93 (MQ=255)

                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgt                                                                     <  1:3775788/93‑1 (MQ=255)

                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgt                                                                     <  1:2391188/93‑1 (MQ=255)

                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtg                                                                    <  2:690988/93‑1 (MQ=255)

                     cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtg                                                                    >  2:3588778/1‑93 (MQ=255)

                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtgg                                                                   <  1:3605916/93‑1 (MQ=255)

                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtgg                                                                   <  2:2380271/93‑1 (MQ=255)

                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtgg                                                                   >  1:2511363/1‑93 (MQ=255)

                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggt                                                                  <  1:1074107/93‑1 (MQ=255)

                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggt                                                                  <  1:4600652/93‑1 (MQ=255)

                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggt                                                                  >  1:4520238/1‑93 (MQ=255)

                        aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggtg                                                                 <  2:3317214/93‑1 (MQ=255)

                        aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggtg                                                                 >  2:3200084/1‑93 (MQ=255)

                        aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCgtggtg                                                                 >  1:5158309/1‑93 (MQ=255)

                         ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTgc                                                                >  2:3381375/1‑93 (MQ=255)

                         ccGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTgc                                                                >  1:986726/1‑93 (MQ=255)

                           ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAAC‑CCCCCACATCAGAAAGCGTGGTgcgc                                                              >  1:4003635/1‑93 (MQ=255)

                            gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGca                                                             >  2:4507280/1‑93 (MQ=255)

                            gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGGGGTGCGca                                                             >  1:3476131/1‑93 (MQ=255)

                             agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGcac                                                            <  1:2326147/93‑1 (MQ=255)

                              gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGcaca                                                           >  2:2582362/1‑93 (MQ=255)

                              gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGcaca                                                           >  1:6276114/1‑93 (MQ=255)

                               ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCacac                                                          >  1:3492697/1‑93 (MQ=255)

                               ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCacac                                                          <  1:6139223/93‑1 (MQ=255)

                                gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                         <  2:487098/93‑1 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                       >  2:3103973/1‑93 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                       >  1:3266757/1‑93 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                       <  1:1134941/93‑1 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                       >  1:4775535/1‑93 (MQ=255)

                                   gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                      >  1:1739692/1‑93 (MQ=255)

                                    aaaaaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                     <  1:5693079/93‑1 (MQ=255)

                                     aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCt                                                    >  1:1364497/1‑93 (MQ=255)

                                       gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                  >  1:386020/1‑93 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGGAAGCGTGGTGCGCACACCGAGCTGCGTa                                               >  1:2536817/1‑93 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                               >  1:4315583/1‑93 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                               >  1:895245/1‑93 (MQ=255)

                                           agagCGACCATTTTGTGGCGGTGGCGCGCAGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                              <  2:37229/93‑1 (MQ=255)

                                              gCGACCATTTTGTGGCGGTGGCGCGCTGAAAAATTCACACACCGT‑AGCCCCCCCCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                           <  2:3621160/93‑1 (MQ=255)

                                              gCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                           <  2:1621057/93‑1 (MQ=255)

                                                 aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                        >  2:3263698/1‑93 (MQ=255)

                                                    aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                     >  2:2168129/1‑93 (MQ=255)

                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                    >  2:4142980/1‑93 (MQ=255)

                                                       ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                  <  2:6279964/93‑1 (MQ=255)

                                                        tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                 >  1:2652098/1‑93 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                <  2:5168012/93‑1 (MQ=255)

                                                         gtggccgtagcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                <  2:5671971/93‑1 (MQ=255)

                                                          tggcggtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                               <  2:6013464/93‑1 (MQ=255)

                                                          tggcggtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                               <  1:1034312/93‑1 (MQ=255)

                                                          tggcggtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                               <  1:22530/93‑1 (MQ=255)

                                                             cggtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGAt                            <  2:3960574/93‑1 (MQ=255)

                                                               gtggcgCGCGGAAAAATTCACACACGGTAAC‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                          >  1:3123988/1‑93 (MQ=255)

                                                               gtggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                          >  1:5641059/1‑93 (MQ=255)

                                                                 ggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                        >  2:3910160/1‑93 (MQ=255)

                                                                 ggcgCGCGGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                        <  2:118079/93‑1 (MQ=255)

                                                                   cgcgcgGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                      <  1:298051/93‑1 (MQ=255)

                                                                    gcgcgGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                     >  2:90230/1‑93 (MQ=255)

                                                                     cgcgGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                    >  2:6135806/1‑93 (MQ=255)

                                                                      gcgGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                   <  1:6383629/93‑1 (MQ=255)

                                                                       cgGAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                  <  2:5352538/93‑1 (MQ=255)

                                                                         gAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                >  2:5010704/1‑93 (MQ=255)

                                                                         gAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                <  1:5115297/93‑1 (MQ=255)

                                                                         gAAAAATTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                <  2:4751888/93‑1 (MQ=255)

                                                                          aaaaaTTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg               <  2:2535739/93‑1 (MQ=255)

                                                                           aaaaTTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc              >  1:548224/1‑93 (MQ=255)

                                                                           aaaaTTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc              <  2:4700066/93‑1 (MQ=255)

                                                                            aaaTTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGca             >  1:1896836/1‑93 (MQ=255)

                                                                             aaTTCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac            <  1:3528863/93‑1 (MQ=255)

                                                                               ttCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca          >  1:5100822/1‑93 (MQ=255)

                                                                               ttCACACACGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca          >  2:953559/1‑93 (MQ=255)

                                                                                  ccacacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt       <  2:1281788/92‑1 (MQ=255)

                                                                                  acacacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt       <  2:1596131/93‑1 (MQ=255)

                                                                                    acacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg     >  2:4398499/1‑93 (MQ=255)

                                                                                    acacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg     >  1:5389788/1‑93 (MQ=255)

                                                                                     cacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc    <  2:1836037/93‑1 (MQ=255)

                                                                                     cacGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgc    >  1:5770979/1‑93 (MQ=255)

                                                                                      acGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAATCAGCACTGGCACCACTTcgcg   <  1:4089285/93‑1 (MQ=255)

                                                                                      acGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg   >  1:1437973/1‑93 (MQ=255)

                                                                                      acGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg   >  2:5419496/1‑93 (MQ=255)

                                                                                       cGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGNACNCCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg  <  1:4903425/93‑1 (MQ=255)

                                                                                       cGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg  <  2:2543300/93‑1 (MQ=255)

                                                                                       cGGT‑AGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg  >  1:4257094/1‑93 (MQ=255)

AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGT‑AG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGG  >  NC_000913/3506221‑3506399

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