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%	542.2 / NA	176	intergenic (‑328/‑287)	purH/rrsE	phosphoribosylaminoimidazolecarboxamideformyltra nsferase/16S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (104/72); total (104/72)

Only 100 of 192 total aligned reads displayed.

TAAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCTCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAG  >  NC_000913/4205348‑4205512

taAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCcccc                                                                          >  2:49345/1‑93 (MQ=255)

taAAGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCcccc                                                                          >  2:4374671/1‑93 (MQ=255)

 aaaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCccccc                                                                         >  2:2031736/1‑92 (MQ=255)

  aaGCGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGCGAACTCCCTATAATGCGCCCCCAt                                                                        >  1:936504/1‑93 (MQ=255)

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

     cGGGCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGa                                                                     <  2:4266595/93‑1 (MQ=255)

      gggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGac                                                                    >  2:2667448/1‑93 (MQ=255)

       ggCTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATTATGCGCCCCCATCGaca                                                                   >  1:3342481/1‑93 (MQ=255)

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

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

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

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

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

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

         cTTTTGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGTGGCCTGCGGAGTACTCCCTATAATGCGCCCCCATCGACAcg                                                                 >  1:275403/1‑93 (MQ=255)

          ttttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGCCAcgg                                                                >  1:1112997/1‑93 (MQ=255)

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

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

            ttGCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcg                                                              >  2:2930708/1‑93 (MQ=255)

              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACAcggcggg                                                            >  2:102101/1‑92 (MQ=255)

              gCTGGGGAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGa                                                            >  1:4176713/1‑93 (MQ=255)

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

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

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

                    gAATGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGaa                                                      <  1:359238/93‑1 (MQ=255)

                      aTGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATc                                                    <  1:12164/93‑1 (MQ=255)

                       tGCAGCAGTCAGTCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAAtta                                                   >  2:2481047/1‑91 (MQ=255)

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

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

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

                           acagtcagtCATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTc                                               <  1:3670819/92‑1 (MQ=255)

                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                             >  2:2527035/1‑93 (MQ=255)

                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                             >  1:2409789/1‑93 (MQ=255)

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

                             agtcagtcATTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcac                                             >  1:4163207/1‑93 (MQ=255)

                              gtcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGGACTCCCTATAATGCGCCCCCCTCCACACCGCGGGTGTGGATGAGTTcaca                                            >  2:4081363/1‑93 (MQ=25)

                              gtcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATTTGAATCACTTcaca                                            >  1:1048458/1‑93 (MQ=255)

                               tcagtcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTcacac                                           <  2:1478489/93‑1 (MQ=255)

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

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

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

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

                                   tcaTTTTTCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACCAAc                                       >  1:3960183/1‑93 (MQ=255)

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

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

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

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

                                       ttttttGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGcc                                   <  2:1639729/93‑1 (MQ=255)

                                       ttttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGcc                                   <  2:2628385/93‑1 (MQ=255)

                                        tttCTGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGGGAATCACTTCACACAAACAGCcg                                  >  1:3752775/1‑93 (MQ=255)

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

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

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

                                            tGCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttc                              >  1:3166883/1‑93 (MQ=255)

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

                                             gCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcg                             >  1:653737/1‑93 (MQ=255)

                                             gCAATTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcg                             <  1:2897669/93‑1 (MQ=255)

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

                                               aaTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggt                           >  2:1323746/1‑93 (MQ=255)

                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                          <  2:269028/93‑1 (MQ=255)

                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                          <  1:2098546/93‑1 (MQ=255)

                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCcggttcggtt                          <  2:119325/93‑1 (MQ=255)

                                                aTTTTTCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATAACTTCACACAAACAGCcggttcggtt                          >  1:3617770/1‑93 (MQ=255)

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

                                                  ttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGa                        >  1:1793449/1‑93 (MQ=255)

                                                  ttttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGa                        >  1:3095704/1‑93 (MQ=255)

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

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

                                                    ttCTATTGCGGCCTGCGGAGAGCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                      <  1:4231034/93‑1 (MQ=255)

                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                      >  2:1755229/1‑93 (MQ=255)

                                                    ttCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAag                      >  2:1679110/1‑93 (MQ=255)

                                                     tCTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAaga                     <  2:4180406/93‑1 (MQ=255)

                                                      cTATTGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAagag                    >  2:337356/1‑93 (MQ=255)

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

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

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

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

                                                        aTTGCAGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaa                  >  2:3055186/1‑93 (MQ=255)

                                                         ttGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaaa                 <  2:3115090/93‑1 (MQ=255)

                                                          tGCGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGaaaa                <  2:409067/93‑1 (MQ=255)

                                                            cGGCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAAt              <  1:528605/93‑1 (MQ=255)

                                                             ggCCTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACCAACAGCCGGTTCGGTTGCAGAGAAAACTc             >  2:94710/1‑93 (MQ=14)

                                                              gCCTGCGGAGATCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAATCAGCCGGTTCGGTTGAAGAGAAAAATcc            <  1:4607114/93‑1 (MQ=255)

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

                                                                cTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTg          >  2:2212611/1‑93 (MQ=255)

                                                                cTGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGATCGGTTGAAGAGAAAAATCCTg          >  2:2763040/1‑93 (MQ=255)

                                                                 tGCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATATGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGa         >  1:325250/1‑93 (MQ=255)

                                                                  gCGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaa        <  2:2684564/93‑1 (MQ=255)

                                                                   cGGAGAACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGaaa       <  1:2043574/93‑1 (MQ=37)

                                                                     gagaACTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAAtt     <  1:1169549/93‑1 (MQ=21)

                                                                      agaCCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACCCAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTc    <  2:1353389/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                        aaCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  >  1:2503266/1‑93 (MQ=14)

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

                                                                        aaCTCCCTATAATGCGCCCCCATCGACACGGCGGATGTGAATCACTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAg  >  1:767850/1‑93 (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` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`