BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	747,952	A→G	V85A (GTA→GCA)	ybgO ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	747,952	0	A	G	99.2%	418.6 / ‑4.9	131	V85A (GTA→GCA)	ybgO	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (53/77); total (54/77)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 137 total aligned reads displayed.

CGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGAT  >  NC_000913/747862‑748040

   cTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGCTTACCCAg                                                                                     >  1:4411876/1‑93 (MQ=255)

   cTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAAAGATCCTGCACCCCCGGATAGGGGTTTGTCCAg                                                                                     >  2:4229918/1‑93 (MQ=255)

    tttttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAgg                                                                                    <  1:2622542/93‑1 (MQ=255)

    tttttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAgg                                                                                    <  1:2622541/93‑1 (MQ=255)

      tttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCg                                                                                  <  1:4679355/93‑1 (MQ=255)

       ttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGt                                                                                 <  2:2477121/93‑1 (MQ=255)

       ttCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGt                                                                                 >  2:188845/1‑93 (MQ=255)

        tCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTa                                                                                <  2:590883/93‑1 (MQ=255)

         cccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAg                                                                               <  2:4391383/93‑1 (MQ=255)

         cccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAg                                                                               >  1:2127985/1‑93 (MQ=255)

          ccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGa                                                                              <  1:1330057/93‑1 (MQ=255)

          ccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGa                                                                              <  2:955205/93‑1 (MQ=255)

          ccAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGa                                                                              <  1:2970752/93‑1 (MQ=255)

           cAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGAc                                                                             <  1:3604787/93‑1 (MQ=255)

            aGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACa                                                                            <  1:2687037/93‑1 (MQ=255)

            aGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACa                                                                            <  1:310453/93‑1 (MQ=255)

            aGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACa                                                                            <  1:1727262/93‑1 (MQ=255)

            aGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACa                                                                            <  1:3303537/93‑1 (MQ=255)

               cTGCCATAGTAATCGACGCCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGc                                                                         <  1:3076285/93‑1 (MQ=255)

                 ggCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCTCGGATAGGGGTTTGCCCAGGCGTAGACATGCTc                                                                       >  1:4565446/1‑93 (MQ=255)

                  gCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCg                                                                      >  1:2127148/1‑93 (MQ=255)

                   cATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGc                                                                     <  2:3288991/93‑1 (MQ=255)

                    aTCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCt                                                                    <  2:4662770/93‑1 (MQ=255)

                    aTCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCt                                                                    >  1:2310688/1‑93 (MQ=255)

                     tCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCtt                                                                   <  2:406093/93‑1 (MQ=255)

                        tAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCa                                                                >  1:4694843/1‑93 (MQ=255)

                         aaTCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCaa                                                               >  2:4096167/1‑93 (MQ=255)

                           tCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAg                                                             >  2:714183/1‑93 (MQ=255)

                            cGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGt                                                            <  2:1056510/93‑1 (MQ=255)

                             gACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGtt                                                           <  1:2750626/93‑1 (MQ=255)

                             gACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGtt                                                           <  2:306840/93‑1 (MQ=255)

                              acacCTTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTg                                                          >  2:3729362/1‑93 (MQ=255)

                              acacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTg                                                          >  2:3247996/1‑93 (MQ=255)

                               cacCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGc                                                         >  2:2804408/1‑93 (MQ=255)

                                 ccGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCg                                                       <  2:4606494/93‑1 (MQ=255)

                                  cGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGt                                                      <  2:2780120/93‑1 (MQ=255)

                                   gTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGtt                                                     <  2:141988/93‑1 (MQ=255)

                                     tGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTgg                                                   >  1:192676/1‑93 (MQ=255)

                                     tGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTgg                                                   <  2:838113/93‑1 (MQ=255)

                                      gTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGt                                                  >  1:1202568/1‑93 (MQ=255)

                                      gTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGt                                                  <  1:1627427/93‑1 (MQ=255)

                                       tATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTa                                                 <  1:3120605/93‑1 (MQ=255)

                                        acgTCACGCCGACTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAt                                                <  2:2139491/91‑1 (MQ=255)

                                         tGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAtt                                               >  2:2595377/1‑93 (MQ=255)

                                         tGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAtt                                               >  1:1792205/1‑93 (MQ=255)

                                         tGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAtt                                               >  2:675036/1‑93 (MQ=255)

                                         tGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAtt                                               <  1:3827873/93‑1 (MQ=255)

                                            cACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTAttgtt                                            <  1:2943426/93‑1 (MQ=255)

                                               gCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTc                                         <  2:371554/93‑1 (MQ=255)

                                                ccGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTca                                        <  2:2917484/93‑1 (MQ=255)

                                                 cGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTcac                                       >  1:4740268/1‑93 (MQ=255)

                                                  gAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGCAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTcaca                                      <  2:1764143/93‑1 (MQ=255)

                                                   aGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACaa                                     <  2:4692335/93‑1 (MQ=255)

                                                    gTTGATAATAACGATTCTGCACCCCCCGATAGGGGTTTGCCCAGGCGTAGGCATGCTCTCTTTCAAGCTTGCCGTTGGTATTGTTGTCATAAt                                    >  2:4028731/1‑93 (MQ=255)

                                                    gTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAAt                                    >  2:4099335/1‑93 (MQ=255)

                                                    gTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAAt                                    >  2:3742530/1‑93 (MQ=255)

                                                     ttGATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATa                                   <  2:4376996/93‑1 (MQ=255)

                                                       gATAATAACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGa                                 >  1:2977145/1‑93 (MQ=255)

                                                        ataataACGATCCTTCACCCGCGGATAGGGGGTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGAc                                <  1:746151/93‑1 (MQ=255)

                                                        ataataACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGAc                                <  1:3987334/93‑1 (MQ=255)

                                                        ataataACGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGAc                                >  1:4193740/1‑93 (MQ=255)

                                                         taataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACg                               >  1:4588713/1‑93 (MQ=255)

                                                          aataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACgg                              <  1:2778825/93‑1 (MQ=255)

                                                          aataaCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACgg                              <  2:2949860/93‑1 (MQ=255)

                                                              aCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGCAt                          >  1:3396632/1‑93 (MQ=255)

                                                              aCGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAt                          <  2:2824275/93‑1 (MQ=255)

                                                               cGATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAtt                         <  2:1289032/93‑1 (MQ=255)

                                                                gATCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAttt                        <  2:1226584/93‑1 (MQ=255)

                                                                 aTCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAAtttt                       >  1:3752687/1‑93 (MQ=255)

                                                                  tCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTa                      >  1:1794031/1‑93 (MQ=255)

                                                                  tCCTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTa                      >  1:4229297/1‑93 (MQ=255)

                                                                   ccTGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTaa                     <  2:3167960/93‑1 (MQ=255)

                                                                    cTGCCCCCCCGGATAGGGGTTTGCACAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAat                    <  2:4344383/93‑1 (MQ=255)

                                                                     tGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAata                   >  1:2619891/1‑93 (MQ=255)

                                                                     tGCACCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAata                   >  1:594086/1‑93 (MQ=255)

                                                                     tGCACCCCCGGATAGGGGGTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAata                   >  1:3604251/1‑93 (MQ=255)

                                                                       cACCCCCGGATGGGGGTTTGCCCAGGCGTAGACCTGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTATCACCATAGACGGGAATTTTAATATc                 <  2:4221127/93‑1 (MQ=255)

                                                                        aCCCCCGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCg                <  1:3193275/93‑1 (MQ=255)

                                                                         cccccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGt               >  2:2529814/1‑93 (MQ=255)

                                                                          ccccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATTGACGGGAATTTTAATATCGtg              >  2:1102806/1‑92 (MQ=255)

                                                                          ccccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTc              <  1:4467173/93‑1 (MQ=255)

                                                                          ccccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTc              <  2:1329138/93‑1 (MQ=255)

                                                                           cccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCa             >  2:2176971/1‑93 (MQ=255)

                                                                           cccGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAACGTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCa             <  1:2731446/93‑1 (MQ=255)

                                                                             cGGATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGa           <  2:2927232/93‑1 (MQ=255)

                                                                              ggATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGAt          <  1:743903/93‑1 (MQ=255)

                                                                               gATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGAtt         <  1:2671882/93‑1 (MQ=255)

                                                                               gATAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGAtt         >  2:372450/1‑93 (MQ=255)

                                                                                 tAGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTACTATCGTCAGCTTcc       >  1:839012/1‑93 (MQ=255)

                                                                                  aGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTccc      <  2:4714597/93‑1 (MQ=255)

                                                                                  aGGGGTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTccc      <  1:741284/93‑1 (MQ=255)

                                                                                   ggggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCa     >  1:4085847/1‑93 (MQ=255)

                                                                                   ggggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCa     >  1:3766251/1‑93 (MQ=255)

                                                                                   ggggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCa     >  1:2343607/1‑93 (MQ=255)

                                                                                    gggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAg    >  2:4137754/1‑93 (MQ=255)

                                                                                    gggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAg    >  2:4617143/1‑93 (MQ=255)

                                                                                     ggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGa   >  2:32364/1‑93 (MQ=255)

                                                                                     ggTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGa   >  1:606275/1‑93 (MQ=255)

                                                                                      gTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGAt  <  2:1712050/93‑1 (MQ=255)

                                                                                      gTTTGCCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGAt  >  2:1552046/1‑93 (MQ=255)

CGACTTTTTCCCAGACTGGCATCGTAATCGACACCGTTGTATGTCACGCCGAGTTGATAATAACGATCCTGCACCCCCGGATAGGGGTTTACCCAGGCGTAGACATGCTCGCTTTCAAAGTTGCCGTTGGTATTGTTGTCACAATAGACGGGAATTTTAATATCGTCAGATTCCCAGAT  >  NC_000913/747862‑748040

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