BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	707,634	A→G	A26A (GCA→GCG)	ybfM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	707,634	0	A	G	99.5%	671.3 / ‑4.6	211	A26A (GCA→GCG)	ybfM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (116/94); total (116/95)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

Only 100 of 220 total aligned reads displayed.

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         >  2:4519305/1‑93 (MQ=255)

agagGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGAt                                                                                         <  1:642041/93‑1 (MQ=255)

  agGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGc                                                                                       <  1:1577369/93‑1 (MQ=255)

    gATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCg                                                                                     >  2:5049665/1‑93 (MQ=255)

    gATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCg                                                                                     <  2:1107647/93‑1 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGa                                                                                   >  1:6386167/1‑93 (MQ=255)

      ttAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGa                                                                                   <  2:2470402/93‑1 (MQ=255)

        aaCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAgg                                                                                 >  1:857618/1‑93 (MQ=255)

         aCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGc                                                                                <  1:4507506/93‑1 (MQ=255)

         aCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGc                                                                                >  1:3472699/1‑93 (MQ=255)

         aCCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCCATGCCGGAGGc                                                                                >  1:819025/1‑93 (MQ=255)

          cccATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCt                                                                               <  1:5654985/93‑1 (MQ=255)

            cATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcg                                                                             >  1:5464831/1‑93 (MQ=255)

             aTGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcgc                                                                            >  2:72859/1‑93 (MQ=255)

             aTGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcgc                                                                            >  1:633192/1‑93 (MQ=255)

             aTGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcgc                                                                            >  1:6729937/1‑93 (MQ=255)

              tGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTcccg                                                                           >  1:1705899/1‑93 (MQ=255)

               gCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGATGCTCgcgc                                                                          <  1:6610060/93‑1 (MQ=255)

                cGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGcc                                                                         <  1:3879091/93‑1 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGa                                                                       >  1:2022971/1‑93 (MQ=255)

                  tACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGa                                                                       <  1:6128912/93‑1 (MQ=255)

                   aCGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGaa                                                                      <  2:3454664/93‑1 (MQ=255)

                   aCGTTTAGTGGCAAACGCAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGATCGCGCCGaa                                                                      >  1:4211274/1‑93 (MQ=255)

                    cGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAg                                                                     >  1:1107982/1‑93 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    <  1:712996/93‑1 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    >  2:4942578/1‑93 (MQ=255)

                     gTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAgg                                                                    >  2:6711311/1‑93 (MQ=255)

                      tttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGa                                                                   >  1:1062400/1‑93 (MQ=255)

                       ttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGAt                                                                  >  1:4071467/1‑93 (MQ=255)

                       ttAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGAt                                                                  <  1:536705/93‑1 (MQ=255)

                         aGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTc                                                                <  1:279593/93‑1 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  2:2778906/93‑1 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               >  2:4621229/1‑93 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               >  1:3259592/1‑93 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  2:6495263/93‑1 (MQ=255)

                          gTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCa                                                               <  1:5044336/93‑1 (MQ=255)

                           tGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCAt                                                              >  1:2222115/1‑93 (MQ=255)

                              cAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcga                                                           <  2:4999852/93‑1 (MQ=255)

                               aaaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgac                                                          >  2:2376341/1‑93 (MQ=255)

                               aaaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgac                                                          >  1:6482197/1‑93 (MQ=255)

                                aaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgacg                                                         >  1:4447009/1‑93 (MQ=255)

                                aaCGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATcgacg                                                         >  2:4313144/1‑93 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGGAGGATTGATCGACGAtt                                                      >  2:4485776/1‑93 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      >  1:1818771/1‑93 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      <  2:5957234/93‑1 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      >  2:5807057/1‑93 (MQ=255)

                                   gtagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGAtt                                                      >  2:4444154/1‑93 (MQ=255)

                                    tagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTc                                                     >  1:349888/1‑93 (MQ=255)

                                    tagtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTc                                                     >  1:3392991/1‑93 (MQ=255)

                                     agtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCa                                                    <  2:1795115/93‑1 (MQ=255)

                                      gtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCaa                                                   <  1:2157494/93‑1 (MQ=255)

                                      gtaCGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCCACGATTCaa                                                   >  2:3390951/1‑93 (MQ=255)

                                         cgctggcgctggCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCt                                                <  2:2013977/93‑1 (MQ=255)

                                           ctggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTa                                              >  2:3234415/1‑93 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             <  1:2425531/93‑1 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             >  2:3582241/1‑93 (MQ=255)

                                            tggcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTaa                                             >  1:5212319/1‑93 (MQ=255)

                                              gcgctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGACGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAAcc                                           <  2:3152289/93‑1 (MQ=255)

                                                gctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcgg                                         <  2:6227132/93‑1 (MQ=255)

                                                gctggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcgg                                         <  2:5146204/93‑1 (MQ=255)

                                                   ggcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACcggcgg                                      >  1:3660440/1‑93 (MQ=255)

                                                    gcTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGt                                     <  2:545427/93‑1 (MQ=255)

                                                     gtATCGCGGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCCCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTa                                    <  1:4633197/92‑1 (MQ=255)

                                                     cTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTa                                    <  1:4787229/93‑1 (MQ=255)

                                                      tATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTAt                                   <  2:2977316/93‑1 (MQ=255)

                                                       aTCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGTGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATc                                  <  2:4702248/93‑1 (MQ=255)

                                                       aTCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATc                                  <  2:2281200/93‑1 (MQ=255)

                                                        tCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCt                                 >  2:40966/1‑93 (MQ=255)

                                                          gCCGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCACCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAt                               <  1:3585631/93‑1 (MQ=255)

                                                            cGGTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTa                             <  2:4226638/93‑1 (MQ=255)

                                                             ggTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTAc                            >  2:1674991/1‑93 (MQ=255)

                                                             ggTGTTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTAc                            >  1:1131492/1‑93 (MQ=255)

                                                              gtgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACt                           <  1:6500318/93‑1 (MQ=255)

                                                              gtgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACt                           >  1:2958857/1‑93 (MQ=255)

                                                               tgtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTg                          <  1:1033344/93‑1 (MQ=255)

                                                                gtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTgg                         >  1:1686994/1‑93 (MQ=255)

                                                                gtTACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTgg                         >  2:3714615/1‑93 (MQ=255)

                                                                 ttACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGc                        >  1:4586660/1‑93 (MQ=255)

                                                                  tACAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTTGCa                       >  1:5783378/1‑93 (MQ=255)

                                                                    cAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGc                     >  2:3773093/1‑93 (MQ=255)

                                                                    cAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGc                     <  2:6409042/93‑1 (MQ=255)

                                                                    cAGCAATGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGc                     >  1:2652018/1‑93 (MQ=255)

                                                                         aTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                >  2:3317094/1‑93 (MQ=255)

                                                                         aTGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGaa                >  1:5005311/1‑93 (MQ=255)

                                                                          tGTCGGGCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAc               <  2:2834038/93‑1 (MQ=255)

                                                                               ggCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGAc          <  2:2152612/93‑1 (MQ=255)

                                                                               ggCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTAATACTGGCAGCGTGAACGCGAc          >  2:5154601/1‑93 (MQ=255)

                                                                                gCTTTATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGGACGCGAcc         >  2:1388628/1‑93 (MQ=255)

                                                                                  tttATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAAAGCGACCga       >  1:1288730/1‑92 (MQ=255)

                                                                                   ttATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTa      >  1:2340346/1‑93 (MQ=255)

                                                                                    tATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTaa     >  2:5700938/1‑93 (MQ=255)

                                                                                    tATGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGGaa     >  2:2947286/1‑93 (MQ=255)

                                                                                     aTGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTaaa    >  1:3017036/1‑93 (MQ=255)

                                                                                      tGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   >  1:6401168/1‑93 (MQ=255)

                                                                                      tGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   <  1:5389549/93‑1 (MQ=255)

                                                                                      tGGCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAg   >  1:1872587/1‑93 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGCCGATTCAACCTTAACCCGCGGTATCTATTACTTGCAGCGTGAACGCGACCGTAAAGa  >  1:2894504/1‑93 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGa  <  1:5588076/93‑1 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGa  <  1:372756/93‑1 (MQ=255)

                                                                                       ggCGATGCCGGAGGCTCGCGCCGAACGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGa  <  2:541664/93‑1 (MQ=255)

AGAGGATTAACCCATGCGTACGTTTAGTGGCAAACGTAGTACGCTGGCGCTGGCTATCGCCGGTGTTACAGCAATGTCGGGCTTTATGGCAATGCCGGAGGCTCGCGCCGAAGGATTCATCGACGATTCAACCTTAACCGGCGGTATCTATTACTGGCAGCGTGAACGCGACCGTAAAGA  >  NC_000913/707544‑707723

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