BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,113,766	A→G	V346A (GTT→GCT)	glpK ←	glycerol kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,113,766	0	A	G	82.8%	657.1 / 140.9	308	V346A (GTT→GCT)	glpK	glycerol kinase
Reads supporting (aligned to +/- strand): ref base A (24/29); new base G (105/150); total (129/179)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 308 total aligned reads displayed.

TGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCC  >  NC_000913/4113676‑4113856

  gTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACcac                                                                                        <  2:2819377/93‑1 (MQ=255)

   ttAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCcaca                                                                                       >  2:4226689/1‑93 (MQ=255)

      gCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATac                                                                                    <  1:1780700/93‑1 (MQ=255)

       cGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATaca                                                                                   <  1:2241963/93‑1 (MQ=255)

       cGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATaca                                                                                   <  2:6012772/93‑1 (MQ=255)

        gTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATacac                                                                                  <  2:394687/93‑1 (MQ=255)

        gTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATacac                                                                                  >  2:5613163/1‑93 (MQ=255)

        gTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATacac                                                                                  >  2:6559227/1‑93 (MQ=255)

            aCCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCa                                                                              >  2:4480058/1‑93 (MQ=255)

             ccccACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCAt                                                                             <  2:6879709/93‑1 (MQ=255)

             ccccACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCAt                                                                             <  1:324684/93‑1 (MQ=255)

               ccACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTTCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTg                                                                           <  2:1480807/93‑1 (MQ=255)

               ccACGAGTCAGACCAAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTg                                                                           >  2:7489438/1‑93 (MQ=255)

                cACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTgg                                                                          >  1:4388096/1‑93 (MQ=255)

                    aGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGtt                                                                      >  1:2496434/1‑93 (MQ=255)

                    aGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGtt                                                                      <  1:5575557/93‑1 (MQ=255)

                     gTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGttt                                                                     >  1:2710097/1‑93 (MQ=255)

                     gTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGttt                                                                     >  2:5309439/1‑93 (MQ=255)

                     gTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGttt                                                                     <  1:3986571/93‑1 (MQ=255)

                     gTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGttt                                                                     <  1:5964515/93‑1 (MQ=255)

                      tCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGtttt                                                                    <  1:6284386/93‑1 (MQ=255)

                       cAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTg                                                                   <  2:4758117/93‑1 (MQ=255)

                       cAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTg                                                                   >  2:1110875/1‑93 (MQ=255)

                         gACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCa                                                                 >  2:700448/1‑93 (MQ=255)

                         gACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCa                                                                 <  2:3847125/93‑1 (MQ=255)

                          aCCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCAc                                                                <  2:4329881/93‑1 (MQ=255)

                          aCCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCAc                                                                <  1:2026512/93‑1 (MQ=255)

                            cGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACtt                                                              <  2:6293993/93‑1 (MQ=255)

                            cGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACtt                                                              <  1:6608262/93‑1 (MQ=255)

                            cGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACtt                                                              >  1:2324688/1‑93 (MQ=255)

                             gAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACttt                                                             <  1:1642592/93‑1 (MQ=255)

                                aaTCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTtggt                                                          >  1:7155090/1‑93 (MQ=255)

                                aaTCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTtggt                                                          <  1:841105/93‑1 (MQ=255)

                                 aTCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTtggtg                                                         >  2:5810825/1‑93 (MQ=255)

                                  tCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTtggtgg                                                        <  1:7365066/93‑1 (MQ=255)

                                   cGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGc                                                       <  2:2239874/93‑1 (MQ=255)

                                   cGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGc                                                       <  2:1293877/93‑1 (MQ=255)

                                     ccccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGa                                                     >  1:5843451/1‑93 (MQ=255)

                                       ccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGaaa                                                   <  1:5013399/93‑1 (MQ=255)

                                       ccGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGaaa                                                   <  1:4176931/93‑1 (MQ=255)

                                          cgcgcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAatat                                                <  2:5024617/93‑1 (MQ=255)

                                          cgcgcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAatat                                                <  1:7323883/93‑1 (MQ=255)

                                           gcgcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATAtt                                               <  1:1951997/93‑1 (MQ=255)

                                           gcgcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATAtt                                               >  1:2941357/1‑93 (MQ=255)

                                           gcgcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATAtt                                               <  1:6017569/93‑1 (MQ=255)

                                             gcATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCg                                             <  1:6492546/93‑1 (MQ=255)

                                               aTACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGa                                           <  1:3827801/93‑1 (MQ=255)

                                                tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGaa                                          >  2:1560401/1‑93 (MQ=255)

                                                tACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGaa                                          >  2:1979869/1‑93 (MQ=255)

                                                 aCGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAAt                                         >  2:2053309/1‑93 (MQ=255)

                                                 aCGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAAt                                         >  2:2402434/1‑93 (MQ=255)

                                                 aCGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAAt                                         >  2:5225928/1‑93 (MQ=255)

                                                  cGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATc                                        <  1:5663339/93‑1 (MQ=255)

                                                     gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTa                                     >  2:249387/1‑93 (MQ=255)

                                                     gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTa                                     <  2:3398581/93‑1 (MQ=255)

                                                     gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTa                                     >  1:2626851/1‑93 (MQ=255)

                                                     gTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTa                                     >  2:6142765/1‑93 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    >  1:1700246/1‑93 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    <  1:3726334/93‑1 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    <  1:2430177/93‑1 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    >  1:4260746/1‑93 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    <  1:1348648/93‑1 (MQ=255)

                                                      tCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAg                                    >  1:3030331/1‑93 (MQ=255)

                                                        ccAGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGc                                  >  2:3159842/1‑93 (MQ=255)

                                                          aGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGcgt                                <  2:7283175/93‑1 (MQ=255)

                                                          aGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGcgt                                >  1:1566466/1‑93 (MQ=255)

                                                          aGTACGGCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGcgt                                <  2:7244895/93‑1 (MQ=255)

                                                               ggCGCACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTa                           <  1:3616283/93‑1 (MQ=255)

                                                                gcgcACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTaa                          >  2:3644615/1‑93 (MQ=255)

                                                                gcgcACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTaa                          <  2:3747493/93‑1 (MQ=255)

                                                                  gcACCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATc                        <  1:4664649/93‑1 (MQ=255)

                                                                    aCCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCaa                      <  2:2592001/93‑1 (MQ=255)

                                                                    aCCCAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCaa                      >  1:5140466/1‑93 (MQ=255)

                                                                    aCCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCaa                      >  2:4842736/1‑93 (MQ=255)

                                                                     cccAGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAAc                     <  1:4430956/93‑1 (MQ=255)

                                                                        aGCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttc                  <  2:825872/93‑1 (MQ=255)

                                                                         gCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttca                 <  2:406564/93‑1 (MQ=255)

                                                                         gCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttca                 <  2:56627/93‑1 (MQ=255)

                                                                         gCCCGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttca                 >  1:6727875/1‑93 (MQ=255)

                                                                          cccGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcat                <  2:7200513/93‑1 (MQ=255)

                                                                          cccGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcat                <  1:1157064/93‑1 (MQ=255)

                                                                           ccGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatt               >  1:7090146/1‑93 (MQ=255)

                                                                           ccGGTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatt               <  1:6242488/93‑1 (MQ=255)

                                                                            cGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcattt              <  2:3883982/93‑1 (MQ=255)

                                                                            cGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcattt              >  1:4547377/1‑93 (MQ=255)

                                                                             ggTAAATGCCGGAGCCTCATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatttc             <  2:5922889/93‑1 (MQ=255)

                                                                             ggTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatttc             <  2:3536698/93‑1 (MQ=255)

                                                                              gTAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatttca            >  2:3023163/1‑93 (MQ=255)

                                                                               tAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatttcat           >  1:3399923/1‑93 (MQ=255)

                                                                               tAAATGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACttcatttcat           <  1:4892324/93‑1 (MQ=255)

                                                                                aaaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATc          <  1:6861199/93‑1 (MQ=255)

                                                                                aaaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATc          <  2:6732367/93‑1 (MQ=255)

                                                                                aaaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATc          <  2:6260010/93‑1 (MQ=255)

                                                                                aaaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATc          <  1:6984666/93‑1 (MQ=255)

                                                                                 aaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATcg         >  1:7447415/1‑93 (MQ=255)

                                                                                 aaTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATcg         >  2:1711414/1‑93 (MQ=255)

                                                                                  aTGCCGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATcgc        <  1:1120151/93‑1 (MQ=255)

                                                                                     ccGGAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCa     >  2:4689534/1‑93 (MQ=255)

                                                                                       ggAGCCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGc   >  1:6851340/1‑93 (MQ=255)

                                                                                        gAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGcc  >  2:5520266/1‑93 (MQ=255)

TGGTTAGCGTTCACCCCACGAGTCAGACCGAAAATCGCCCCGCGCGCATACGGGTCCCAGTACGGCGCACCCAGCCCGGTAAATGCCGGAACCACATACACACCATTGGTGTTTTGCACTTTGGTGGCGAAATATTCGGAATCGTAGGCGTCGTTAATCAACTTCATTTCATCGCGCAGCC  >  NC_000913/4113676‑4113856

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