BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,239,132	T→C	H186R (CAC→CGC)	galS ←	mgl repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,239,132	0	T	C	99.6%	909.2 / ‑4.9	260	H186R (CAC→CGC)	galS	mgl repressor
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (151/108); minor base A (1/0); total (152/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

Only 100 of 269 total aligned reads displayed.

GCGCCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGTGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAGGCA  >  NC_000913/2239042‑2239220

 cgcCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGc                                                                                       <  1:1046444/93‑1 (MQ=255)

   ccAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTg                                                                                     >  1:200056/1‑93 (MQ=255)

    cAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCACGGCTgg                                                                                    >  1:3377051/1‑93 (MQ=255)

     aaTCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGga                                                                                   <  2:2722909/93‑1 (MQ=255)

      aTCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaa                                                                                  >  1:3055668/1‑93 (MQ=255)

       tCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGTTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaag                                                                                 <  1:3554446/93‑1 (MQ=255)

       tCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaag                                                                                 <  2:3253431/93‑1 (MQ=255)

        ccAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaaga                                                                                <  2:3615627/93‑1 (MQ=255)

         cAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaagaa                                                                               <  1:978196/93‑1 (MQ=255)

         cAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGgaagaa                                                                               >  2:6068687/1‑93 (MQ=255)

          aGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaa                                                                              >  1:4083337/1‑93 (MQ=255)

          aGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaa                                                                              >  2:3909886/1‑93 (MQ=255)

            cTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGaagaaaga                                                                            <  1:3233524/93‑1 (MQ=255)

             tttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGAt                                                                           >  1:1491042/1‑93 (MQ=255)

              ttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATa                                                                          <  2:5452146/93‑1 (MQ=255)

              ttCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATa                                                                          <  1:5755001/93‑1 (MQ=255)

               tCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATaa                                                                         >  2:2489473/1‑93 (MQ=255)

                ccGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAAc                                                                        <  1:584078/93‑1 (MQ=255)

                  ggcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCa                                                                      <  1:2958815/93‑1 (MQ=255)

                  ggcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAAACa                                                                      >  2:3436585/1‑93 (MQ=255)

                   gcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCaa                                                                     >  2:556231/1‑93 (MQ=255)

                   gcggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCaa                                                                     <  2:1328107/93‑1 (MQ=255)

                    cggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAAt                                                                    <  1:3594548/93‑1 (MQ=255)

                    cggAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAAt                                                                    <  2:4071317/93‑1 (MQ=255)

                      gAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATAc                                                                  <  2:6265613/93‑1 (MQ=255)

                      gAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATAc                                                                  <  1:2332432/93‑1 (MQ=255)

                       aataatATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACg                                                                 >  2:2412971/1‑93 (MQ=255)

                        ataataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGt                                                                <  1:127156/93‑1 (MQ=255)

                        ataataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGt                                                                <  2:4406724/93‑1 (MQ=255)

                        ataataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGt                                                                <  2:4770862/93‑1 (MQ=255)

                         taataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGtt                                                               >  1:2615028/1‑93 (MQ=255)

                         taataTCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGtt                                                               <  1:4574843/93‑1 (MQ=255)

                           atatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtga                                                             >  2:5260392/1‑93 (MQ=255)

                            tatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATTGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgat                                                            >  1:6145432/1‑93 (MQ=255)

                            tatCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgat                                                            <  1:6170384/93‑1 (MQ=255)

                             atCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTtgatg                                                           <  2:90231/93‑1 (MQ=255)

                               ccTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGAc                                                         <  2:5340780/93‑1 (MQ=255)

                               ccTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGAc                                                         >  2:5784486/1‑93 (MQ=255)

                                cTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGAcc                                                        >  2:5100336/1‑93 (MQ=255)

                                 tGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCg                                                       <  2:4971018/93‑1 (MQ=255)

                                 tGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCg                                                       >  1:2286110/1‑93 (MQ=255)

                                 tGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCg                                                       <  1:4809300/93‑1 (MQ=255)

                                  gCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGt                                                      >  2:796106/1‑93 (MQ=255)

                                    tctTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTa                                                    >  2:6429319/1‑93 (MQ=255)

                                     ctTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTAt                                                   >  2:1866578/1‑93 (MQ=255)

                                     ctTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTAt                                                   <  1:4955152/93‑1 (MQ=255)

                                       ttCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTc                                                 <  1:4381997/93‑1 (MQ=255)

                                        tCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCa                                                >  2:758502/1‑93 (MQ=255)

                                         cAACTCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAg                                               <  2:2940673/93‑1 (MQ=255)

                                           aCGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCa                                             >  2:1603666/1‑93 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCCGCAACAt                                         >  1:2141336/1‑93 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACAt                                         <  2:3802389/93‑1 (MQ=255)

                                               aCTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACAt                                         <  2:167278/93‑1 (MQ=255)

                                                cTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATg                                        <  1:904113/93‑1 (MQ=255)

                                                cTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGAACGTTATTCAGCAACATg                                        >  1:1354579/1‑93 (MQ=255)

                                                 tCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATgc                                       <  2:5918491/93‑1 (MQ=255)

                                                    tCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAAGATGcgcg                                    >  1:2452828/1‑93 (MQ=255)

                                                    tCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGcgcg                                    >  1:4871712/1‑93 (MQ=255)

                                                     ccAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGt                                   >  1:326581/1‑93 (MQ=255)

                                                     ccAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGt                                   >  2:1095212/1‑93 (MQ=255)

                                                     ccAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGt                                   <  2:3686655/93‑1 (MQ=255)

                                                      cAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTc                                  <  2:3178629/93‑1 (MQ=255)

                                                      cAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTc                                  >  2:1462588/1‑93 (MQ=255)

                                                      cAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTc                                  >  2:2801239/1‑93 (MQ=255)

                                                        gCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGc                                <  2:2909407/93‑1 (MQ=255)

                                                        gCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGc                                >  2:4543883/1‑93 (MQ=255)

                                                        gCCTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGc                                <  2:1796699/93‑1 (MQ=255)

                                                         ccTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGcc                               <  1:2872664/93‑1 (MQ=255)

                                                         ccTGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGcc                               >  1:3825663/1‑93 (MQ=255)

                                                           tGCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAt                             <  1:5219578/93‑1 (MQ=255)

                                                            gCTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCAtt                            >  1:4668240/1‑93 (MQ=255)

                                                             cTTTACGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTc                           <  2:1074745/93‑1 (MQ=255)

                                                                  cGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCa                      >  2:2646556/1‑93 (MQ=255)

                                                                  cGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCa                      >  2:5225192/1‑93 (MQ=255)

                                                                  cGCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCa                      <  2:2777425/93‑1 (MQ=255)

                                                                   gCATGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCAc                     >  2:4467933/1‑93 (MQ=255)

                                                                      tGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGc                  <  1:1738071/93‑1 (MQ=255)

                                                                      tGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGc                  <  2:5518327/93‑1 (MQ=255)

                                                                      tGGCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGc                  >  2:4039185/1‑93 (MQ=255)

                                                                      tGGCGTCATCTTCAATGCCGCGGCTGGAACAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGc                  <  1:6354713/93‑1 (MQ=255)

                                                                       ggCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCt                 >  2:732291/1‑93 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                <  1:240963/93‑1 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                >  2:2686769/1‑93 (MQ=255)

                                                                        gCGTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTg                <  1:6419405/93‑1 (MQ=255)

                                                                          gTCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTgag              >  1:4988184/1‑93 (MQ=255)

                                                                           tCATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTgaga             >  2:1203433/1‑93 (MQ=255)

                                                                            cATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGAt            >  2:6245707/1‑93 (MQ=255)

                                                                            cATCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGAt            >  2:1984218/1‑93 (MQ=255)

                                                                             aTCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGAtt           >  1:1097010/1‑93 (MQ=255)

                                                                             aTCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATNCAGCAACATGCGCGTCNCCATTCGGGCACCNCNGAGNtt           >  2:2086312/1‑93 (MQ=255)

                                                                              tCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTa          <  1:1905207/93‑1 (MQ=255)

                                                                              tCTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTa          >  2:4278683/1‑93 (MQ=255)

                                                                               cTTCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTAt         >  2:6018030/1‑93 (MQ=255)

                                                                                ttCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATc        >  1:41896/1‑93 (MQ=255)

                                                                                ttCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATc        >  1:5644419/1‑93 (MQ=255)

                                                                                ttCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATc        >  2:2203933/1‑93 (MQ=255)

                                                                                 tCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATcc       <  1:1210873/93‑1 (MQ=255)

                                                                                 tCAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGAAACGCTGAGATTATcc       >  1:2738372/1‑93 (MQ=255)

                                                                                  cAATGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCa      >  1:5099023/1‑93 (MQ=255)

                                                                                     tGCCGCGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAgga   >  2:1249850/1‑92 (MQ=255)

GCGCCAATCCAGCTTTCCGGCGGAATAATATCCTGCTCTTTCAACGCACTCATCCAGCCTGCTTTACGCATGGCGTCATCTTCAATGCCGTGGCTGGAAGAAAGATAACCAATACGTTGATGACCGTTATTCAGCAACATGCGCGTCGCCATTCGGGCACCGCTGAGATTATCCAGGCA  >  NC_000913/2239042‑2239220

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