BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,804,740	A→G	D117G (GAT→GGT)	pfkB →	6‑phosphofructokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,804,740	0	A	G	81.5%	467.9 / 96.7	200	D117G (GAT→GGT)	pfkB	6‑phosphofructokinase II
Reads supporting (aligned to +/- strand): ref base A (21/16); new base G (82/81); total (103/97)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 202 total aligned reads displayed.

AAAGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGAAGCT  >  NC_000913/1804649‑1804830

aaaGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGAt                                                                                           >  1:5810252/1‑93 (MQ=255)

 aaGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTg                                                                                          <  1:6249425/93‑1 (MQ=255)

  ataCTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTAAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGa                                                                                         <  1:2347434/91‑1 (MQ=255)

  aGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGa                                                                                         >  1:103230/1‑93 (MQ=255)

  aGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGa                                                                                         >  2:1547743/1‑93 (MQ=255)

     cTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGtt                                                                                      <  2:6057088/93‑1 (MQ=255)

      tGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGttt                                                                                     >  2:1450251/1‑93 (MQ=255)

      tGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGttt                                                                                     >  1:4378765/1‑93 (MQ=255)

      tGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGttt                                                                                     <  2:5112536/93‑1 (MQ=255)

      tGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGttt                                                                                     <  2:6300885/93‑1 (MQ=255)

       ggACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTc                                                                                    <  1:5104079/93‑1 (MQ=255)

       ggACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTc                                                                                    <  2:3816284/93‑1 (MQ=255)

        gACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCg                                                                                   >  1:3427765/1‑93 (MQ=255)

         aCCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGc                                                                                  >  2:1314342/1‑93 (MQ=255)

             ggCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGc                                                                              >  1:2714143/1‑93 (MQ=255)

             ggCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGc                                                                              <  2:2317230/93‑1 (MQ=255)

              gCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCt                                                                             <  1:1435614/93‑1 (MQ=255)

               cAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGGATTTAATGAAGGTGAGTTTGGCCAGCtt                                                                            >  2:1948578/1‑93 (MQ=255)

               cAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCtt                                                                            <  1:5199893/93‑1 (MQ=255)

               cAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCtt                                                                            <  2:4012748/93‑1 (MQ=255)

               cAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGAATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGGGTTTCGCCAGCtc                                                                            >  1:336487/1‑92 (MQ=255)

                aGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTg                                                                           <  2:2382656/93‑1 (MQ=255)

                 gAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGa                                                                          >  1:1902472/1‑93 (MQ=255)

                 gAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGa                                                                          >  2:3892579/1‑93 (MQ=255)

                  aaTTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGaa                                                                         >  2:4484930/1‑93 (MQ=255)

                   aTTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAag                                                                        >  1:3134149/1‑93 (MQ=255)

                    tttACACGTACATGTGGAAGCAAGCGGTGAGCATTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAaga                                                                       <  1:181554/93‑1 (MQ=255)

                    tttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAaga                                                                       >  2:3701044/1‑93 (MQ=255)

                    tttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAaga                                                                       <  2:5993827/93‑1 (MQ=255)

                    tttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAaga                                                                       <  2:2448551/93‑1 (MQ=255)

                    tttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAaga                                                                       >  1:2224755/1‑93 (MQ=255)

                     ttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAagag                                                                      <  2:824897/93‑1 (MQ=255)

                     ttACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAagag                                                                      <  2:6243690/93‑1 (MQ=255)

                      tACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGc                                                                     >  1:5465355/1‑93 (MQ=255)

                       acacGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCa                                                                    >  2:2184933/1‑93 (MQ=255)

                       acacGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCa                                                                    >  2:5035046/1‑93 (MQ=255)

                           gTACATGTTGAAGCAAGCGGTGAGCAGTAGCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGtt                                                                <  1:1462319/93‑1 (MQ=255)

                           gTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGtt                                                                >  1:6389066/1‑93 (MQ=255)

                            tACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTc                                                               <  2:5414718/93‑1 (MQ=255)

                             aCATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCt                                                              >  1:1562357/1‑93 (MQ=255)

                             aCATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCt                                                              >  1:6044690/1‑93 (MQ=255)

                                tgtgGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGa                                                           <  1:4626604/93‑1 (MQ=255)

                                 gtgGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGaa                                                          >  1:2857769/1‑93 (MQ=255)

                                 gtgGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGaa                                                          >  1:6570134/1‑93 (MQ=255)

                                   ggAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAAt                                                        <  2:1885048/93‑1 (MQ=255)

                                   ggAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAAt                                                        >  1:4696758/1‑93 (MQ=255)

                                    gAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAAtt                                                       >  1:4893840/1‑93 (MQ=255)

                                     aagcaagcGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTg                                                      >  1:5905374/1‑93 (MQ=255)

                                        caagcGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAAt                                                   >  1:2150665/1‑93 (MQ=255)

                                         aagcGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATc                                                  <  2:3513930/93‑1 (MQ=255)

                                         aagcGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATc                                                  >  2:3336690/1‑93 (MQ=255)

                                           gcGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCg                                                >  1:5249032/1‑93 (MQ=255)

                                              gTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCgggg                                             >  2:5435628/1‑93 (MQ=255)

                                               tGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGc                                            >  1:1711806/1‑93 (MQ=255)

                                                gAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGcc                                           <  1:5893322/93‑1 (MQ=255)

                                                 aGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCa                                          <  1:2108101/93‑1 (MQ=255)

                                                   cAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATc                                        >  1:2944477/1‑93 (MQ=255)

                                                     gTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCt                                      <  1:6098359/93‑1 (MQ=255)

                                                        tCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGt                                   <  2:5951507/93‑1 (MQ=255)

                                                        tCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGt                                   >  2:2684099/1‑93 (MQ=255)

                                                            tttGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATa                               >  2:2990716/1‑93 (MQ=255)

                                                            tttGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATa                               <  1:3860127/93‑1 (MQ=255)

                                                             ttGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATaa                              <  1:1661241/93‑1 (MQ=255)

                                                             ttGTTATGCCAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATaa                              >  2:5485045/1‑93 (MQ=255)

                                                              tGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAg                             >  2:4962012/1‑93 (MQ=255)

                                                                     ccAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGc                      >  2:724826/1‑93 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATGCTGGTCATAAGCGGAAGcc                     >  2:416802/1‑93 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGcc                     >  2:3145507/1‑93 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGcc                     <  2:3165747/93‑1 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGcc                     <  2:2622024/93‑1 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGcc                     >  2:691278/1‑93 (MQ=255)

                                                                      cAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGcc                     <  1:1635469/93‑1 (MQ=255)

                                                                       aGGCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCt                    >  2:2874414/1‑93 (MQ=255)

                                                                        ggCGCGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGACTg                   >  1:1319108/1‑93 (MQ=255)

                                                                         gcgcGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgc                  <  1:409853/93‑1 (MQ=255)

                                                                         gcgcGGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgc                  >  1:3486157/1‑93 (MQ=255)

                                                                           gcgGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgccg                >  2:5808118/1‑93 (MQ=255)

                                                                            cgGCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgccgc               <  2:344475/93‑1 (MQ=255)

                                                                             ggCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgccgcc              >  2:3856617/1‑93 (MQ=255)

                                                                             ggCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgccgcc              >  2:361973/1‑93 (MQ=255)

                                                                             ggCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTgccgcc              >  1:4751056/1‑93 (MQ=255)

                                                                              gCATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCa             >  2:6621353/1‑93 (MQ=255)

                                                                               cATTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGc                      >  1:127453/1‑83 (MQ=255)

                                                                               cATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAg            >  2:122264/1‑93 (MQ=255)

                                                                                aTTAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAgg           >  2:2613112/1‑93 (MQ=255)

                                                                                 ttAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgt          >  2:4401103/1‑93 (MQ=255)

                                                                                  tAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtg         <  2:1802338/93‑1 (MQ=255)

                                                                                  tAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtg         <  2:5649399/93‑1 (MQ=255)

                                                                                  tAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtg         <  2:4713812/93‑1 (MQ=255)

                                                                                  tAAATGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtg         >  2:367531/1‑93 (MQ=255)

                                                                                   aaaTGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtgt        >  2:2502988/1‑93 (MQ=255)

                                                                                   aaaTGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGgtgt        <  1:4930522/93‑1 (MQ=255)

                                                                                    aaTGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGtgtg       <  2:4685601/93‑1 (MQ=255)

                                                                                    aaTGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGtgtg       <  1:3939470/93‑1 (MQ=255)

                                                                                      tGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGaa     <  2:5855232/93‑1 (MQ=255)

                                                                                      tGAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGaa     <  1:1144453/93‑1 (MQ=255)

                                                                                      tGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGaa     <  2:4229418/93‑1 (MQ=255)

                                                                                       gAAGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGAAg    <  1:2532306/93‑1 (MQ=255)

                                                                                        aaGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGAAGc   >  2:5581445/1‑93 (MQ=255)

                                                                                        aaGGTGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGAAGc   <  1:1876835/93‑1 (MQ=255)

AAAGACTGGACCCGGCAGAATTTACACGTACATGTGGAAGCAAGCGGTGAGCAGTATCGTTTTGTTATGCCAGGCGCGGCATTAAATGAAGATGAGTTTCGCCAGCTTGAAGAGCAAGTTCTGGAAATTGAATCCGGGGCCATCCTGGTCATAAGCGGAAGCCTGCCGCCAGGTGTGAAGCT  >  NC_000913/1804649‑1804830

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