BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,621,342	C→T	G93G (GGC→GGT)	ppk →	polyphosphate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,621,342	0	C	T	87.5%	695.8 / 66.3	239	G93G (GGC→GGT)	ppk	polyphosphate kinase
Reads supporting (aligned to +/- strand): ref base C (18/12); new base T (103/106); total (121/118)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.55e-01

Only 100 of 247 total aligned reads displayed.

ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG  >  NC_000913/2621250‑2621434

attaGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGc                                                                                              >  1:5612695/1‑93 (MQ=255)

  taGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTct                                                                                            <  1:3529133/93‑1 (MQ=255)

   aGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGtctc                                                                                           >  2:2678307/1‑93 (MQ=255)

    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAACATTCAGTCCCGGGTGCTGAAAGCCGATCAGGTATTCGACGGTctct                                                                                          >  1:1088832/1‑93 (MQ=255)

    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTctct                                                                                          >  2:2820286/1‑93 (MQ=255)

    gCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTctct                                                                                          <  1:4731522/93‑1 (MQ=255)

     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         <  1:114724/93‑1 (MQ=255)

     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         >  2:1160878/1‑93 (MQ=255)

     cGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTa                                                                                         <  1:570987/93‑1 (MQ=255)

      gaagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTAc                                                                                        >  2:2512900/1‑93 (MQ=255)

       aagaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACa                                                                                       >  2:4684322/1‑93 (MQ=255)

        agaaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACaa                                                                                      <  2:3473227/93‑1 (MQ=255)

         gaaCAAGGCTCCAACTCTCATTCCCGCCATGTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAAc                                                                                     <  2:787284/93‑1 (MQ=255)

          aaCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACg                                                                                    <  2:5064025/93‑1 (MQ=255)

           aCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGa                                                                                   <  1:1968440/93‑1 (MQ=255)

           aCAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGa                                                                                   >  1:4488442/1‑93 (MQ=255)

              aGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCt                                                                                >  1:3070862/1‑93 (MQ=255)

              aGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCt                                                                                >  2:2016198/1‑93 (MQ=255)

               ggCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTa                                                                               <  1:4777943/93‑1 (MQ=255)

               agCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTa                                                                               <  2:1427808/92‑1 (MQ=255)

                gCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTAt                                                                              <  1:579255/93‑1 (MQ=255)

                  tCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTg                                                                            >  2:1910257/1‑93 (MQ=255)

                   ccAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGc                                                                           >  1:1119110/1‑93 (MQ=255)

                     aaCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTg                                                                         >  2:3074675/1‑93 (MQ=255)

                      aCTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTgg                                                                        >  1:5661545/1‑93 (MQ=255)

                       ctctCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGga                                                                       <  1:556701/93‑1 (MQ=255)

                         ctcATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGgaga                                                                     >  2:2421556/1‑93 (MQ=255)

                          tcATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGAt                                                                    <  1:6245852/93‑1 (MQ=255)

                           cATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATg                                                                   <  2:4048433/93‑1 (MQ=255)

                            aTTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATgg                                                                  <  2:3427137/93‑1 (MQ=255)

                            aTTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATgg                                                                  >  2:45894/1‑93 (MQ=255)

                             ttCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgc                                                                 >  2:5187668/1‑93 (MQ=255)

                              tCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGgcg                                                                >  2:42688/1‑93 (MQ=255)

                              gcccGCCATTGACTGGTCAAAAGTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcg                                                                <  1:186203/92‑1 (MQ=255)

                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              >  1:6063185/1‑93 (MQ=255)

                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              >  2:2962397/1‑93 (MQ=255)

                                ccGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGgcgcg                                                              <  1:897583/93‑1 (MQ=255)

                                  gCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCa                                                            >  1:2223731/1‑93 (MQ=255)

                                   ccATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCaa                                                           >  2:3153939/1‑93 (MQ=255)

                                    cATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAc                                                          >  1:5217308/1‑93 (MQ=255)

                                    cATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAc                                                          <  1:3584997/93‑1 (MQ=255)

                                     aTTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAcc                                                         <  1:6447721/93‑1 (MQ=255)

                                     aTTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAAcc                                                         >  2:3632241/1‑93 (MQ=255)

                                      tttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCa                                                        >  1:2360715/1‑93 (MQ=255)

                                      tttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCa                                                        >  2:2581633/1‑93 (MQ=255)

                                       ttACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAg                                                       <  1:4631330/93‑1 (MQ=255)

                                        tACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGa                                                      >  1:651319/1‑93 (MQ=255)

                                         aCTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGAt                                                     <  1:497153/93‑1 (MQ=255)

                                            gggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCtt                                                  >  2:102017/1‑93 (MQ=255)

                                            gggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCtt                                                  <  2:3205628/93‑1 (MQ=255)

                                             ggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                 <  1:5870543/93‑1 (MQ=255)

                                             ggCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTc                                                 <  2:1658880/93‑1 (MQ=255)

                                               cAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCt                                               <  1:5933950/93‑1 (MQ=255)

                                               cAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCt                                               <  1:134994/93‑1 (MQ=255)

                                                  aaTTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGCCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGAt                                            >  2:5560725/1‑93 (MQ=255)

                                                    ttCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTa                                          >  1:2922992/1‑93 (MQ=255)

                                                     tCAGTCCCGGTTGCTAAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTaa                                         <  2:6101394/93‑1 (MQ=255)

                                                      cAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAAt                                        >  2:5284541/1‑93 (MQ=255)

                                                        gTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGa                                      <  2:4455825/93‑1 (MQ=255)

                                                        gTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGa                                      <  2:6256920/93‑1 (MQ=255)

                                                          cccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAAc                                    <  2:5375144/93‑1 (MQ=255)

                                                           ccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACg                                   >  1:1611026/1‑93 (MQ=255)

                                                           ccGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACg                                   >  2:5954324/1‑93 (MQ=255)

                                                            cGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                  >  1:6200863/1‑93 (MQ=255)

                                                            cGGGTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGc                                  <  1:3875261/93‑1 (MQ=255)

                                                             gggTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGcc                                 >  1:3946554/1‑93 (MQ=255)

                                                              ggTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCa                                >  1:538977/1‑93 (MQ=255)

                                                               gTGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAg                               <  2:3366874/93‑1 (MQ=255)

                                                                tGCTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGc                              >  2:728008/1‑93 (MQ=255)

                                                                  cTGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctc                            >  2:1335331/1‑93 (MQ=255)

                                                                   tGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctct                           >  2:1949526/1‑93 (MQ=255)

                                                                   tGAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGctct                           >  1:810787/1‑93 (MQ=255)

                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          >  1:3763073/1‑93 (MQ=255)

                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          >  1:6356130/1‑93 (MQ=255)

                                                                    gAAAGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCtctc                          <  1:3661365/93‑1 (MQ=255)

                                                                     aaaGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTcc                         <  2:6318735/93‑1 (MQ=255)

                                                                     aaaGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTcc                         <  2:4322485/93‑1 (MQ=255)

                                                                       aGCCGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGt                       >  1:2501334/1‑93 (MQ=255)

                                                                       aGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCCGCTCTCCGt                       >  1:4015648/1‑93 (MQ=255)

                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  1:2992860/93‑1 (MQ=255)

                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  1:4289726/93‑1 (MQ=255)

                                                                          cGATCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaa                    <  2:569528/93‑1 (MQ=255)

                                                                            aTCAGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGtcaatc                  <  2:2307150/93‑1 (MQ=255)

                                                                               aGGAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaac               >  2:2199570/1‑93 (MQ=255)

                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              >  2:4893309/1‑93 (MQ=255)

                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              <  1:4867024/93‑1 (MQ=255)

                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              >  2:6500772/1‑93 (MQ=255)

                                                                                ggAATTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATcaaca              <  2:6591122/93‑1 (MQ=255)

                                                                                   aTTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACaaaa           >  1:3500606/1‑93 (MQ=255)

                                                                                   aTTCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACaaaa           <  1:1647344/93‑1 (MQ=255)

                                                                                    ttCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAc          <  1:6587756/93‑1 (MQ=255)

                                                                                    ttCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAAc          <  2:1603982/93‑1 (MQ=255)

                                                                                     tCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACt         <  1:4565954/93‑1 (MQ=255)

                                                                                     tCGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACt         <  1:2612494/93‑1 (MQ=255)

                                                                                      cGACGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTg        <  2:2102601/93‑1 (MQ=255)

                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:14702/93‑1 (MQ=255)

                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:3300903/93‑1 (MQ=255)

                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      >  1:3156469/1‑93 (MQ=255)

                                                                                        aCGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGc      <  1:69875/93‑1 (MQ=255)

                                                                                         cGGTCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCt     <  1:4249818/93‑1 (MQ=255)

ATTAGCGAAGAACAAGGCTCCAACTCTCATTCCCGCCATTTACTGGGCAAAATTCAGTCCCGGGTGCTGAAAGCCGATCAGGAATTCGACGGCCTCTACAACGAGCTATTGCTGGAGATGGCGCGCAACCAGATCTTCCTGATTAATGAACGCCAGCTCTCCGTCAATCAACAAAACTGGCTGCG  >  NC_000913/2621250‑2621434

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