BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,567,248	G→C	A287G (GCG→GGG)	eutH ←	ethanolamine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,567,248	0	G	C	99.5%	632.4 / ‑4.6	183	A287G (GCG→GGG)	eutH	ethanolamine transporter
Reads supporting (aligned to +/- strand): ref base G (0/1); new base C (88/94); total (88/95)
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.47e-01

Only 100 of 186 total aligned reads displayed.

CATGCCGGCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACGCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAG  >  NC_000913/2567151‑2567342

cctGCTGGCTCCCGCGAGGTTGTGCTTGGTGAGTTCTGTACCGACGCTCATCAGCGGTTTTTCAAAGCAGCGAGTCAGCAGCAGCACCATCGGATAcccc                                                                                              <  2:4129482/98‑1 (MQ=255)

caTGCCGGCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATAcccc                                                                                              <  1:505749/100‑1 (MQ=255)

   gCCGGCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCt                                                                                           <  2:3509407/100‑1 (MQ=255)

      ggCTGCCGCGATGTTGTTCATATTCAGTCCTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAAc                                                                                        <  2:1047040/100‑1 (MQ=255)

       gCTGCCTCGATGTTGTCCATATTCAGTCCTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACa                                                                                       <  2:4993067/100‑1 (MQ=255)

       gCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACa                                                                                       <  2:917752/100‑1 (MQ=255)

       gCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACa                                                                                       <  2:1559346/100‑1 (MQ=255)

        cTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAg                                                                                      <  2:2855074/100‑1 (MQ=255)

          gCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGaa                                                                                    >  1:453231/1‑100 (MQ=255)

           ccGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAAc                                                                                   >  1:1512106/1‑100 (MQ=255)

            cgcgATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACg                                                                                  <  1:1703701/100‑1 (MQ=255)

            cgcgATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACg                                                                                  <  1:3636000/100‑1 (MQ=255)

            cgcgATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACg                                                                                  >  1:3104799/1‑100 (MQ=255)

             gcgATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGc                                                                                 <  2:3905292/100‑1 (MQ=255)

               gatcttgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAg                                                                               <  1:784734/96‑1 (MQ=255)

               gATGTTGATCATATTCAGTACTTTACCGACGCTCATCAGCGGTGTTTCAAACCAGCGAGTCAGCAGAAGCACCATGGGATAGCCCCGTAACAGAACGcgg                                                                               >  2:146398/1‑98 (MQ=255)

                 tgttgtTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGga                                                                             >  1:2937928/1‑100 (MQ=255)

                  gttgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTTTTTTTTTTTATTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGgag                                                                            >  2:1414430/1‑100 (MQ=255)

                  gttgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGgag                                                                            <  2:4414378/100‑1 (MQ=255)

                    tgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGAt                                                                          >  2:2417068/1‑100 (MQ=255)

                    tgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGAt                                                                          >  2:4719124/1‑100 (MQ=255)

                    tgttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGAt                                                                          <  1:926285/100‑1 (MQ=255)

                     gttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATa                                                                         >  1:3282478/1‑100 (MQ=255)

                      ttCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAg                                                                        <  2:2875971/100‑1 (MQ=255)

                        cATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGaa                                                                      >  1:4377967/1‑100 (MQ=255)

                        cATATTCAGCACTTTACCGCCGCTCATCGGCGGTTTTTCCACCCAGCGAGTCAGCCTCAGCACCCTCCGATACCCCCCTAAAAGAACGCAGGAGATAGaa                                                                      >  2:2112790/1‑100 (MQ=255)

                         atatTCAGTACTTTACCTACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAAc                                                                     <  1:3022491/100‑1 (MQ=255)

                         atatTCAGTACTTTACCAACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGaga                                                                     >  2:674037/1‑98 (MQ=255)

                          tatTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAAcc                                                                    >  1:79696/1‑100 (MQ=255)

                            ttCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGa                                                                  <  1:2262437/100‑1 (MQ=255)

                             tCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGAt                                                                 >  1:905819/1‑100 (MQ=255)

                             tCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGAt                                                                 >  1:2613570/1‑100 (MQ=255)

                              cAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATa                                                                >  1:564472/1‑100 (MQ=255)

                              cAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATa                                                                <  1:2700080/100‑1 (MQ=255)

                               aGTACTTTACCGACGCTCATCAGCGGTTTTTTTTTTTTTTTTTTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATaa                                                               >  2:2234716/1‑100 (MQ=255)

                                ggaCTTTACCGACGCTCATCAGCGGTTTTTCATGCGCGCGAGGCAGCAGCAGCACCATCGGACACCCGCCCAACAGATCGCAGGAGATAGGACCGATAAc                                                              >  2:4814861/3‑100 (MQ=255)

                                 tACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACGCCCCTAACAGAACGCAGGAGATAGAACCTATAACt                                                             <  2:3618551/100‑1 (MQ=255)

                                 tACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACt                                                             <  2:749154/100‑1 (MQ=255)

                                  aCTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACtt                                                            >  1:2530283/1‑100 (MQ=255)

                                  aCTTTACCGACGCTCAGCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACtt                                                            >  2:4910144/1‑100 (MQ=255)

                                    tttACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCa                                                          <  1:3674780/100‑1 (MQ=255)

                                     ttACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCaa                                                         >  1:3369971/1‑100 (MQ=255)

                                     ttACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCaa                                                         >  1:115779/1‑100 (MQ=255)

                                      tACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAAt                                                        <  1:1836987/100‑1 (MQ=255)

                                      tACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAAt                                                        <  2:837182/100‑1 (MQ=255)

                                       aCCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATg                                                       >  1:395184/1‑100 (MQ=255)

                                        ccGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATgg                                                      <  1:1031143/100‑1 (MQ=255)

                                          gACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGgcg                                                    <  1:1476795/100‑1 (MQ=255)

                                           aCGCTCATCAGCGGTTTTTCAAACAAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGgcgc                                                   <  2:4470594/100‑1 (MQ=255)

                                            cGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGgcgcg                                                  >  1:2134734/1‑100 (MQ=255)

                                              cTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCa                                                <  1:932832/100‑1 (MQ=255)

                                              cTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCa                                                >  1:4338086/1‑100 (MQ=255)

                                              cTCATCAGCGGTTTTTCAAACCAGCGAGTCACCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCa                                                >  1:4134808/1‑100 (MQ=255)

                                                catcaGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATATAACCTATAACTTCAATGGCGCGCAtc                                              >  2:3683873/1‑100 (MQ=255)

                                                catcaGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtc                                              >  1:1339438/1‑100 (MQ=255)

                                                  tcaGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcac                                            <  2:4954455/100‑1 (MQ=255)

                                                  tcaGCGGTTGTTCACACCAGCGAGTCAGCAGCCGCACCATCGGCTACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcac                                            <  1:1521168/100‑1 (MQ=255)

                                                   caGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacc                                           <  1:1250868/100‑1 (MQ=255)

                                                    aGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacct                                          <  1:4011212/100‑1 (MQ=255)

                                                     gCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCTCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcaccta                                         >  2:2618920/1‑99 (MQ=255)

                                                     gCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacctc                                         <  1:1489689/100‑1 (MQ=255)

                                                     gCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacctc                                         <  2:3801348/100‑1 (MQ=255)

                                                     gCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacctc                                         >  2:1422520/1‑100 (MQ=255)

                                                     gCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCAtcacctc                                         >  1:3602181/1‑100 (MQ=255)

                                                        gTTTCTCAAACCAGCGAGTCAGCATCAGCACCATCGGAGACCCCCCTAACAGAACGCAGTACATAGGACCGATAACGTTGACGTCGCGCAtcacctgacc                                      >  2:4308984/1‑100 (MQ=255)

                                                          ttttCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCgg                                    <  2:3070036/100‑1 (MQ=255)

                                                           tttCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCggg                                   >  1:1965721/1‑100 (MQ=255)

                                                            ttCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGt                                  >  1:4563354/1‑100 (MQ=255)

                                                             tCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGtt                                 <  2:2971119/100‑1 (MQ=255)

                                                             tCAAACCAGCGAGTCAGCAGCAGCACCAGCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGtt                                 <  1:2938746/100‑1 (MQ=255)

                                                              cAAACCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGttt                                >  2:814106/1‑100 (MQ=255)

                                                               aaaCCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTa                               >  1:1753561/1‑100 (MQ=255)

                                                               aaaCCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTa                               >  1:843117/1‑100 (MQ=255)

                                                                 aCCAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATc                             <  2:2392486/100‑1 (MQ=255)

                                                                  ccAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCg                            >  2:4365464/1‑100 (MQ=255)

                                                                   cAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGc                           >  1:4813831/1‑100 (MQ=255)

                                                                   cAGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGc                           <  2:3691900/100‑1 (MQ=255)

                                                                    aGCGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGcc                          <  1:460301/100‑1 (MQ=255)

                                                                      cGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAg                        >  2:255860/1‑100 (MQ=255)

                                                                      cGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAg                        <  2:4560492/100‑1 (MQ=255)

                                                                      cGAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGAGTTTATCGccgg                        <  1:4437286/100‑3 (MQ=255)

                                                                       gAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGACCGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTGTATCGCCAgg                       <  2:781042/100‑1 (MQ=255)

                                                                       gAGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAgg                       <  2:1490334/100‑1 (MQ=255)

                                                                        aGTCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAggg                      <  2:1730949/100‑1 (MQ=255)

                                                                          tCAGCAGCAGCACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGc                    >  2:116941/1‑100 (MQ=255)

                                                                            agcagcagcACCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCa                  <  2:1482514/100‑1 (MQ=255)

                                                                                gcagcaCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATaaa              <  2:3340126/100‑1 (MQ=255)

                                                                                 cagcaCCATCGGATGCCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAg             <  1:3391279/100‑1 (MQ=255)

                                                                                 cagcaCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCGCCTCACCGGGTTTATCGCCAGGGGCCATAAAg             <  2:2189336/100‑1 (MQ=255)

                                                                                 cagcaCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAg             <  1:3382453/100‑1 (MQ=255)

                                                                                 cagcaCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAg             >  1:2107940/1‑100 (MQ=255)

                                                                                  agcaCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGa            <  1:3393892/100‑1 (MQ=255)

                                                                                    caCCATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATa          <  1:1021569/100‑1 (MQ=255)

                                                                                      ccATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAgg        >  2:3931673/1‑100 (MQ=255)

                                                                                      ccATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAgg        <  1:3569666/100‑1 (MQ=255)

                                                                                      ccATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAgg        >  1:2141743/1‑100 (MQ=255)

                                                                                       cATCGGATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGa       >  2:2863066/1‑100 (MQ=255)

                                                                                          gagATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATcc    <  1:2871093/98‑1 (MQ=255)

                                                                                            gATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAg  <  2:3001877/100‑1 (MQ=255)

                                                                                            gATACCCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAg  <  2:2905148/100‑1 (MQ=255)

CATGCCGGCTGCCGCGATGTTGTTCATATTCAGTACTTTACCGACGCTCATCAGCGGTTTTTCAAACCAGCGAGTCAGCAGCAGCACCATCGGATACGCCCCTAACAGAACGCAGGAGATAGAACCGATAACTTCAATGGCGCGCATCACCTCACCGGGTTTATCGCCAGGGGCCATAAAGATAGGATCCAG  >  NC_000913/2567151‑2567342

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