BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,838,231	A→T	17.0%	intergenic (‑126/+23)	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,231	0	A	T	17.0%	331.8 / 55.2	135	intergenic (‑126/+23)	hydN/ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand): ref base A (59/52); new base T (12/11); total (72/63)
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.47e-01

Only 100 of 137 total aligned reads displayed.

CGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGC  >  NC_000913/2838133‑2838330

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

   cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGAt                                                                                                 >  2:744941/1‑100 (MQ=255)

    gATACATTCACAAATCATGCTAGTTTGTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                <  1:865399/100‑1 (MQ=255)

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

    gATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                >  2:1866103/1‑100 (MQ=255)

    gATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                >  1:1383779/1‑100 (MQ=255)

    gATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                <  1:587534/100‑1 (MQ=255)

       aCATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACa                                                                                             <  2:1453510/100‑1 (MQ=255)

        cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGAcgt                                                                                            <  2:1970075/100‑3 (MQ=255)

        cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAt                                                                                            <  2:2401431/100‑1 (MQ=255)

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

          ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtcg                                                                                          <  2:2078159/100‑3 (MQ=255)

          ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtcg                                                                                          <  2:698767/100‑3 (MQ=255)

          ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                          >  1:2597384/1‑100 (MQ=255)

          ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGCAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                          <  2:1677470/100‑1 (MQ=255)

          ttCACAAATCATGCCAGTTTTTAATGTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtcg                                                                                          <  1:1392191/100‑3 (MQ=38)

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

            cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtcg                                                                                          <  2:1351341/98‑3 (MQ=255)

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

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

            cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGCGACGATGACATGtga                                                                                        >  1:2272126/1‑100 (MQ=255)

             acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgat                                                                                       <  2:2736529/100‑1 (MQ=255)

             acaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgat                                                                                       <  1:822335/100‑1 (MQ=255)

               aaaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                     <  1:2557978/100‑1 (MQ=255)

                  tCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACaa                                                                                  >  1:697238/1‑100 (MQ=255)

                   cATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAGTTTCTGTAATGCAGGAAAAACGATGTCATCGACAGTAGTGACGATGACATGTGATGACa                                                                                   >  1:1801483/1‑98 (MQ=255)

                   cATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAAt                                                                                 >  2:44670/1‑100 (MQ=255)

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

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

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

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

                        cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGTAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTa                                                                            <  1:2780416/100‑1 (MQ=255)

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

                            ttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgc                                                                        <  2:2500127/100‑1 (MQ=255)

                             tttAATTTACTGTTATTTAAGGAAATTATTTTCTGTAATAAAAAAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcgcg                                                                       <  1:631758/100‑1 (MQ=255)

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

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

                                aaTTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAg                                                                    <  1:1632602/100‑1 (MQ=255)

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

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

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

                                  tgtACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATTGAGACTAGTGACGATGACATGTGATGACAAGGTTTATCGCGAAGGa                                                                  <  1:2725079/98‑1 (MQ=255)

                                      cTGTTATTTAAGGCAATTAATATCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGCGCaa                                                              >  1:74683/1‑100 (MQ=255)

                                      cTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGCCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                              >  1:2279036/1‑100 (MQ=255)

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

                                       tGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAt                                                             <  2:974742/100‑1 (MQ=255)

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

                                          tATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                          >  2:353194/1‑100 (MQ=255)

                                           aTTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                         >  1:770840/1‑100 (MQ=255)

                                           aTTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagt                                                         <  2:1522820/100‑1 (MQ=255)

                                            tttAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtg                                                        >  1:75150/1‑100 (MQ=255)

                                               aaGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATgagtgagt                                                     >  2:1382351/1‑100 (MQ=255)

                                                   aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTcgcg                                                 <  2:1117258/100‑1 (MQ=255)

                                                        aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGa                                            <  2:2612109/100‑1 (MQ=255)

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

                                                         attTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                           >  2:503965/1‑100 (MQ=255)

                                                         attTCTGTAATGCAGGAAAAACGATGTCATCGACACGAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                           <  2:1838418/100‑1 (MQ=255)

                                                          tttCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                          >  2:2714529/1‑100 (MQ=255)

                                                           ttCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                         <  1:2444632/100‑1 (MQ=255)

                                                            tCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                        <  1:918321/100‑1 (MQ=255)

                                                            tCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAg                                        <  1:1397050/100‑1 (MQ=255)

                                                              tGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtt                                      >  2:1722839/1‑100 (MQ=255)

                                                               gTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGttt                                     <  2:1925540/100‑1 (MQ=255)

                                                                tAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtttt                                    >  1:917492/1‑100 (MQ=255)

                                                                  aTGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                  >  2:1508207/1‑100 (MQ=255)

                                                                  aTGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                  >  1:1938615/1‑100 (MQ=255)

                                                                  aTGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                  >  2:2295844/1‑100 (MQ=255)

                                                                   ttcAGGAAAAACTATGTCATCGACACTAGTGACGATGGCATGTGATAACGATGTTTATCGCGAAGGAACAATGAGTGAGTCGCGGCGGCTCAGTTTTcc                                  >  1:1615614/3‑99 (MQ=255)

                                                                     ctggAAAAACGATATCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCt                               >  2:1977247/3‑100 (MQ=255)

                                                                     cAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCt                               >  1:1012651/1‑100 (MQ=255)

                                                                      aGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTg                              <  1:1632249/100‑1 (MQ=255)

                                                                       ggAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGa                             >  2:393369/1‑100 (MQ=255)

                                                                         ttccgatctGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                           <  1:2818419/92‑1 (MQ=38)

                                                                         tcatcacaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                           >  1:2442855/8‑100 (MQ=38)

                                                                           atcacaTGTCATCGTCACTGGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGGGTGAGTCGCGGGGGATGGGTGTTCGGATGAAgg                          >  1:2521880/6‑99 (MQ=21)

                                                                              acaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAACGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGtttt                      >  2:1241292/3‑100 (MQ=38)

                                                                               caTGTCATCGTCACTAGTGTCGATGACATGTGATGACATTGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTc                     >  1:1075783/2‑100 (MQ=38)

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

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

                                                                               caTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCAGCTGAAGGTTTTc                     >  1:862894/2‑100 (MQ=255)

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

                                                                                aTGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTGTTCCGCTGAAGGTTTTcg                    <  1:945136/100‑1 (MQ=255)

                                                                                 tGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGttt                       >  1:1105312/1‑96 (MQ=255)

                                                                                 tGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGATGGTTTTcgg                   <  2:51613/100‑1 (MQ=255)

                                                                                  ntCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTGTTCCGCTGAAGGTTTTcggc                  >  1:721994/2‑100 (MQ=255)

                                                                                  gTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAgagatcggaa                  >  2:2818437/1‑91 (MQ=38)

                                                                                   tCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGATCAATGAGGGAGTCGCGGCGGATCAGTTTTCCGCTGAAGATTTTcggcg                 >  1:602600/1‑100 (MQ=255)

                                                                                    cATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAACGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg                <  1:1241287/100‑1 (MQ=255)

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

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

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

                                                                                     aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggcgg                >  1:353227/1‑99 (MQ=255)

                                                                                     aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGt               >  1:2069973/1‑100 (MQ=255)

                                                                                     aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGt               >  2:1317726/1‑100 (MQ=255)

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

                                                                                               gTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATccc     >  2:2370809/1‑100 (MQ=255)

                                                                                                tGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc    >  2:114279/1‑100 (MQ=255)

                                                                                                 gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg   >  1:2416305/1‑100 (MQ=255)

                                                                                                  aCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGc  <  1:2030795/100‑1 (MQ=255)

                                                                                                  aCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCACTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGc  <  1:1790005/100‑1 (MQ=255)

CGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGC  >  NC_000913/2838133‑2838330

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: `‑`