BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,079,836	T→C	V600V (GTA→GTG)	tktA ←	transketolase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,079,836	0	T	C	100.0%	1386.2 / NA	329	V600V (GTA→GTG)	tktA	transketolase 1
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (136/193); total (136/193)

Only 100 of 330 total aligned reads displayed.

GTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGT  >  NC_000913/3079738‑3079933

gtCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCAc                                                                                                 <  1:5188487/101‑1 (MQ=255)

gtCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCAc                                                                                                 <  2:5189017/101‑1 (MQ=255)

gtCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCAc                                                                                                 <  1:5188954/101‑1 (MQ=255)

 tCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACg                                                                                                <  2:5188854/101‑1 (MQ=255)

 tCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACg                                                                                                >  1:5189318/1‑101 (MQ=255)

  cATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACgg                                                                                               >  2:5189319/1‑101 (MQ=255)

  cATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACgg                                                                                               <  2:5188603/101‑1 (MQ=255)

  cATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACgg                                                                                               <  1:5188762/101‑1 (MQ=255)

   aTACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGa                                                                                              <  2:5188881/101‑1 (MQ=255)

   aTACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGa                                                                                              <  1:5189211/101‑1 (MQ=255)

    tACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGAt                                                                                             <  1:5188861/101‑1 (MQ=255)

    tACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGAt                                                                                             <  1:5188016/101‑1 (MQ=255)

      ccGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTc                                                                                           <  1:5188675/101‑1 (MQ=255)

       cgacgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCa                                                                                          <  1:5188821/101‑1 (MQ=255)

       cgacgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCa                                                                                          >  2:5189326/1‑101 (MQ=255)

       cgacgaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCa                                                                                          >  1:5189327/1‑101 (MQ=255)

           gaTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGt                                                                                      >  1:5189332/1‑101 (MQ=255)

            aTAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTa                                                                                     <  1:5188584/101‑1 (MQ=255)

             tAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTaa                                                                                    <  2:5188832/101‑1 (MQ=255)

             tAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTaa                                                                                    <  2:5188892/101‑1 (MQ=255)

              aGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAag                                                                                   >  1:5189336/1‑101 (MQ=255)

                cACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAagca                                                                                 >  2:5189337/1‑101 (MQ=255)

                 aCCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAagcag                                                                                <  1:5189064/101‑1 (MQ=255)

                    gTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCAt                                                                             <  2:5188652/101‑1 (MQ=255)

                    gTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGAAt                                                                             >  2:5189341/1‑101 (MQ=255)

                      tCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATcc                                                                           <  2:5188685/101‑1 (MQ=255)

                        aGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTg                                                                         <  1:5188709/101‑1 (MQ=255)

                         ggCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGc                                                                        <  2:5188948/101‑1 (MQ=255)

                         ggCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGc                                                                        <  1:5188830/101‑1 (MQ=255)

                          gCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCt                                                                       >  1:5189345/1‑101 (MQ=255)

                          gCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCt                                                                       >  2:5189347/1‑101 (MQ=255)

                            cAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                     <  1:5188803/101‑1 (MQ=255)

                            cAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTg                                                                     >  1:5189353/1‑101 (MQ=255)

                             aaCATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGt                                                                    <  2:5188503/101‑1 (MQ=255)

                              aCATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTc                                                                   >  1:5189356/1‑101 (MQ=255)

                               cATACTTGTACCAGTAGTCAGCAATACCCGCTTCTAGAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCa                                                                  <  2:5189026/101‑1 (MQ=255)

                               cATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCa                                                                  <  1:5188701/101‑1 (MQ=255)

                               cATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCa                                                                  <  2:5188558/101‑1 (MQ=255)

                                aTACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCaa                                                                 >  2:5189358/1‑101 (MQ=255)

                                 tACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCaaa                                                                >  1:5189360/1‑101 (MQ=255)

                                 tACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCaaa                                                                <  2:5188761/101‑1 (MQ=255)

                                 tACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCaaa                                                                <  2:5189136/101‑1 (MQ=255)

                                  aCTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                               <  1:5188999/101‑1 (MQ=255)

                                  aCTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAAt                                                               <  1:5188741/101‑1 (MQ=255)

                                   cTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATg                                                              <  1:5188985/101‑1 (MQ=255)

                                    ttGTACCAGTAGTCAGCAATCCCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGc                                                             <  2:5189067/101‑1 (MQ=255)

                                    ttGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGc                                                             <  1:5188872/101‑1 (MQ=255)

                                     tGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCg                                                            >  1:5189365/1‑101 (MQ=255)

                                      gTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGt                                                           <  2:5188689/101‑1 (MQ=255)

                                       tACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTc                                                          >  1:5189369/1‑101 (MQ=255)

                                       tACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTc                                                          <  1:5188944/101‑1 (MQ=255)

                                        aCCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCg                                                         >  2:5189370/1‑101 (MQ=255)

                                        aCCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCg                                                         <  2:5188906/101‑1 (MQ=255)

                                          cAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATg                                                              <  1:5188655/94‑1 (MQ=255)

                                          cAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATTCGTCGGt                                                       >  1:5189373/1‑101 (MQ=255)

                                             gagtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGa                                                    >  1:5189379/2‑101 (MQ=255)

                                              agtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGAc                                                   >  2:5189380/1‑101 (MQ=255)

                                               gtCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACg                                                  <  2:5189003/101‑1 (MQ=255)

                                                tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACgg                                                 <  2:5189172/101‑1 (MQ=255)

                                                tCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACgg                                                 <  2:5188810/101‑1 (MQ=255)

                                                 cAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGc                                                <  2:5188963/101‑1 (MQ=255)

                                                  aGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCa                                               >  1:5189383/1‑101 (MQ=255)

                                                   gCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCAt                                              <  1:5189025/101‑1 (MQ=255)

                                                     aaTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATgg                                            >  1:5189387/1‑101 (MQ=255)

                                                      aTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                           <  1:5188913/101‑1 (MQ=255)

                                                      aTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                           <  2:5188820/101‑1 (MQ=255)

                                                      aTACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGa                                           >  2:5189391/1‑101 (MQ=255)

                                                        aCCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGaca                                         >  1:5189394/1‑101 (MQ=255)

                                                         cccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAag                                                                                   <  1:5188846/58‑1 (MQ=255)

                                                          ccGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcacc                                       >  2:5189396/1‑101 (MQ=255)

                                                           cGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAgcagca                                                                              >  2:5189399/1‑61 (MQ=255)

                                                            gCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGAcaccac                                     <  2:5188924/101‑1 (MQ=255)

                                                             cTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcg                                    >  2:5189405/1‑101 (MQ=255)

                                                              ttCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACg                                                  <  2:5188743/86‑1 (MQ=255)

                                                              ttCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgc                                   >  1:5189407/1‑101 (MQ=255)

                                                               tCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGAt                                                                                             <  1:5188728/42‑1 (MQ=255)

                                                               tCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcg                                  <  1:5189142/101‑1 (MQ=255)

                                                                cTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCAcgcgc                                 <  1:5188788/101‑1 (MQ=255)

                                                                  aCAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACgcgcgc                               >  1:5189411/1‑101 (MQ=255)

                                                                  aCAGCAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACgcgcgc                               <  2:5188680/101‑1 (MQ=255)

                                                                      cAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                           >  1:5189416/1‑101 (MQ=255)

                                                                      cAACGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTc                           >  2:5189415/1‑101 (MQ=255)

                                                                        aCGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCAc                         >  2:5189424/1‑101 (MQ=255)

                                                                        aCGCGTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCAc                         >  1:5189421/1‑101 (MQ=255)

                                                                          gcgTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGc                       <  1:5188879/101‑1 (MQ=255)

                                                                            gTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                     <  2:5188901/101‑1 (MQ=255)

                                                                            gTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                     >  2:5189430/1‑101 (MQ=255)

                                                                            gTGCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCt                     >  1:5189428/1‑101 (MQ=255)

                                                                              gCAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTc                   >  1:5189435/1‑101 (MQ=255)

                                                                               cAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                  <  2:5189095/101‑1 (MQ=255)

                                                                               cAGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCg                  >  2:5189436/1‑101 (MQ=255)

                                                                                aGTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCgg                 <  2:5188794/101‑1 (MQ=255)

                                                                                 gTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                <  2:5189079/101‑1 (MQ=255)

                                                                                 gTAACCGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGc                >  2:5189441/1‑101 (MQ=255)

                                                                                      cGCTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtca           <  2:5188703/101‑1 (MQ=255)

                                                                                        cTTTCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGcagtcagt         >  1:5189443/1‑101 (MQ=255)

                                                                                         tttCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGTCACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt        <  1:5189057/101‑1 (MQ=255)

                                                                                         tttCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtt        >  2:5189448/1‑101 (MQ=255)

                                                                                          ttCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGttt       >  1:5189451/1‑101 (MQ=255)

                                                                                           tCGGCAGCACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGtttt      <  1:5188930/101‑1 (MQ=255)

GTCATACCGACGATAGCACCGTTCAGGCCAACATACTTGTACCAGTAGTCAGCAATACCCGCTTCTACAGCAACGCGTGCAGTAACCGCTTTCGGCAGTACGGATTCACGGTAAGCAGCATCCTGCTTGTCAAATGCGTCGGTAGACGGCATGGACACCACGCGCGCTTTCACGCCTTCGGCAGTCAGTTTTTCGT  >  NC_000913/3079738‑3079933

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