BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,938,052	A→G	Q135R (CAG→CGG)	pykA →	pyruvate kinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,938,052	0	A	G	100.0%	924.7 / NA	224	Q135R (CAG→CGG)	pykA	pyruvate kinase II
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (113/111); total (113/111)

Only 100 of 224 total aligned reads displayed.

AGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGC  >  NC_000913/1937954‑1938149

 ggTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGtt                                                                                                >  2:2525815/1‑101 (MQ=255)

 ggTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGtt                                                                                                <  2:2525539/101‑1 (MQ=255)

  gTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTa                                                                                               >  2:2525817/1‑101 (MQ=255)

    gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGt                                                                                                 >  1:2525821/1‑97 (MQ=255)

    gAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTaaa                                                                                             >  2:2525820/1‑101 (MQ=255)

      aGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAg                                                                                           >  2:2525824/1‑101 (MQ=255)

      aGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAg                                                                                           <  1:2525513/101‑1 (MQ=255)

       ggCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGt                                                                                          >  1:2525826/1‑101 (MQ=255)

        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTa                                                                                         >  2:2525829/1‑101 (MQ=255)

        gCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTa                                                                                         <  1:2525441/101‑1 (MQ=255)

         cGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTAc                                                                                        >  2:2525831/1‑101 (MQ=255)

          gACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACt                                                                                       >  1:2525833/1‑101 (MQ=255)

           aCAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTg                                                                                      <  2:2525512/101‑1 (MQ=255)

            cAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTgg                                                                                     >  1:2525836/1‑101 (MQ=255)

             aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTa                                                                                         <  2:2525549/96‑1 (MQ=255)

             aaaGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGa                                                                                    <  1:2525311/101‑1 (MQ=255)

             aaaCAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGa                                                                                    >  2:2525838/1‑101 (MQ=255)

               aGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAg                                                                                  >  2:2525842/1‑101 (MQ=255)

                   aaaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCa                                                                              <  2:2525706/101‑1 (MQ=255)

                    aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAg                                                                             >  1:2525846/1‑101 (MQ=255)

                    aaGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAg                                                                             >  2:2525845/1‑101 (MQ=255)

                      gTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAggg                                                                           >  2:2525848/1‑101 (MQ=255)

                       tCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGc                                                                          <  2:2525606/101‑1 (MQ=255)

                        cGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGt                                                                                                 <  1:2525427/77‑1 (MQ=255)

                         ggTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCAt                                                                        <  1:2525205/101‑1 (MQ=255)

                           tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGa                                                                      >  1:2525852/1‑101 (MQ=255)

                           tATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGa                                                                      <  2:2525566/101‑1 (MQ=255)

                            aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                     <  1:2525372/101‑1 (MQ=255)

                            aTCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGaa                                                                     >  2:2525853/1‑101 (MQ=255)

                             tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                                    <  2:2525600/101‑1 (MQ=255)

                             tCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGaaa                                                                    <  1:2525415/101‑1 (MQ=255)

                              cGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAg                                                                   <  2:2525398/101‑1 (MQ=255)

                                 cTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtgt                                                                <  2:2525520/101‑1 (MQ=255)

                                 cTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtgt                                                                <  2:2525581/101‑1 (MQ=255)

                                 cTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAgtgt                                                                <  1:2525578/101‑1 (MQ=255)

                                  tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGtt                                                               >  1:2525857/1‑101 (MQ=255)

                                  tACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGtt                                                               >  2:2525859/1‑101 (MQ=255)

                                     aaaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTg                                                                                      >  2:2525862/1‑75 (MQ=255)

                                     aaaGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCAc                                                            <  1:2525563/101‑1 (MQ=255)

                                       aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCg                                                          <  1:2525616/101‑1 (MQ=255)

                                       aGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCg                                                          >  1:2525863/1‑101 (MQ=255)

                                        ggCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGa                                                         <  1:2525626/101‑1 (MQ=255)

                                        ggCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGa                                                         <  2:2525462/101‑1 (MQ=255)

                                          cctgcctgCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAg                                                       >  1:2525865/1‑101 (MQ=255)

                                             gcctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                    >  1:2525868/1‑101 (MQ=255)

                                              cctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAc                                                   <  1:2525485/101‑1 (MQ=255)

                                              cctgcTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAc                                                   >  2:2525869/1‑101 (MQ=255)

                                                tgctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCg                                                 <  1:2525537/101‑1 (MQ=255)

                                                 gctgACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGt                                                >  1:2525871/1‑101 (MQ=255)

                                                    gccgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCgg                                             >  2:2525875/3‑101 (MQ=255)

                                                    gACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCgg                                             <  1:2525343/101‑1 (MQ=255)

                                                     aCGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggt                                            <  2:2525543/101‑1 (MQ=255)

                                                      cgtcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtg                                           <  1:2525516/101‑1 (MQ=255)

                                                        tcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtggt                                         >  2:2525879/1‑101 (MQ=255)

                                                        tcgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtggt                                         <  1:2525505/101‑1 (MQ=255)

                                                         cgtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTc                                        >  1:2525881/1‑101 (MQ=255)

                                                          gtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCAc                                                   >  2:2525884/1‑89 (MQ=255)

                                                          gtGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTcc                                       >  1:2525883/1‑101 (MQ=255)

                                                           tGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTccc                                      <  2:2525523/101‑1 (MQ=255)

                                                            gCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCg                                                 >  1:2525886/1‑89 (MQ=255)

                                                             ccTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCa                                                    <  2:2525529/85‑1 (MQ=255)

                                                              cTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctc                                   >  1:2525889/1‑101 (MQ=255)

                                                               tGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                  <  2:2525647/101‑1 (MQ=255)

                                                               tGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctct                                  <  2:2525490/101‑1 (MQ=255)

                                                                ggTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCtctc                                 <  1:2525495/101‑1 (MQ=255)

                                                                 gTGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCctc                                   <  1:2525643/98‑1 (MQ=255)

                                                                  tGACATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCca                               >  2:2525892/1‑101 (MQ=255)

                                                                     cATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaaca                            >  2:2525894/1‑101 (MQ=255)

                                                                     cATCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCcaaca                            <  1:2525391/101‑1 (MQ=255)

                                                                       tCCTGCTGCTGGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAAt                          >  1:2525896/1‑101 (MQ=255)

                                                                           gctgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                      <  2:2525653/101‑1 (MQ=255)

                                                                           gctgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAg                      >  2:2525899/1‑101 (MQ=255)

                                                                            ctgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                     >  1:2525903/1‑101 (MQ=255)

                                                                            ctgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                     >  1:2525900/1‑101 (MQ=255)

                                                                            ctgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAgg                     <  1:2525547/101‑1 (MQ=255)

                                                                             tgctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGt                    >  2:2525904/1‑101 (MQ=255)

                                                                              gctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTa                   <  2:2525388/101‑1 (MQ=255)

                                                                              gctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTa                   <  1:2525592/101‑1 (MQ=255)

                                                                               ctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTAt                  >  2:2525908/1‑101 (MQ=255)

                                                                               ctgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTAt                  >  1:2525907/1‑101 (MQ=255)

                                                                                tgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACACTAAAGGTATc                 <  1:2525685/101‑1 (MQ=255)

                                                                                tgGACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGCGGTCCCCTCTCCAACAATAAAGGTAt                  >  2:2525910/1‑100 (MQ=255)

                                                                                 ggACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCCCCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                <  1:2525555/101‑1 (MQ=255)

                                                                                 ggACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATca                <  1:2525760/101‑1 (MQ=255)

                                                                                  gACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa               <  2:2525737/101‑1 (MQ=255)

                                                                                  gACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa               >  2:2525914/1‑101 (MQ=255)

                                                                                  gACGATGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATcaa               <  1:2525631/101‑1 (MQ=255)

                                                                                      aTGGTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACaaa           >  1:2525917/1‑101 (MQ=255)

                                                                                        ggTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt         <  2:2525668/101‑1 (MQ=255)

                                                                                        ggTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACt         >  1:2525919/1‑101 (MQ=255)

                                                                                        ggTCGCGTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTATCCAACAATAAAGGTATCAACAAACt         >  2:2525920/1‑101 (MQ=255)

                                                                                           cgcgTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTgg      <  2:2525456/101‑1 (MQ=255)

                                                                                           cgcgTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTgg      <  2:2525635/101‑1 (MQ=255)

                                                                                            gcgTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggc     <  1:2525597/101‑1 (MQ=255)

                                                                                            gcgTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggc     >  2:2525923/1‑101 (MQ=255)

                                                                                             cgTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggcg    <  1:2525532/101‑1 (MQ=255)

                                                                                              gTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACtt        >  1:2525929/1‑96 (MQ=255)

                                                                                              gTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggcgg   >  1:2525924/1‑101 (MQ=255)

                                                                                              gTCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTggcgg   >  2:2525926/1‑101 (MQ=255)

                                                                                               tCCGGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCggtgg                                          <  2:2525444/61‑1 (MQ=255)

AGGTGAAGGCGACAAAGAAAAAGTCGGTATCGACTACAAAGGCCTGCCTGCTGACGTCGTGCCTGGTGACATCCTGCTGCTGGACGATGGTCGCGTCCAGTTAAAAGTACTGGAAGTTCAGGGCATGAAAGTGTTCACCGAAGTCACCGTCGGTGGTCCCCTCTCCAACAATAAAGGTATCAACAAACTTGGCGGC  >  NC_000913/1937954‑1938149

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