BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,838,222	A→T	15.3%	intergenic (‑117/+32)	hydN ← / ← ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,222	0	A	T	15.3%	359.2 / 53.7	144	intergenic (‑117/+32)	hydN/ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand): ref base A (62/60); new base T (9/13); total (71/73)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 152 total aligned reads displayed.

ACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGA  >  NC_000913/2838123‑2838321

aCAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGa                                                                                                     >  2:1411639/1‑100 (MQ=255)

 cAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGac                                                                                                    <  2:3608687/100‑1 (MQ=255)

 cAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGac                                                                                                    <  2:3288587/100‑1 (MQ=255)

 cAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGac                                                                                                    >  2:2831205/1‑100 (MQ=255)

 cAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGac                                                                                                    >  2:2829950/1‑100 (MQ=255)

    ttAATGCGCCGATACATTCACAAATCATTCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACt                                                                                                 <  1:1789533/100‑1 (MQ=255)

     tAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTa                                                                                                >  2:2726462/1‑100 (MQ=255)

     tAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTa                                                                                                >  2:1585238/1‑100 (MQ=255)

          cgcCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGAc                                                                                           >  1:1612603/1‑100 (MQ=255)

           gcCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACg                                                                                          >  2:2358686/1‑100 (MQ=255)

            tcGATACATGCACAAATCATGCCAGTTTTTAATTTACTGTTATTTTAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGAc                                                                                           <  1:1982976/97‑1 (MQ=255)

              gATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATg                                                                                       <  1:2732808/100‑1 (MQ=255)

               aTACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGa                                                                                      <  2:3576132/100‑1 (MQ=255)

               aTACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGa                                                                                      <  2:1593400/100‑1 (MQ=255)

               aTACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGa                                                                                      <  2:1728342/100‑1 (MQ=255)

                tACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGAc                                                                                     <  2:2173669/100‑1 (MQ=255)

                  cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAt                                                                                   <  1:1290209/100‑1 (MQ=255)

                  cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAt                                                                                   >  2:2360933/1‑100 (MQ=255)

                  cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAt                                                                                   >  1:3723866/1‑100 (MQ=255)

                   aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtc                                                                                  <  2:486098/100‑2 (MQ=255)

                   aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtg                                                                                  >  2:3208724/1‑100 (MQ=255)

                   aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGAAGTCATCGACACTAGTGACAATGACAtg                                                                                  >  1:1389292/1‑100 (MQ=255)

                   aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGAATGAAAAACGATGTCATCGACACTAGTGACGATGACAtg                                                                                  >  1:3415537/1‑100 (MQ=255)

                    ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                 <  1:2123948/100‑1 (MQ=255)

                     tCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                <  1:2089729/100‑1 (MQ=255)

                      cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                               >  2:2279296/1‑100 (MQ=255)

                      cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                               <  1:2418887/100‑1 (MQ=255)

                      cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                               >  1:182948/1‑100 (MQ=255)

                       accaaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGCAAAAACGATCTCATCGACCCTAGTGACGATGCCATGtga                                                                               >  1:2997426/4‑99 (MQ=255)

                          aaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGAc                                                                           >  1:114769/1‑100 (MQ=255)

                            tCATGCACGTTATTAATTTTCTGTTATATGATGAAATTAATTTCTGTAATGCAGGAAAAACGTTGTAATCGTCACTAGTGTCGCTGCCAtcgtttttc                                                                           <  1:1710768/98‑9 (MQ=16)

                              ntGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTATTTTAGGTAAAACGATGTAATCTACACTAGTGACGCTGACATATGATGACAAt                                                                        >  1:3739408/2‑99 (MQ=255)

                              aTGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATg                                                                       >  2:198104/1‑100 (MQ=255)

                              aTGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATg                                                                       <  2:1334574/100‑1 (MQ=255)

                               tGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                      >  2:2363082/1‑100 (MQ=255)

                                gCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGtt                                                                     >  2:2988460/1‑100 (MQ=255)

                                gCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGtt                                                                     <  2:668407/100‑1 (MQ=255)

                                 ccAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                    >  1:1299413/1‑100 (MQ=255)

                                 ccAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTCACGATGACATGTGATGACAATGt                                                                      >  1:2933424/1‑98 (MQ=255)

                                  cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                   <  2:445832/100‑1 (MQ=255)

                                   aGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGtt                                                                     <  1:587765/97‑1 (MQ=255)

                                    gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATgagatcg                                                                 >  1:122481/1‑94 (MQ=255)

                                    gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATc                                                                 <  2:594552/100‑1 (MQ=255)

                                    gTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATc                                                                 <  1:2643812/100‑1 (MQ=255)

                                    gTTTTTAATTTACTGTTATTGAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAt                                                                  >  1:502534/1‑99 (MQ=255)

                                     tttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcg                                                                <  1:1541293/100‑1 (MQ=255)

                                      ttttAATTTACTGTTATTTAATGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgc                                                               <  1:3434077/100‑1 (MQ=255)

                                      ttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgc                                                               >  1:3049287/1‑100 (MQ=255)

                                       tttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgcg                                                              <  1:3452100/100‑1 (MQ=255)

                                        ttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGa                                                             <  1:1541868/100‑1 (MQ=255)

                                           aTTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAgg                                                          <  2:3507413/100‑1 (MQ=255)

                                           aTTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAgg                                                          <  1:2886424/100‑1 (MQ=255)

                                             ttACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAg                                                        <  1:1428967/100‑1 (MQ=255)

                                              tACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGc                                                       <  2:1861947/100‑1 (MQ=255)

                                               aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGCCAATGTTTATCGCGAAGGAGCa                                                      <  1:3787653/100‑1 (MQ=255)

                                               aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCa                                                      <  2:3727074/100‑1 (MQ=255)

                                                cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                     >  1:2541511/1‑100 (MQ=255)

                                                  gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtg                                                   <  1:3032015/100‑1 (MQ=255)

                                                  gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtg                                                   <  2:182952/100‑1 (MQ=255)

                                                  gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtg                                                   <  1:2402882/100‑1 (MQ=255)

                                                  gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTAATCGACACTAGTGACGATGACATGTGATGACAGTGGTTATCGCGAAGGATCAAtg                                                   >  1:1670728/1‑100 (MQ=255)

                                                        tAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTAATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                             >  1:782653/1‑100 (MQ=255)

                                                        naaGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgag                                             >  1:3532944/2‑100 (MQ=255)

                                                           ggAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTcg                                          >  2:3640500/1‑100 (MQ=255)

                                                            gAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTcgc                                         >  1:2891794/1‑100 (MQ=255)

                                                                ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcggcg                                     >  2:1372613/1‑100 (MQ=255)

                                                                ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCgcggcg                                     <  1:3192998/100‑1 (MQ=255)

                                                                  aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGa                                   >  2:2895125/1‑100 (MQ=255)

                                                                    ttttTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                 <  1:3018420/100‑1 (MQ=255)

                                                                    tttCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGTGTCGCGGCGGATc                                 <  1:2362941/100‑1 (MQ=255)

                                                                      tCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                               >  1:2744607/1‑100 (MQ=255)

                                                                       cTGTAATGCAGGAAAAACGATGCCATCGACACTAGTGACGATGACATGTGATGACAATGTTGATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                               <  1:2536282/99‑1 (MQ=255)

                                                                         gTAATGCAGGCAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGttt                            <  1:2507608/100‑1 (MQ=255)

                                                                             tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCg                        >  2:3805134/1‑100 (MQ=255)

                                                                             tGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCg                        >  2:919893/1‑100 (MQ=255)

                                                                              gCAGGAAAAACGATGTTATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGc                       <  1:2225658/100‑1 (MQ=255)

                                                                              gCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGc                       <  1:549221/100‑1 (MQ=255)

                                                                               cAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGCCAATGTTTATCGCGAAGGAGCAATGAGTGACCCGCGGCGGATCGGTTTTCCGc                       <  1:2664420/99‑1 (MQ=255)

                                                                                aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCACTGAGTGAGTCGCGGCGGATCAGTTTTCCGCTg                     <  2:3626393/100‑1 (MQ=255)

                                                                                 ggAAAAACGATGTCATCGACACTAGTGACGATGACATGAGATGACAATGTTTATAGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                    >  1:1876115/1‑100 (MQ=255)

                                                                                   aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                  >  2:1920504/1‑100 (MQ=255)

                                                                                   aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                  <  2:2494997/100‑1 (MQ=255)

                                                                                   aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                  >  2:1916768/1‑100 (MQ=255)

                                                                                   aaaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                  >  2:1772011/1‑100 (MQ=255)

                                                                                      aaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtt               >  1:1422553/1‑100 (MQ=255)

                                                                                         gATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc            <  2:1299421/100‑1 (MQ=255)

                                                                                         caTGTCATCGTCACTAGTGTCGATGACATGTGATGACACTGTTTATAGCGAAGAAGCAATGAGTGAGTCGAGGCGGATTAGTTTTCCGCTGAAGGTTTTc            >  2:2220360/2‑100 (MQ=37)

                                                                                         caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc            >  1:67884/2‑100 (MQ=255)

                                                                                         caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc            >  1:3239339/2‑100 (MQ=255)

                                                                                         caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc            >  2:2648341/2‑100 (MQ=255)

                                                                                         caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTACTGACCGATCAATTTTCAGCACAACGtttt             >  1:3072575/2‑99 (MQ=19)

                                                                                          aTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcg           >  1:499659/1‑100 (MQ=255)

                                                                                          aTGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcg           <  1:1612551/100‑1 (MQ=255)

                                                                                            gTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggc         >  1:716781/1‑100 (MQ=255)

                                                                                              cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg       >  2:1998743/1‑100 (MQ=255)

                                                                                              cATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg       >  1:3139116/1‑100 (MQ=255)

                                                                                               aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGt      <  1:135849/100‑1 (MQ=255)

                                                                                                gccacacTCGTCACAATGACATGTGATGACATTGTTAATCGCGAACGAGCAATAAGTGAGACGCGCCACACCACTCTTCCGCTGAATGTTTGcgga         <  1:3547896/93‑2 (MQ=255)

                                                                                                  gACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgag   <  2:351974/100‑1 (MQ=255)

                                                                                                   acacTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgaga  >  2:1034225/1‑100 (MQ=255)

ACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGA  >  NC_000913/2838123‑2838321

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