BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,799	G→A	T650T (ACG→ACA)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,799	0	G	A	100.0%	533.0 / NA	134	T650T (ACG→ACA)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (54/80); total (54/80)

Only 100 of 134 total aligned reads displayed.

GAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062652‑3062946

      gATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCg                                                                                                                                            <  1:124638/151‑1 (MQ=255)

         cGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGa                                                                                                                                          <  1:2053831/150‑1 (MQ=255)

            cAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTg                                                                                                                                      <  1:1456911/151‑1 (MQ=255)

               cgcCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTc                                                                                                                                   >  1:101308/1‑151 (MQ=255)

                       ccGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCt                                                                                                                           >  2:1921976/1‑151 (MQ=255)

                          aTCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGa                                                                                                                         >  1:129345/1‑150 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGaa                                                                                                                                <  1:1217450/139‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCt                                                                                                                           <  2:1121614/144‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                          <  2:576105/145‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTg                                                                                                                          <  2:218207/145‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaaaa                                                                                                                      <  2:208703/149‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaa                                                                                                                        <  1:941696/147‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGaa                                                                                                                        <  1:250967/147‑1 (MQ=255)

                              cGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGa                                                                                                                         <  2:206555/146‑1 (MQ=255)

                               ggTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAAt                                                                                                                   >  1:111265/1‑151 (MQ=255)

                                        gTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGAcgcg                                                                                                          >  2:1653611/1‑151 (MQ=255)

                                            aTTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGa                                                                                                         >  1:918918/1‑148 (MQ=255)

                                            aTTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGa                                                                                                      <  1:1939391/151‑1 (MQ=255)

                                              ttAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGat                                                                                                     >  2:1026082/1‑150 (MQ=255)

                                               tAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGata                                                                                                    <  2:551610/150‑1 (MQ=255)

                                                aaGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGatat                                                                                                   <  1:241182/150‑1 (MQ=255)

                                                aaGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATc                                                                                                  >  1:1982360/1‑151 (MQ=255)

                                                aaGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATc                                                                                                  >  2:953592/1‑151 (MQ=255)

                                                aaGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATc                                                                                                  >  2:2005128/1‑151 (MQ=255)

                                                 aGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATc                                                                                                  >  2:1782269/1‑150 (MQ=255)

                                                  gCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTg                                                                                                >  1:1757062/1‑151 (MQ=255)

                                                  gCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTg                                                                                                <  2:1306518/151‑1 (MQ=255)

                                                   cccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGc                                                                                               >  1:531959/1‑151 (MQ=255)

                                                   cccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGc                                                                                               <  2:39425/151‑1 (MQ=255)

                                                   cccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGc                                                                                               <  2:151355/151‑1 (MQ=255)

                                                   cccGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGc                                                                                               >  2:620836/1‑151 (MQ=255)

                                                    ccGATGTTCTCTACTCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCg                                                                                              <  2:1496551/151‑1 (MQ=255)

                                                      gATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTg                                                                                            <  2:258485/151‑1 (MQ=255)

                                                       aTGTTCTCTACCCCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTgg                                                                                           <  1:872358/151‑1 (MQ=255)

                                                         gTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGt                                                                                          <  2:785056/150‑1 (MQ=255)

                                                         gTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTc                                                                                         >  2:2006328/1‑151 (MQ=255)

                                                            ctctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTcg                                                                                        <  2:109417/149‑1 (MQ=255)

                                                              ctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTcgcg                                                                                      <  2:225831/149‑1 (MQ=255)

                                                              ctACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTcgcg                                                                                      <  2:593717/149‑1 (MQ=255)

                                                               tACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGt                                                                                   >  2:1922962/1‑151 (MQ=255)

                                                                acacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTg                                                                                  <  1:1908987/151‑1 (MQ=255)

                                                                 cacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTgg                                                                                 <  2:1089773/151‑1 (MQ=255)

                                                                 cacCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTg                                                                                  <  2:2027390/150‑1 (MQ=255)

                                                                   ccTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGc                                                                                <  1:693571/150‑1 (MQ=255)

                                                                     tGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTc                                                                             >  1:1711997/1‑151 (MQ=255)

                                                                        agagATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCAtt                                                                          <  1:2169606/151‑1 (MQ=255)

                                                                          agaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTcc                                                                        <  1:6359/151‑1 (MQ=255)

                                                                          agaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTcc                                                                        <  2:2105846/151‑1 (MQ=255)

                                                                           gaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccg                                                                       <  1:1261978/151‑1 (MQ=255)

                                                                           gaTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccg                                                                       <  1:115786/151‑1 (MQ=255)

                                                                            aTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgc                                                                      >  1:698032/1‑151 (MQ=255)

                                                                            aTCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTccgc                                                                      <  2:379672/151‑1 (MQ=255)

                                                                                cccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCa                                                                  <  1:1601701/151‑1 (MQ=255)

                                                                                 ccgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAg                                                                 <  2:1007845/151‑1 (MQ=255)

                                                                                  cgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAgg                                                                >  2:1460893/1‑151 (MQ=255)

                                                                                  cgccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCAAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCcgccgc                                                                   <  2:433341/148‑1 (MQ=255)

                                                                                   gccTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAgg                                                                >  2:1173404/1‑150 (MQ=255)

                                                                                     cTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAtt                                                             >  1:359220/1‑151 (MQ=255)

                                                                                       ggCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGAtt                                                             <  1:153188/149‑1 (MQ=255)

                                                                                       ggCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTAc                                                           >  2:1185921/1‑151 (MQ=255)

                                                                                        gCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACg                                                          <  1:1203366/151‑1 (MQ=255)

                                                                                          cGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGcc                                                        <  1:1795247/151‑1 (MQ=255)

                                                                                          cGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGcc                                                        >  2:574986/1‑151 (MQ=255)

                                                                                          cGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGcc                                                        >  1:770613/1‑151 (MQ=255)

                                                                                           gTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCt                                                       <  1:693606/151‑1 (MQ=255)

                                                                                           gTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCt                                                       <  2:504751/151‑1 (MQ=255)

                                                                                            tAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCtt                                                      >  1:140216/1‑151 (MQ=255)

                                                                                              gAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTcc                                                    <  2:16564/151‑1 (MQ=255)

                                                                                               aaaaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCt                                                   <  2:1329491/151‑1 (MQ=255)

                                                                                                aaaCGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTg                                                  <  1:1709946/151‑1 (MQ=255)

                                                                                                  acgacgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCa                                                <  2:326681/151‑1 (MQ=255)

                                                                                                     acgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAaga                                             <  1:2116933/151‑1 (MQ=255)

                                                                                                      cgTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAag                                              >  2:35396/1‑149 (MQ=255)

                                                                                                       gTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgc                                           <  1:684608/151‑1 (MQ=255)

                                                                                                         tCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAgcg                                          <  1:1633645/150‑1 (MQ=255)

                                                                                                            cGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCgcggc                                      <  2:1426086/151‑1 (MQ=255)

                                                                                                               aGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTg                                   <  1:1922011/151‑1 (MQ=255)

                                                                                                               aGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTg                                   >  1:1359770/1‑151 (MQ=255)

                                                                                                                gTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTgg                                  <  2:648516/151‑1 (MQ=255)

                                                                                                                     cgcATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGgcggcg                             >  1:414727/1‑151 (MQ=255)

                                                                                                                     cgcATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGgcggcg                             >  2:260737/1‑151 (MQ=255)

                                                                                                                      gcATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGgcggcg                             >  2:263979/1‑150 (MQ=255)

                                                                                                                       cATCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAt                           >  1:200688/1‑151 (MQ=255)

                                                                                                                        aTCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtt                          <  2:1770335/151‑1 (MQ=255)

                                                                                                                        aTCCTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACCCGAAGATATCTGCGTGGTCGCGGGGGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAtt                          >  2:1962094/1‑151 (MQ=255)

                                                                                                                          ccTCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGAttt                         <  2:1093873/150‑1 (MQ=255)

                                                                                                                            tCACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATg                      >  2:439697/1‑151 (MQ=255)

                                                                                                                             cACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATgg                     >  2:2152335/1‑151 (MQ=255)

                                                                                                                             cACTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATg                      <  2:1741300/150‑1 (MQ=255)

                                                                                                                              aCTGGCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGt                    >  1:441082/1‑151 (MQ=255)

                                                                                                                                 ggCTGCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCa                 <  1:2045501/151‑1 (MQ=255)

                                                                                                                                     gCCGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAgg               <  2:40134/149‑1 (MQ=255)

                                                                                                                                       cGGTCATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTac            >  1:1817834/1‑150 (MQ=255)

                                                                                                                                           cATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCt        >  1:766884/1‑150 (MQ=255)

                                                                                                                                           cATAAAACACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTa       <  1:1072062/151‑1 (MQ=255)

                                                                                                                                              aaaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTCTGGTCCAGGTACACCTATg     >  2:610228/1‑150 (MQ=255)

                                                                                                                                              aaaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGc    <  1:1529226/151‑1 (MQ=255)

                                                                                                                                              aaaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGc    >  1:941635/1‑151 (MQ=255)

                                                                                                                                               aaaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCt   >  2:1316909/1‑151 (MQ=255)

                                                                                                                                                aaCACTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTc  <  1:1899015/151‑1 (MQ=255)

GAAAGTGATCGCCAGCGCCTATTCCGATCTCGGTTTCGACGTAGATTTAAGCCCGATGTTCTCTACACCTGAAGAGATCGCCCGCCTGGCCGTAGAAAACGACGTTCACGTAGTGGGCGCATCCTCACTGGCTGCCGGTCATAAAACGCTGATCCCGGAACTGGTCGAAGCGCTGAAAAAATGGGGACGCGAAGATATCTGCGTGGTCGCGGGTGGCGTCATTCCGCCGCAGGATTACGCCTTCCTGCAAGAGCGCGGCGTGGCGGCGATTTATGGTCCAGGTACACCTATGCTC  >  NC_000913/3062652‑3062946

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