BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,205,438	T→C	intergenic (‑328/‑287)	purH ← / → rrsE	phosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,205,438	0	T	C	100.0%	534.3 / NA	170	intergenic (‑328/‑287)	purH/rrsE	phosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (97/73); total (97/73)

Only 100 of 179 total aligned reads displayed.

TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAG  >  NC_000913/4205348‑4205512

 aaaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCa                                                                         >  2:442313/1‑93 (MQ=255)

  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCAt                                                                        >  1:4852183/1‑93 (MQ=255)

  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCAt                                                                        <  1:5631518/93‑1 (MQ=255)

   aGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATc                                                                       >  2:5786467/1‑93 (MQ=255)

    gCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCg                                                                      >  1:1811781/1‑93 (MQ=255)

     cGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGa                                                                     >  2:2494025/1‑93 (MQ=255)

       ggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGaca                                                                   >  1:92206/1‑93 (MQ=255)

        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGCGAACTCCCTATAATGCGCCCCCATCGacac                                                                  >  1:3420215/1‑93 (MQ=255)

        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                  >  2:4848683/1‑93 (MQ=255)

        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                  <  1:1648495/93‑1 (MQ=255)

        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                  <  1:3562654/93‑1 (MQ=255)

        gCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGacac                                                                  <  2:1026808/93‑1 (MQ=255)

          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                >  2:1759986/1‑93 (MQ=255)

          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                >  2:1621134/1‑93 (MQ=255)

          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                >  2:2174708/1‑93 (MQ=255)

          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcgg                                                                >  2:730668/1‑93 (MQ=255)

           tttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGACCTCCCTATAATGCGCCCCCATCGACAcggc                                                               <  1:5549703/93‑1 (MQ=255)

           tttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggc                                                               >  1:5159658/1‑93 (MQ=255)

           tttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggc                                                               <  2:1406430/93‑1 (MQ=255)

             tGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcgg                                                             >  1:1518193/1‑93 (MQ=255)

              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcggt                                                            >  2:1239331/1‑92 (MQ=255)

                tGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtg                                                          >  2:4835552/1‑93 (MQ=255)

                tGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATAGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtg                                                          >  1:3384955/1‑93 (MQ=255)

                  gggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGAtgtg                                                        >  2:4998759/1‑93 (MQ=255)

                   ggAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGa                                                       >  2:3487637/1‑93 (MQ=255)

                    gAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGaa                                                      >  2:3298306/1‑93 (MQ=255)

                      aTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGGAGTGAATc                                                    >  2:4342512/1‑93 (MQ=255)

                      aTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATc                                                    >  2:3640970/1‑93 (MQ=255)

                       tGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCa                                                   <  1:2885511/93‑1 (MQ=255)

                        gcagcaGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCAc                                                  <  1:4341287/93‑1 (MQ=255)

                         cagcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACt                                                 <  1:3826655/93‑1 (MQ=255)

                          agcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACtt                                                >  2:4004059/1‑93 (MQ=255)

                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                               >  1:5030961/1‑93 (MQ=255)

                           gcagTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                               <  1:6264596/93‑1 (MQ=255)

                            cagtcagtCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTca                                              <  1:4083179/93‑1 (MQ=255)

                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGGACTCCCTATAATGCGACCCCATCGACACGGCCGATGTGAATCCCTTcac                                             >  2:5413718/1‑93 (MQ=25)

                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                             <  1:1787125/93‑1 (MQ=255)

                              gtcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcaca                                            <  1:4519718/93‑1 (MQ=255)

                               tcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacac                                           >  1:3188885/1‑93 (MQ=255)

                                cagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacaca                                          >  2:520901/1‑93 (MQ=255)

                                cagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacaca                                          >  1:3304507/1‑93 (MQ=255)

                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                         >  2:16595/1‑93 (MQ=255)

                                 agtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaa                                         <  2:6015582/93‑1 (MQ=255)

                                  gtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACaaa                                        >  1:4189848/1‑93 (MQ=255)

                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                       >  1:337315/1‑93 (MQ=255)

                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                       <  1:2505462/93‑1 (MQ=255)

                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAc                                       <  2:1249534/93‑1 (MQ=255)

                                    caTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACa                                      <  1:6232672/93‑1 (MQ=255)

                                    caTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACa                                      <  2:2947582/93‑1 (MQ=255)

                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                     >  2:2346036/1‑93 (MQ=255)

                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                     >  1:1016623/1‑93 (MQ=255)

                                     aTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAg                                     >  1:3925755/1‑93 (MQ=255)

                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    <  1:2419908/93‑1 (MQ=255)

                                      tttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGc                                    >  2:4910932/1‑93 (MQ=255)

                                       ttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGcc                                   >  1:3472870/1‑93 (MQ=255)

                                        tttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcg                                  >  1:265710/1‑93 (MQ=255)

                                          tCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggt                                >  2:5599369/1‑93 (MQ=255)

                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               <  2:403698/93‑1 (MQ=255)

                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               >  2:2606350/1‑93 (MQ=255)

                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               <  2:2893685/93‑1 (MQ=255)

                                           cTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggtt                               >  2:4304011/1‑93 (MQ=255)

                                            tgtaaTTTTTATTTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCCCTTCACACAAACAGCcggttc                              <  2:2018193/90‑1 (MQ=17)

                                            tGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttc                              <  2:2579569/93‑1 (MQ=255)

                                            tGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttc                              <  1:1853561/93‑1 (MQ=255)

                                             gCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTACAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcg                             >  1:760911/1‑93 (MQ=255)

                                              cAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcgg                            <  1:2327804/93‑1 (MQ=255)

                                               aaTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcgggtcggt                           >  2:5444530/1‑93 (MQ=255)

                                                 tttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTg                         >  2:6252539/1‑93 (MQ=255)

                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                       >  2:2071626/1‑93 (MQ=255)

                                                   tttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGaa                       >  1:5919811/1‑93 (MQ=255)

                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAAAAGCCTGTTCGGTTgtaa                      >  2:591361/1‑90 (MQ=18)

                                                     tCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAaga                     >  2:1571449/1‑93 (MQ=255)

                                                     tCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAaga                     >  2:3699287/1‑93 (MQ=255)

                                                       tATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAgaga                   >  1:1057831/1‑93 (MQ=255)

                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                  >  2:3736859/1‑93 (MQ=255)

                                                        aTTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                  >  2:3287594/1‑93 (MQ=255)

                                                         ttGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaaa                 >  2:4439761/1‑93 (MQ=255)

                                                         ttGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGGGAATCACTTCACACAAACAGCCGGCTCGGTTGAAGAGaaa                 >  1:1671177/1‑93 (MQ=255)

                                                          tGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaaaa                >  2:4741123/1‑93 (MQ=255)

                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATc             >  1:1407808/1‑93 (MQ=255)

                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATc             <  2:2792985/93‑1 (MQ=255)

                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATc             >  1:4118514/1‑93 (MQ=255)

                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc            >  1:4003682/1‑93 (MQ=255)

                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc            >  2:5908401/1‑93 (MQ=255)

                                                              gCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATcc            >  1:1878618/1‑93 (MQ=255)

                                                               ccTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCt           <  2:4526458/93‑1 (MQ=255)

                                                                cTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAACCAGCCGGTTCGGTTGAAGAGAAAAATCCTg          >  1:1546141/1‑93 (MQ=255)

                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa         <  2:3908697/93‑1 (MQ=255)

                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa         >  1:3708750/1‑93 (MQ=255)

                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa         <  2:1363559/93‑1 (MQ=255)

                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt      <  2:4581302/93‑1 (MQ=25)

                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt      >  2:4531476/1‑93 (MQ=25)

                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt      <  1:4720345/93‑1 (MQ=25)

                                                                    ggAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAt      >  2:5648446/1‑93 (MQ=25)

                                                                      agaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    >  1:911295/1‑93 (MQ=18)

                                                                      agaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    <  2:5432290/93‑1 (MQ=18)

                                                                      agaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    <  1:3556691/93‑1 (MQ=18)

                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa   <  2:3607078/93‑1 (MQ=17)

                                                                       gaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCa   >  2:1513531/1‑93 (MQ=17)

                                                                        aaCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  <  1:2062707/93‑1 (MQ=14)

TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAG  >  NC_000913/4205348‑4205512

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