BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,238	G→T	L133I (CTA→ATA)	yjjI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,238	0	G	T	99.7%	1252.0 / ‑4.1	328	L133I (CTA→ATA)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); major base T (172/155); minor base A (1/0); total (173/155)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 329 total aligned reads displayed.

CGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGG  >  NC_000913/4614148‑4614327

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         <  2:6257730/93‑1 (MQ=255)

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         <  1:4572119/93‑1 (MQ=255)

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         >  2:695671/1‑93 (MQ=255)

cGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATaa                                                                                         >  1:5090698/1‑93 (MQ=255)

      tGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGttataattataa                                                                                   >  1:5493448/1‑93 (MQ=255)

       ggCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAAc                                                                                  <  2:2708488/93‑1 (MQ=255)

         cGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                                                <  1:5438351/93‑1 (MQ=255)

          gTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTGTTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                                               <  1:3974/93‑1 (MQ=255)

           tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                              >  2:5242978/1‑93 (MQ=255)

           tGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                                              <  1:5821845/93‑1 (MQ=255)

            gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                                             <  1:5749792/93‑1 (MQ=255)

            gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                                             <  1:2917830/93‑1 (MQ=255)

            gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                                             <  2:2831135/93‑1 (MQ=255)

            gCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                                             <  2:7454833/93‑1 (MQ=255)

             cATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACgg                                                                            <  2:5363230/93‑1 (MQ=255)

             cATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACgg                                                                            >  1:1258370/1‑93 (MQ=255)

               tAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                                          >  2:6657615/1‑93 (MQ=255)

                aaaGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTg                                                                         >  2:849729/1‑93 (MQ=255)

                  aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGca                                                                       >  2:4337735/1‑93 (MQ=255)

                   ggCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaa                                                                      <  2:2160982/93‑1 (MQ=255)

                      ggcTGGCAGGGTTCTGGCGAGGTAACGGCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                                   <  2:2556413/91‑1 (MQ=255)

                         tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                                                >  1:4223218/1‑93 (MQ=255)

                         tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                                                >  2:3336769/1‑93 (MQ=255)

                         tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGc                                                                >  2:3661853/1‑93 (MQ=255)

                           gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                                              >  2:4606539/1‑93 (MQ=255)

                           gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                                              <  1:7013994/93‑1 (MQ=255)

                             aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                            >  2:1353635/1‑93 (MQ=255)

                             aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                            >  1:3276784/1‑93 (MQ=255)

                             aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATcc                                                            >  1:4503455/1‑93 (MQ=255)

                             aGGGATCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTTTGTTTTATTTCTAACATTCGGTTGGAACAACGCATcc                                                            >  1:6869613/1‑93 (MQ=255)

                              gggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                                           >  2:88983/1‑93 (MQ=255)

                               ggTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAg                                                          <  1:6666324/93‑1 (MQ=255)

                                gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                         >  2:3263881/1‑93 (MQ=255)

                                gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                                         <  2:3491821/93‑1 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        <  1:5581471/93‑1 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        >  1:6497523/1‑93 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        >  1:3739499/1‑93 (MQ=255)

                                 ttCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                                        <  1:3060784/93‑1 (MQ=255)

                                  tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                       <  1:4742077/93‑1 (MQ=255)

                                  tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTg                                                       >  1:2268647/1‑93 (MQ=255)

                                   cTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGc                                                      <  2:6862604/93‑1 (MQ=255)

                                     gTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                    <  1:1176172/93‑1 (MQ=255)

                                     gTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                                    >  1:5220259/1‑93 (MQ=255)

                                       cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa                                                  <  2:7356823/93‑1 (MQ=255)

                                        gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                 >  2:6585022/1‑93 (MQ=255)

                                        gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAg                                                 <  1:7115286/93‑1 (MQ=255)

                                         aGGTAACGCCAGAAACGTTTTATTCGTACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                                                >  1:5294903/1‑93 (MQ=255)

                                         aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                                                <  1:6185439/93‑1 (MQ=255)

                                         aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                                                <  2:5571292/93‑1 (MQ=255)

                                          ggTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                                               >  1:6063861/1‑93 (MQ=255)

                                            tAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAg                                             >  2:5406060/1‑93 (MQ=255)

                                             aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGa                                            <  1:6384763/93‑1 (MQ=255)

                                              aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                           >  2:6474737/1‑93 (MQ=255)

                                              aCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAc                                           <  2:490263/93‑1 (MQ=255)

                                                  cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                       >  1:246902/1‑93 (MQ=255)

                                                  cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                                       >  1:3359800/1‑93 (MQ=255)

                                                   aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                                      <  2:164622/93‑1 (MQ=255)

                                                   aGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                                      >  1:4338678/1‑93 (MQ=255)

                                                    gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                                     >  1:2650657/1‑93 (MQ=255)

                                                    gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                                     <  2:451963/93‑1 (MQ=255)

                                                    gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                                     >  1:1313796/1‑93 (MQ=255)

                                                     aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                                    >  1:816419/1‑93 (MQ=255)

                                                      aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                   >  1:4954090/1‑93 (MQ=255)

                                                      aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                                   <  2:3840851/93‑1 (MQ=255)

                                                       aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                                  <  2:1494219/93‑1 (MQ=255)

                                                       aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGa                                  <  2:4215413/93‑1 (MQ=255)

                                                          ttttATTCGTACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                               <  2:6192976/93‑1 (MQ=255)

                                                          ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                               >  1:1106503/1‑93 (MQ=255)

                                                          ttttATTCGAACATCGATCTCGTCTTGTGTTAAAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt                               >  1:1762192/1‑93 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  1:2748475/93‑1 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  2:1777508/93‑1 (MQ=255)

                                                             tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc                            <  2:4069070/93‑1 (MQ=255)

                                                               ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg                          >  2:5631421/1‑93 (MQ=255)

                                                               ttCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCg                          >  2:4787770/1‑93 (MQ=255)

                                                                  gAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                       >  2:4403591/1‑93 (MQ=255)

                                                                  gAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACg                       >  2:1105284/1‑93 (MQ=255)

                                                                      aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTaa                   >  1:2412980/1‑92 (MQ=255)

                                                                      aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                   <  2:7083274/93‑1 (MQ=255)

                                                                      aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTAc                   >  1:4437954/1‑93 (MQ=255)

                                                                        cGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgt                 >  1:1581928/1‑93 (MQ=255)

                                                                          aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg               >  2:188939/1‑93 (MQ=255)

                                                                          aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtgg               >  2:2982944/1‑93 (MQ=255)

                                                                           tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt              >  2:2789418/1‑93 (MQ=255)

                                                                           tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACgtggt              >  2:7166864/1‑93 (MQ=255)

                                                                             tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg            >  2:3069326/1‑93 (MQ=255)

                                                                             tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGtggtgg            <  1:221702/93‑1 (MQ=255)

                                                                               gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTa          >  2:730006/1‑93 (MQ=255)

                                                                                 cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        <  2:915298/93‑1 (MQ=255)

                                                                                 cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        <  1:6298600/93‑1 (MQ=255)

                                                                                 cTTGTGTTATAACTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTaaa        <  2:5510728/93‑1 (MQ=255)

                                                                                  ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAg       >  2:6043887/1‑93 (MQ=255)

                                                                                   tgtgTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa      <  2:6088753/93‑1 (MQ=255)

                                                                                   tgtgTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGa      >  1:1006556/1‑93 (MQ=255)

                                                                                    gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGaa     <  1:923352/93‑1 (MQ=255)

                                                                                     tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt    <  2:1222977/93‑1 (MQ=255)

                                                                                     tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt    <  1:1491021/93‑1 (MQ=255)

                                                                                     tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAAt    >  2:3896989/1‑93 (MQ=255)

                                                                                      gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg   >  1:3823701/1‑93 (MQ=255)

                                                                                      gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATg   <  2:235650/93‑1 (MQ=255)

                                                                                       ttATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATgg  >  1:3639716/1‑93 (MQ=255)

CGATATTGGCGTGCATAAAGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTACGTGGTGGTAAAGAATGG  >  NC_000913/4614148‑4614327

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