BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	791,797	T→C	G405G (GGA→GGG)	modF ←	ModF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	791,797	0	T	C	99.5%	733.5 / ‑4.7	211	G405G (GGA→GGG)	modF	ModF
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (104/106); total (105/106)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.39e-01

Only 100 of 220 total aligned reads displayed.

CAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCG  >  NC_000913/791707‑791886

 aGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGcccc                                                                                        >  2:2979642/1‑93 (MQ=255)

   cGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATTTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCa                                                                                      <  2:5892800/93‑1 (MQ=255)

    ggATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                     <  1:636778/93‑1 (MQ=255)

    ggATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAg                                                                                     >  2:6469099/1‑93 (MQ=255)

     gATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAgg                                                                                    >  2:1445260/1‑93 (MQ=255)

     gATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAgg                                                                                    >  2:3011608/1‑93 (MQ=255)

      aTCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGa                                                                                   >  2:4532610/1‑93 (MQ=255)

         aaGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAg                                                                                >  2:6132493/1‑93 (MQ=255)

         aaGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAg                                                                                >  2:472731/1‑93 (MQ=255)

          aGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGa                                                                               >  1:1014744/1‑93 (MQ=255)

           gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGTGCCAGACGCTGCTGCCCCCAGGAAAGAc                                                                              >  2:6213153/1‑93 (MQ=255)

           gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGAc                                                                              <  2:579065/93‑1 (MQ=255)

           gCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGAc                                                                              <  2:512854/93‑1 (MQ=255)

             cccTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTa                                                                            >  2:11599/1‑93 (MQ=255)

                tGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATgg                                                                         <  1:1735826/93‑1 (MQ=255)

                  tAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGaa                                                                       >  2:3861298/1‑93 (MQ=255)

                  tAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGaa                                                                       <  1:2687685/93‑1 (MQ=255)

                   aGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAAc                                                                      <  2:1818946/93‑1 (MQ=255)

                   aGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAAc                                                                      >  1:2389884/1‑93 (MQ=255)

                   aGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAAc                                                                      <  2:5068936/93‑1 (MQ=255)

                     tGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACgg                                                                    <  1:4113273/93‑1 (MQ=255)

                       gTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAg                                                                  >  2:1937487/1‑93 (MQ=255)

                         tCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCg                                                                <  1:1637738/93‑1 (MQ=255)

                         tCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCg                                                                >  1:5562514/1‑93 (MQ=255)

                              gagaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGc                                                           >  2:2523461/1‑93 (MQ=255)

                               agaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGcc                                                          <  1:1490849/93‑1 (MQ=255)

                               agaATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGcc                                                          <  2:3836955/93‑1 (MQ=255)

                                  aTAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTg                                                       >  2:2388669/1‑93 (MQ=255)

                                   tAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGc                                                      >  1:1659431/1‑93 (MQ=255)

                                    aaGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCg                                                     <  2:4953650/93‑1 (MQ=255)

                                     aGCAACGTCGGATGTTTCATCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                    <  1:5766413/93‑1 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                    <  2:3603934/93‑1 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                    <  2:2130227/93‑1 (MQ=255)

                                     aGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGt                                                    <  1:4663986/93‑1 (MQ=255)

                                      gCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGtt                                                   <  1:5742998/93‑1 (MQ=255)

                                      gCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGtt                                                   >  1:3050946/1‑93 (MQ=255)

                                      gCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGtt                                                   <  1:6189630/93‑1 (MQ=255)

                                       cAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGttt                                                  <  1:3542549/93‑1 (MQ=255)

                                       cAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGttt                                                  <  1:5714806/93‑1 (MQ=255)

                                        aaCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTa                                                 >  1:4742197/1‑93 (MQ=255)

                                        aaCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTa                                                 <  2:4793683/93‑1 (MQ=255)

                                         aCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                                >  2:2588171/1‑93 (MQ=255)

                                         aCGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTAt                                                >  1:2194565/1‑93 (MQ=255)

                                          cGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATc                                               >  2:549199/1‑93 (MQ=255)

                                            tCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCaa                                             <  2:3190096/93‑1 (MQ=255)

                                             cGGATGTTTCACCAGTGCGCGGACAATCAGCGGCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCCAt                                            >  1:4949645/1‑93 (MQ=255)

                                              ggATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATg                                           >  2:749055/1‑93 (MQ=255)

                                              ggATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATg                                           >  1:5966630/1‑93 (MQ=255)

                                               gATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGc                                          >  1:945995/1‑93 (MQ=255)

                                               gATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGc                                          <  2:405184/93‑1 (MQ=255)

                                               gATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGc                                          >  1:5067141/1‑93 (MQ=255)

                                               gATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGc                                          >  2:5319530/1‑93 (MQ=255)

                                                aTGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGcc                                         <  1:2097776/93‑1 (MQ=255)

                                                   tttCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCgag                                      <  1:5896448/93‑1 (MQ=255)

                                                    ttCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCgaga                                     >  1:405667/1‑93 (MQ=255)

                                                      cACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAat                                   <  2:3676948/93‑1 (MQ=255)

                                                        ccAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAatat                                 >  1:6115609/1‑93 (MQ=255)

                                                        ccAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAatat                                 >  2:207554/1‑93 (MQ=255)

                                                        ccAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAatat                                 >  2:3151287/1‑93 (MQ=255)

                                                         cAGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATc                                >  1:3409406/1‑93 (MQ=255)

                                                          aGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATcc                               <  2:2647203/93‑1 (MQ=255)

                                                          aGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATcc                               >  2:5609617/1‑93 (MQ=255)

                                                          aGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATcc                               <  2:5818241/93‑1 (MQ=255)

                                                          aGTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATcc                               <  2:3318524/93‑1 (MQ=255)

                                                           gTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCa                              >  2:4237370/1‑93 (MQ=255)

                                                           gTGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCa                              >  1:6286168/1‑93 (MQ=255)

                                                            tGCGCGGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAg                             <  2:1557320/93‑1 (MQ=255)

                                                                cgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                         <  2:2783689/93‑1 (MQ=255)

                                                                cgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                         <  2:4894395/93‑1 (MQ=255)

                                                                cgGACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAc                         >  2:4447161/1‑93 (MQ=255)

                                                                 ggACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAct                        <  2:1317237/93‑1 (MQ=255)

                                                                 ggACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCAct                        >  2:5563826/1‑93 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                       >  1:3803488/1‑93 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                       >  1:4775567/1‑93 (MQ=255)

                                                                  gACAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGGCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                       >  2:3584324/1‑93 (MQ=255)

                                                                  gACAATAAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActg                       >  1:5363396/1‑93 (MQ=255)

                                                                   aCAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGACActgc                      >  2:1924995/1‑93 (MQ=255)

                                                                    cAATCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgct                     <  2:4427710/93‑1 (MQ=255)

                                                                     aaTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgctg                    <  2:878947/93‑1 (MQ=255)

                                                                     aaTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCActgctg                    <  1:682672/93‑1 (MQ=255)

                                                                      aTCAGCGCCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACtgctgc                   >  2:3517640/1‑93 (MQ=255)

                                                                          gcgcCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCa               >  2:398867/1‑93 (MQ=255)

                                                                          gcgcCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCa               >  1:3593430/1‑93 (MQ=255)

                                                                          gcgcCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCa               >  2:2799442/1‑93 (MQ=255)

                                                                          gccccAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCa               <  1:1761540/90‑1 (MQ=255)

                                                                           cgcCAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAg              <  1:3994861/93‑1 (MQ=255)

                                                                             ccAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATACAGCCACTGCTGCACCAGtt            >  1:2628533/1‑93 (MQ=255)

                                                                              cAGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttt           >  2:917048/1‑93 (MQ=255)

                                                                               aGACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGtttt          <  1:5687747/93‑1 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt         >  2:1708124/1‑93 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt         <  1:4404656/93‑1 (MQ=255)

                                                                                gACGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGttttt         >  1:2113678/1‑93 (MQ=255)

                                                                                  cGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc       <  2:4152776/93‑1 (MQ=255)

                                                                                  cGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc       <  1:3096011/93‑1 (MQ=255)

                                                                                  cGCTGCTGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGc       <  1:2397937/93‑1 (MQ=255)

                                                                                   gctgctGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCt      >  1:2317814/1‑93 (MQ=255)

                                                                                   gctgctGCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCt      <  2:831719/93‑1 (MQ=255)

                                                                                    ctgctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTg     >  1:3852687/1‑93 (MQ=255)

                                                                                      gctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGc   >  2:2229956/1‑93 (MQ=255)

                                                                                       ctgCCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCg  <  1:1302938/93‑1 (MQ=255)

CAGCGGATCAAGCCCCTGTAGTGGTTCATCGAGAATAAGCAACGTCGGATGTTTCACCAGTGCGCGGACAATCAGCGCCAGACGCTGCTGTCCCCAGGAAAGACTATGGAACGGAGCGTCAGCCGTGCGTTTATCAATGCCGAGAATATCCAGCCACTGCTGCACCAGTTTTTGCTGGCG  >  NC_000913/791707‑791886

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