BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,750,491	T→C	Y455C (TAC→TGC)	b1673 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,750,491	0	T	C	98.6%	433.6 / ‑3.7	139	Y455C (TAC→TGC)	b1673	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (69/67); total (71/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 149 total aligned reads displayed.

ACAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGGTAG  >  NC_000913/1750401‑1750580

 cAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCAct                                                                                        <  2:1087091/93‑1 (MQ=255)

  aaTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCActc                                                                                       >  1:3899941/1‑93 (MQ=255)

   aTGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCActct                                                                                      >  1:1037326/1‑93 (MQ=255)

    tGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCtt                                                                                     >  2:1720097/1‑93 (MQ=255)

    tGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCtt                                                                                     >  2:104923/1‑93 (MQ=255)

       aaCCGAACTGAGCTGCCGCTTCATTGGCATGGTGAACAGGATAGCGAAACGGCGACCAGAGTTTATTTTTCGAGTAGCGCCAGCACTCTTCAc                                                                                  >  2:2832858/1‑93 (MQ=255)

         ccGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCAccc                                                                                >  1:4389269/1‑93 (MQ=255)

          cGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCa                                                                               <  2:450161/93‑1 (MQ=255)

           gACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaa                                                                              >  1:2756600/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             >  1:3616686/1‑93 (MQ=255)

            aCCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCaaa                                                                             >  1:2272493/1‑93 (MQ=255)

              cTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAAtt                                                                           >  1:3380178/1‑93 (MQ=255)

               tGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTc                                                                          >  1:234848/1‑93 (MQ=255)

               tGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTc                                                                          >  1:555742/1‑93 (MQ=255)

               tGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTc                                                                          <  1:999298/93‑1 (MQ=255)

                 cgcTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCa                                                                        <  2:779908/93‑1 (MQ=255)

                  gcTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAg                                                                       <  1:4350794/93‑1 (MQ=255)

                   cTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGc                                                                      >  2:533197/1‑93 (MQ=255)

                     gACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGt                                                                    <  1:4301662/93‑1 (MQ=255)

                     gACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGt                                                                    >  1:1364102/1‑93 (MQ=255)

                     gACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGt                                                                    <  1:2812232/93‑1 (MQ=255)

                       cGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTc                                                                  <  1:1093668/93‑1 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCt                                                                 <  2:2087632/93‑1 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 >  1:3305433/1‑93 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 >  2:1454458/1‑93 (MQ=255)

                        gCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCt                                                                 <  1:202976/93‑1 (MQ=255)

                         cTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTg                                                                >  2:905290/1‑93 (MQ=255)

                          ttCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGc                                                               <  2:4627687/93‑1 (MQ=255)

                          ttCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGc                                                               <  1:1614113/93‑1 (MQ=255)

                              ttGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATg                                                           <  2:1497277/93‑1 (MQ=255)

                                ggCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGc                                                         >  1:1622914/1‑93 (MQ=255)

                                  cATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCat                                                       >  2:3062701/1‑93 (MQ=255)

                                   aTGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCata                                                      >  2:3923574/1‑93 (MQ=255)

                                    tGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCatat                                                     <  1:918755/93‑1 (MQ=255)

                                         aaCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAAcc                                                >  1:4538459/1‑93 (MQ=255)

                                         aaCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAAcc                                                <  2:3511843/93‑1 (MQ=255)

                                          aCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCa                                               >  1:3090696/1‑93 (MQ=255)

                                          aCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCa                                               <  1:3672271/93‑1 (MQ=255)

                                          aCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCa                                               <  1:3672394/93‑1 (MQ=255)

                                          aCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCa                                               <  2:3380343/93‑1 (MQ=255)

                                          aCCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGAACTCTTCACCCAAATTCTAGCGTTCTGCGATGGCATATGAAcaa                                               >  2:2490946/1‑91 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCACATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             >  1:4577663/1‑93 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             <  2:4372090/93‑1 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             <  2:1454011/93‑1 (MQ=255)

                                            cGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATc                                             <  1:2273813/93‑1 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            >  2:1580015/1‑93 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            >  2:1343628/1‑93 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            <  1:3769037/93‑1 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            <  2:3581391/93‑1 (MQ=255)

                                             ggATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCa                                            <  1:845098/93‑1 (MQ=255)

                                              gATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAg                                           <  2:2675326/93‑1 (MQ=255)

                                              gATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAg                                           <  2:1373947/93‑1 (MQ=255)

                                                tAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         <  2:2797186/93‑1 (MQ=255)

                                                tAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGcc                                         >  2:2025525/1‑93 (MQ=255)

                                                  gCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGGACTCTTCACTCAAATTCCGGCGTTCGGCGCTGGCATATGAACCATCAGCGAg                                       >  1:2711725/1‑93 (MQ=255)

                                                   ccAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAgg                                      >  1:917898/1‑93 (MQ=255)

                                                    cAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGt                                     <  1:360128/93‑1 (MQ=255)

                                                     aaaCGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTg                                    >  1:996912/1‑93 (MQ=255)

                                                     aaaCGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTg                                    >  1:4427253/1‑93 (MQ=255)

                                                      accGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGa                                   <  1:2572264/91‑1 (MQ=255)

                                                      aaCGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGa                                   >  2:1966937/1‑93 (MQ=255)

                                                      aaCGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGa                                   <  1:2572263/93‑1 (MQ=255)

                                                        cGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACt                                 >  2:2909897/1‑93 (MQ=255)

                                                         ggCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTc                                <  2:2213240/93‑1 (MQ=255)

                                                          gCGACCAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCa                               <  1:2632741/93‑1 (MQ=255)

                                                              ccAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTc                           >  1:481027/1‑93 (MQ=255)

                                                               cAGAGTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCa                          >  2:1148298/1‑93 (MQ=255)

                                                                agagTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAc                         >  2:557447/1‑93 (MQ=255)

                                                                agagTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAc                         >  2:937736/1‑93 (MQ=255)

                                                                 gagTTTATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCAcc                        >  2:2473031/1‑93 (MQ=255)

                                                                    tttATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCcac                     >  2:2270299/1‑93 (MQ=255)

                                                                    tttATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCcac                     <  2:1244617/93‑1 (MQ=255)

                                                                      tATTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCacac                   <  2:3684436/93‑1 (MQ=255)

                                                                       aTTTTTCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACg                  >  2:3277551/1‑93 (MQ=255)

                                                                        tttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGa                 <  2:4583711/93‑1 (MQ=255)

                                                                         ttttCGCGTAGCCCCAGCCCTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGAt                <  2:1025270/93‑1 (MQ=255)

                                                                         ttttCGCGTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGAt                <  2:604358/93‑1 (MQ=255)

                                                                              gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           >  1:370457/1‑93 (MQ=255)

                                                                              gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           >  2:186314/1‑93 (MQ=255)

                                                                              gcgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGcc           >  1:2338599/1‑93 (MQ=255)

                                                                               cgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCa          <  2:2310796/93‑1 (MQ=255)

                                                                               cgTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCa          >  2:2662378/1‑93 (MQ=255)

                                                                                gTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAg         >  1:650668/1‑93 (MQ=255)

                                                                                gTAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAg         <  1:2005207/93‑1 (MQ=255)

                                                                                 tAGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGa        <  1:3066230/93‑1 (MQ=255)

                                                                                  aGCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGAc       >  2:2334643/1‑93 (MQ=255)

                                                                                   gCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACg      <  1:929412/93‑1 (MQ=255)

                                                                                   gCCCCAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACg      <  1:1129038/93‑1 (MQ=255)

                                                                                    ccccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATAAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACgg     <  1:1542645/93‑1 (MQ=255)

                                                                                     cgcAGCACTCTTCACCCAAATCCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACGCGATGCGCCAGACGgt    <  2:325042/91‑1 (MQ=255)

                                                                                     cccAGCACTCTTCCCCCAAATTCCAGCGTTCTGCGATGGCATATGAGCCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgt    <  2:383281/93‑1 (MQ=255)

                                                                                     cccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgt    <  2:1787630/93‑1 (MQ=255)

                                                                                     cccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgt    >  1:2592422/1‑93 (MQ=255)

                                                                                      ccAGCACTCTTCGCCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   <  2:1651558/93‑1 (MQ=255)

                                                                                      ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgtg   >  1:2870038/1‑92 (MQ=255)

                                                                                      ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   >  2:3583792/1‑93 (MQ=255)

                                                                                      ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   <  1:716469/93‑1 (MQ=255)

                                                                                      ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   >  1:3859374/1‑93 (MQ=255)

                                                                                      ccAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgta   >  2:1983597/1‑93 (MQ=255)

                                                                                       cAGCACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGgtag  >  2:340231/1‑93 (MQ=255)

ACAATGGAACCGACCTGCGCTGACGCTTCATTGGCATGGTGAACCGGATAGCCAAACGGCGACCAGAGTTTATTTTTCGCGTAGCCCCAGTACTCTTCACCCAAATTCCAGCGTTCTGCGATGGCATATGAACCATCAGCCAGGTGACTCAGCTCACCCACACGATGCGCCAGACGGTAG  >  NC_000913/1750401‑1750580

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