BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,298,707	T→C	T71A (ACC→GCC)	yhbT ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,298,707	0	T	C	99.5%	676.9 / NA	216	T71A (ACC→GCC)	yhbT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (103/112); minor base G (1/0); total (104/112)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

CAATCATCAGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGTAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAGCGCCTGGCGGA  >  NC_000913/3298617‑3298796

cAATCATCAGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCaa                                                                                         >  1:2236745/1‑93 (MQ=255)

 aaTCATCAGCAGATCGCTGGCCTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCaaa                                                                                        >  2:2983794/1‑93 (MQ=255)

  atcatcAGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAAc                                                                                       >  2:2914148/1‑93 (MQ=255)

  atcatcAGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAAc                                                                                       <  1:3907143/93‑1 (MQ=255)

    catcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCa                                                                                     <  1:4362589/93‑1 (MQ=255)

    catcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCa                                                                                     <  1:2536360/93‑1 (MQ=255)

     atcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCAt                                                                                    >  1:1294693/1‑93 (MQ=255)

     atcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCAt                                                                                    <  1:1533734/93‑1 (MQ=255)

      tcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCAtt                                                                                   <  1:1368773/93‑1 (MQ=255)

      tcaGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCGTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCAtt                                                                                   <  1:275256/93‑1 (MQ=255)

       cagcagATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTg                                                                                  <  2:684008/93‑1 (MQ=255)

         gcagATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCa                                                                                <  1:750572/93‑1 (MQ=255)

         gcagATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCa                                                                                <  2:3151702/93‑1 (MQ=255)

         gcagATCGCTGGCGTCGGCACGAAAACTCACAGCAGCGTGCGCTTTCTGGCTAACGGCCAGTTTTCCATTCACCACCGAGGCAAACCATTGCa                                                                                <  2:3135200/93‑1 (MQ=255)

          cagATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAg                                                                               <  1:2775821/93‑1 (MQ=255)

           agATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAgg                                                                              <  1:4309721/93‑1 (MQ=255)

           agATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAgg                                                                              <  2:412233/93‑1 (MQ=255)

             aTCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTc                                                                            <  2:3810412/93‑1 (MQ=255)

             aTCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTc                                                                            >  1:1140579/1‑93 (MQ=255)

                gCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAat                                                                         >  2:4107974/1‑93 (MQ=255)

                  tGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAatat                                                                       >  1:4490920/1‑93 (MQ=255)

                   ggCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATc                                                                      <  1:874520/93‑1 (MQ=255)

                   ggCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATc                                                                      >  1:3759099/1‑93 (MQ=255)

                    gcgcCGTCACTAAAACTCGCATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCa                                                                     <  1:2018672/93‑1 (MQ=255)

                    gCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCa                                                                     <  1:602695/93‑1 (MQ=255)

                    gCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCa                                                                     <  2:4526200/93‑1 (MQ=255)

                    gCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCa                                                                     <  1:2398159/93‑1 (MQ=255)

                     cGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCAc                                                                    >  1:2169923/1‑93 (MQ=255)

                       tCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGc                                                                  <  2:3720585/93‑1 (MQ=255)

                       gcGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGc                                                                  <  1:2865492/92‑1 (MQ=255)

                        cGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGca                                                                 <  2:3235385/93‑1 (MQ=255)

                         ggCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAACCCATTGCAGGTCAATATCACGcac                                                                <  2:4520164/93‑1 (MQ=255)

                         ggCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGcac                                                                >  2:4615443/1‑93 (MQ=255)

                          gCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGcaca                                                               <  2:1347281/93‑1 (MQ=255)

                          gCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGcaca                                                               <  1:442896/93‑1 (MQ=255)

                          gCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGcaca                                                               <  1:3218235/93‑1 (MQ=255)

                           cACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACAt                                                              >  2:732565/1‑93 (MQ=255)

                           cACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACAt                                                              >  1:445575/1‑93 (MQ=255)

                           cACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACAt                                                              >  1:2350350/1‑93 (MQ=255)

                           cACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTAACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACAt                                                              >  2:2462495/1‑93 (MQ=255)

                            aCTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATg                                                             <  1:3811707/93‑1 (MQ=255)

                             cTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGa                                                            >  1:3647880/1‑93 (MQ=255)

                              tAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGaa                                                           <  1:971384/93‑1 (MQ=255)

                                aaaCTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATa                                                         <  2:47793/93‑1 (MQ=255)

                                  aCTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACt                                                       <  1:785340/93‑1 (MQ=255)

                                  aCTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACt                                                       <  1:3165174/93‑1 (MQ=255)

                                  aCTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACt                                                       >  1:2898087/1‑93 (MQ=255)

                                    tCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTa                                                     >  1:4447721/1‑93 (MQ=255)

                                    tCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTa                                                     >  2:1184675/1‑93 (MQ=255)

                                     cacaTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTaa                                                    <  1:3552552/93‑1 (MQ=255)

                                     cacaTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTaa                                                    <  2:3531249/93‑1 (MQ=255)

                                     cacaTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTaa                                                    <  2:569307/93‑1 (MQ=255)

                                     cacaTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTaa                                                    <  2:1045721/93‑1 (MQ=255)

                                      acaTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGGCAATATCACGCACATGAATACTTGAc                                                   >  2:1832969/1‑93 (MQ=255)

                                       caTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAAcc                                                  <  2:911900/93‑1 (MQ=255)

                                       caTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAAcc                                                  <  2:2230580/93‑1 (MQ=255)

                                       caTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAAcc                                                  <  2:4082614/93‑1 (MQ=255)

                                       caTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAAcc                                                  <  1:3013960/93‑1 (MQ=255)

                                       caTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAAcc                                                  <  2:1554926/93‑1 (MQ=255)

                                        aTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCa                                                 <  1:1039226/93‑1 (MQ=255)

                                        aTCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCa                                                 <  2:3696651/93‑1 (MQ=255)

                                          cAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGc                                               <  2:4445637/93‑1 (MQ=255)

                                          cAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGc                                               >  1:1916658/1‑93 (MQ=255)

                                          cAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGc                                               <  1:1844978/93‑1 (MQ=255)

                                           aGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCg                                              >  1:3070740/1‑93 (MQ=255)

                                           aGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCg                                              >  2:282275/1‑93 (MQ=255)

                                             cTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGc                                            <  2:778688/93‑1 (MQ=255)

                                               tGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCt                                          >  2:656877/1‑93 (MQ=255)

                                               tGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCt                                          >  1:3431565/1‑93 (MQ=255)

                                                gcgcGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCtt                                         <  1:4056165/93‑1 (MQ=255)

                                                  gcgTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCa                                       <  2:2549158/93‑1 (MQ=255)

                                                    gTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAg                                     >  1:1204491/1‑93 (MQ=255)

                                                    gTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAg                                     <  1:906628/93‑1 (MQ=255)

                                                     ttCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGa                                    >  1:561980/1‑93 (MQ=255)

                                                     ttCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGa                                    >  1:516522/1‑93 (MQ=255)

                                                     ttCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGa                                    >  1:4164600/1‑93 (MQ=255)

                                                     ttCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGa                                    >  2:4212974/1‑93 (MQ=255)

                                                     ttCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGa                                    >  2:2716009/1‑93 (MQ=255)

                                                      tCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAAAATCACGCACATGAATACTTAACCAGCGGCCTTCAAGaa                                   >  1:1494531/1‑93 (MQ=255)

                                                          gCTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTc                               >  1:3539368/1‑93 (MQ=255)

                                                           cTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTcc                              <  1:3979934/93‑1 (MQ=255)

                                                           cTAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTcc                              <  2:1380345/93‑1 (MQ=255)

                                                            tAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCa                             >  1:1641519/1‑93 (MQ=255)

                                                            tAACGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCa                             >  2:2859221/1‑93 (MQ=255)

                                                              aCGACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGc                           <  2:3374668/93‑1 (MQ=255)

                                                                gACCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTc                         >  2:2931725/1‑93 (MQ=255)

                                                                 aCCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGGATACTTAACCAGCGGCCTTCAAGCAACTCCAGcgcg                        >  1:4094425/1‑93 (MQ=255)

                                                                 aCCAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCg                        >  2:1922181/1‑93 (MQ=255)

                                                                   cAGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGcc                      >  2:3548367/1‑93 (MQ=255)

                                                                    aGTTTGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCca                     <  2:221234/93‑1 (MQ=255)

                                                                       ttGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCcatca                  >  2:3989010/1‑93 (MQ=255)

                                                                        tGCCATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCcatcat                 <  2:2698345/93‑1 (MQ=255)

                                                                         tccATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCatcatc                <  1:813449/92‑1 (MQ=255)

                                                                          ccATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATcc               <  2:1532499/93‑1 (MQ=255)

                                                                           cATTCACCACCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCa              <  2:1718281/93‑1 (MQ=255)

                                                                               caccacCGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAgcgc          >  2:2752121/1‑93 (MQ=255)

                                                                                 ccaccGAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAGCGCCt        <  2:3319187/93‑1 (MQ=255)

                                                                                   accGAGGCGAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAGCGCCTgg      >  2:196948/1‑93 (MQ=255)

                                                                                      gAGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAGCGCCTGGCgg   >  1:658497/1‑93 (MQ=255)

                                                                                       aGGCAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCTAGCTCGCCATCATCCAGCGCCTGGCGGa  >  2:2117193/1‑93 (MQ=255)

CAATCATCAGCAGATCGCTGGCGTCGGCACTAAAACTCACATCAGCTTGCGCGTTCTGGCTAACGACCAGTTTGCCATTCACCACCGAGGTAAACCATTGCAGGTCAATATCACGCACATGAATACTTAACCAGCGGCCTTCAAGAAACTCCAGCTCGCCATCATCCAGCGCCTGGCGGA  >  NC_000913/3298617‑3298796

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