BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,754,821	T→C	T212A (ACT→GCT)	yiaY ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,754,821	0	T	C	99.6%	869.3 / ‑4.8	268	T212A (ACT→GCT)	yiaY	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (129/138); total (129/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

Only 100 of 287 total aligned reads displayed.

CTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATAC  >  NC_000913/3754731‑3754910

cTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGa                                                                                         <  1:6376567/93‑1 (MQ=255)

 ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGAt                                                                                        >  1:1068890/1‑93 (MQ=255)

 ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGAt                                                                                        >  2:6620241/1‑93 (MQ=255)

 ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGAt                                                                                        <  2:5829568/93‑1 (MQ=255)

 ttCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGAt                                                                                        >  1:1887203/1‑93 (MQ=255)

  tCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATc                                                                                       <  1:5726492/93‑1 (MQ=255)

  tCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATc                                                                                       <  2:4479234/93‑1 (MQ=255)

   cgcgcgCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCg                                                                                      <  2:1972096/91‑1 (MQ=255)

   cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCg                                                                                      <  2:2780088/93‑1 (MQ=255)

   cACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCg                                                                                      <  2:5154188/93‑1 (MQ=255)

    aCGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCgg                                                                                     <  1:2500596/93‑1 (MQ=255)

    aCGCGCTTTCGCATGACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGGCAGCGATCgg                                                                                     <  1:3768579/93‑1 (MQ=255)

     cgcgCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCggg                                                                                    >  2:3747205/1‑92 (MQ=255)

      gcgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCg                                                                                   <  1:1657150/93‑1 (MQ=255)

       cgcTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGt                                                                                  <  1:1402725/93‑1 (MQ=255)

         cTTTCGCATTTCTGCCATGTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTgg                                                                                <  2:5463135/93‑1 (MQ=255)

         cTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTgg                                                                                >  2:6381157/1‑93 (MQ=255)

         cTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTgg                                                                                <  2:2208411/93‑1 (MQ=255)

          tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggc                                                                               <  2:1143722/93‑1 (MQ=255)

          tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggc                                                                               >  1:810737/1‑93 (MQ=255)

          tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggc                                                                               >  2:6406657/1‑93 (MQ=255)

          tttCGCATTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggg                                                                               >  1:3377330/1‑92 (MQ=255)

            tCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggcgg                                                                             <  2:4700346/93‑1 (MQ=255)

             cGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTggcggc                                                                            >  2:411164/1‑93 (MQ=255)

               cATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCaa                                                                          >  2:964692/1‑93 (MQ=255)

                aTTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAAt                                                                         <  2:6468399/93‑1 (MQ=255)

                aTTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAAt                                                                         >  1:4263594/1‑93 (MQ=255)

                aTTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAAt                                                                         <  1:5774653/93‑1 (MQ=255)

                aTTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAAt                                                                         >  1:5080979/1‑93 (MQ=255)

                aTTACTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGGAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCCAt                                                                         >  2:3551179/1‑93 (MQ=255)

                 ttACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATa                                                                        <  2:2718770/93‑1 (MQ=255)

                 ttACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATa                                                                        >  2:522860/1‑93 (MQ=255)

                 ttACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATa                                                                        <  2:6643958/93‑1 (MQ=255)

                 ttACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATa                                                                        <  1:886944/93‑1 (MQ=255)

                 ttACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCCGCGTGGCGGCAATa                                                                        >  1:23693/1‑93 (MQ=255)

                  tACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAg                                                                       >  2:6710291/1‑93 (MQ=255)

                   aCTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGa                                                                      <  1:4850278/93‑1 (MQ=255)

                   aCTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGa                                                                      <  1:1854571/93‑1 (MQ=255)

                   aCTGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGa                                                                      >  1:214839/1‑93 (MQ=255)

                    cTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGaa                                                                     <  2:5938762/93‑1 (MQ=255)

                     tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCCGCGATCGGCGTGGCGGCAATAGaaa                                                                    >  1:3899618/1‑93 (MQ=255)

                     tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGaaa                                                                    <  2:4187915/93‑1 (MQ=255)

                     tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGaaa                                                                    <  1:2253121/93‑1 (MQ=255)

                     tGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGaaa                                                                    >  2:2498788/1‑93 (MQ=255)

                     tGCCATCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGaaa                                                                    >  1:3465914/1‑93 (MQ=255)

                      gCCATCTTCAAGGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAAc                                                                   <  1:3861351/93‑1 (MQ=255)

                       ccATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACa                                                                  >  2:2644004/1‑93 (MQ=255)

                        cATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCCCAAGCGTCAGCGATCGGCGCGGCGGCAATAGaaaaca                                                                 >  2:2369227/1‑90 (MQ=255)

                        cATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACaa                                                                 >  2:5653392/1‑92 (MQ=255)

                         aTCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACata                                                                <  1:6274832/93‑1 (MQ=255)

                         aTCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACata                                                                <  1:6085178/93‑1 (MQ=255)

                          tCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACatat                                                               >  2:2933537/1‑93 (MQ=255)

                          tCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACatat                                                               >  2:2449692/1‑93 (MQ=255)

                          tCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACatat                                                               <  2:5345149/93‑1 (MQ=255)

                          tCTTCAACGGCTAACGGCAGGGTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACatat                                                               >  2:5280811/1‑93 (MQ=255)

                           cTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATg                                                              >  1:1594997/1‑93 (MQ=255)

                            ttCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGc                                                             >  1:3627443/1‑93 (MQ=255)

                             tCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCt                                                            <  2:1292394/93‑1 (MQ=255)

                              cAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCtt                                                           >  1:2049584/1‑93 (MQ=255)

                              cAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCtt                                                           <  1:5472326/93‑1 (MQ=255)

                               aaCGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCTGCAATAGAAACATATGCTTc                                                          <  1:6570271/93‑1 (MQ=255)

                               aaCGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTc                                                          >  1:5296551/1‑93 (MQ=255)

                                aCGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCg                                                         >  2:1735101/1‑93 (MQ=255)

                                    cTAACGGCAGGTTTTTGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAg                                                     >  1:6159653/1‑93 (MQ=255)

                                    cTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAAAATATGCTTCGATAg                                                     >  2:1512845/1‑93 (MQ=255)

                                     tAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAgc                                                    >  2:5547282/1‑93 (MQ=255)

                                     tAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAgc                                                    <  1:2675541/93‑1 (MQ=255)

                                            aGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGtt                                             <  1:1328593/93‑1 (MQ=255)

                                             ggTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTTGATAGCGTGCGTTa                                            >  1:4799437/1‑93 (MQ=255)

                                              gtgttCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGGCAGCAATCGGCGTGGCGGCAATAGAAACATATGCTCCGATAGCGTGCGTTaa                                           <  1:3144150/90‑1 (MQ=255)

                                              gTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTaa                                           >  1:985378/1‑93 (MQ=255)

                                               ttttCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAg                                          <  2:6343782/93‑1 (MQ=255)

                                               ttttCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAg                                          <  1:2929300/93‑1 (MQ=255)

                                                    ggCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATc                                     <  1:1484312/93‑1 (MQ=255)

                                                     gCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATcc                                    >  2:5380268/1‑93 (MQ=255)

                                                       aaTCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATCGAAACATATGCTTCGATAGCGTGCGTTAAGGCATccct                                  >  1:3109089/1‑91 (MQ=255)

                                                        aTCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATa                                 >  1:2684007/1‑93 (MQ=255)

                                                        aTCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATa                                 <  2:3480118/93‑1 (MQ=255)

                                                        aTCATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATa                                 >  1:5133736/1‑93 (MQ=255)

                                                         gcATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAAGAGAAACTTATGCTTCGATAGCGTGCGTTAAGGCATCCATAc                                <  1:5936144/92‑1 (MQ=255)

                                                          cATGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATAcc                               <  1:6477228/93‑1 (MQ=255)

                                                          cATGGTCACGGCTTACAGTGCACAAGCGTCAGCTATCGGCGTGGCAGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATAcc                               <  1:6314020/93‑1 (MQ=255)

                                                           aTGGTCACGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATAccc                              >  1:2429266/1‑93 (MQ=255)

                                                                 aCGGCTTTCAGTGCACAAGCGTCAGCGATCGGCGCGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTgcg                        >  1:6099604/1‑93 (MQ=255)

                                                                 aCGGCTTTCAGGGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAACGCATCCATACCCGCTgcg                        >  1:3791650/1‑93 (MQ=255)

                                                                   ggCTTTCAGTGCACAAGCGTCAGCGATCTGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTgcggc                      >  1:5003517/1‑93 (MQ=255)

                                                                          aGTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGt               <  2:48816/93‑1 (MQ=255)

                                                                           gTGCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTg              >  1:5647503/1‑93 (MQ=255)

                                                                            ggCACAAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGa             <  2:2107362/92‑1 (MQ=255)

                                                                              cacaAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACt           >  2:914775/1‑93 (MQ=255)

                                                                              cacaAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACt           >  2:3943581/1‑93 (MQ=255)

                                                                                caAGCGTCAGCGTTCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTc         <  2:3228953/93‑1 (MQ=255)

                                                                                caAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTc         <  2:4836637/93‑1 (MQ=255)

                                                                                caAGCGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTc         <  2:3627430/93‑1 (MQ=255)

                                                                                   ccGTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATTGCGTGCGTTAAGCCATCCATACCCGTTGCGGCGGTCAGTGACCTCGGc      <  1:1177665/92‑1 (MQ=255)

                                                                                     gTCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGGTAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCAt    <  1:4159409/93‑1 (MQ=255)

                                                                                      tCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCGGTGACTTTGGCATa   >  2:3291892/1‑93 (MQ=255)

                                                                                      tCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATa   <  1:3718724/93‑1 (MQ=255)

                                                                                      tCAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATa   <  2:1658569/93‑1 (MQ=255)

                                                                                       cAGCGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATAc  <  2:4428932/93‑1 (MQ=255)

CTTCACGCGCTTTCGCATTACTGCCATCTTCAACGGCTAACGGCAGGTTTTCGGCAATCATGGTCACGGCTTTCAGTGCACAAGCGTCAGTGATCGGCGTGGCGGCAATAGAAACATATGCTTCGATAGCGTGCGTTAAGGCATCCATACCCGTTGCGGCGGTCAGTGACTTCGGCATAC  >  NC_000913/3754731‑3754910

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