BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,989,527	T→C	Q101R (CAG→CGG)	hemD ←	uroporphyrinogen‑III synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,989,527	0	T	C	100.0%	411.1 / NA	108	Q101R (CAG→CGG)	hemD	uroporphyrinogen‑III synthase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (35/73); total (35/73)

Only 100 of 109 total aligned reads displayed.

CGTCAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCTGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGCAAACA  >  NC_000913/3989391‑3989662

cGTCAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTc                                                                                                                                      >  1:189716/1‑140 (MQ=255)

cGTCAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTc                                                                                                                                      <  1:1381250/140‑1 (MQ=255)

cGTCAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTc                                                                                                                                      <  1:1310845/140‑1 (MQ=255)

   cAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCAc                                                                                                                                   <  2:1108138/140‑1 (MQ=255)

   cAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCAc                                                                                                                                   <  2:189411/140‑1 (MQ=255)

   cAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCAc                                                                                                                                   <  2:812242/140‑1 (MQ=255)

   cAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCAc                                                                                                                                   <  1:32750/140‑1 (MQ=255)

     gggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACtt                                                                                                                                 <  1:1035414/140‑1 (MQ=255)

     gggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACtt                                                                                                                                 <  1:891661/140‑1 (MQ=255)

     gggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACtt                                                                                                                                 <  1:1328341/140‑1 (MQ=255)

     gggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACtt                                                                                                                                 <  1:115336/140‑1 (MQ=255)

     gggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACtt                                                                                                                                 <  2:655804/140‑1 (MQ=255)

      ggTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTa                                                                                                                                <  1:206843/140‑1 (MQ=255)

       gTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTAc                                                                                                                               <  2:655061/140‑1 (MQ=255)

        tATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACg                                                                                                                              <  1:918841/140‑1 (MQ=255)

        tATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACg                                                                                                                              <  2:21494/140‑1 (MQ=255)

        tATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACg                                                                                                                              <  2:400458/140‑1 (MQ=255)

        tATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACg                                                                                                                              <  2:464792/140‑1 (MQ=255)

         aTCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACgg                                                                                                                             <  1:93080/140‑1 (MQ=255)

         aTCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACgg                                                                                                                             <  2:281356/140‑1 (MQ=255)

         aTCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACgg                                                                                                                             <  1:324856/140‑1 (MQ=255)

             ccAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATg                                                                                                                         <  2:1361898/140‑1 (MQ=255)

             ccAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATg                                                                                                                         <  1:336519/140‑1 (MQ=255)

                 ttAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGt                                                                                                                     <  2:496414/140‑1 (MQ=255)

                   aGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGc                                                                                                                   <  2:424339/140‑1 (MQ=255)

                   aGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGc                                                                                                                   <  1:816182/140‑1 (MQ=255)

                    gCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGcc                                                                                                                  >  1:975283/1‑140 (MQ=255)

                    gCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGcc                                                                                                                  >  1:356254/1‑140 (MQ=255)

                    gCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGcc                                                                                                                  >  1:1432483/1‑140 (MQ=255)

                     cTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCa                                                                                                                 <  2:1310556/140‑1 (MQ=255)

                     cTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCa                                                                                                                 >  2:635909/1‑140 (MQ=255)

                        aCGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTg                                                                                                              <  1:1598205/140‑1 (MQ=255)

                          gACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCg                                                                                                            >  2:50814/1‑140 (MQ=255)

                           accaccATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCgg                                                                                                           <  1:387610/140‑1 (MQ=255)

                             caccaTTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCggtg                                                                                                         <  2:146270/140‑1 (MQ=255)

                               ccaTTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCggtggt                                                                                                       <  1:684363/140‑1 (MQ=255)

                                    gCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTc                                                                                                  <  2:1451753/140‑1 (MQ=255)

                                      cACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCa                                                                                                >  2:265897/1‑140 (MQ=255)

                                      cACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCa                                                                                                >  2:1124859/1‑140 (MQ=255)

                                         gTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATg                                                                                             <  2:271600/140‑1 (MQ=255)

                                         gTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATg                                                                                             <  2:1251384/140‑1 (MQ=255)

                                          tAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATgg                                                                                            >  1:1237061/1‑140 (MQ=255)

                                           aaTATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGc                                                                                           <  1:947529/140‑1 (MQ=255)

                                            atatCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCg                                                                                          <  2:161387/140‑1 (MQ=255)

                                            atatCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCg                                                                                          <  2:528926/140‑1 (MQ=255)

                                             tatCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGa                                                                                         <  2:622854/140‑1 (MQ=255)

                                             tatCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGa                                                                                         <  1:905474/140‑1 (MQ=255)

                                                     cACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCa                                                                                 >  1:1154574/1‑140 (MQ=255)

                                                      aCGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAg                                                                                <  1:1439952/140‑1 (MQ=255)

                                                       cGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAgg                                                                               >  1:806714/1‑140 (MQ=255)

                                                       cGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAgg                                                                               <  1:723603/140‑1 (MQ=255)

                                                       cGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAgg                                                                               <  2:31150/140‑1 (MQ=255)

                                                           tGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGt                                                                           <  2:1065755/140‑1 (MQ=255)

                                                             cccGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCg                                                                         >  2:921907/1‑140 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:297263/140‑1 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:976711/140‑1 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  1:1332526/140‑1 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:373108/140‑1 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:148359/140‑1 (MQ=255)

                                                                gCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:1223711/140‑1 (MQ=255)

                                                                gCAATATTTTGGAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCgggg                                                                      <  2:1348111/140‑1 (MQ=255)

                                                                 cAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGc                                                                     <  2:331061/140‑1 (MQ=255)

                                                                      ttttGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCAttt                                                                <  1:635909/140‑1 (MQ=255)

                                                                       tttGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTa                                                               <  2:129704/140‑1 (MQ=255)

                                                                           tAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGAt                                                           <  1:872564/140‑1 (MQ=255)

                                                                                cAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGc                                                      >  2:1564026/1‑140 (MQ=255)

                                                                                cAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGc                                                      >  2:1046686/1‑140 (MQ=255)

                                                                                   gTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGa                                                   <  1:996187/140‑1 (MQ=255)

                                                                                   gTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGa                                                   >  1:1362871/1‑140 (MQ=255)

                                                                                    tAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGAt                                                  <  2:5214/140‑1 (MQ=255)

                                                                                     aaTTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATg                                                 <  1:720694/140‑1 (MQ=255)

                                                                                      aTTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGc                                                <  1:50814/140‑1 (MQ=255)

                                                                                           aGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCtg                                           >  1:1265953/1‑140 (MQ=255)

                                                                                            gCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCtgt                                          >  2:861263/1‑140 (MQ=255)

                                                                                             cAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCtgtg                                         >  2:1084540/1‑140 (MQ=255)

                                                                                                gACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGAGAtt                                      >  2:198200/1‑140 (MQ=255)

                                                                                                 aCTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTg                                     >  1:264409/1‑140 (MQ=255)

                                                                                                  cTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTgg                                    >  1:1185055/1‑140 (MQ=255)

                                                                                                  cTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTgg                                    >  2:1359408/1‑140 (MQ=255)

                                                                                                    tCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGc                                  >  1:1449708/1‑140 (MQ=255)

                                                                                                     cGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCa                                 >  1:1023259/1‑140 (MQ=255)

                                                                                                       cTGATTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCaaa                               <  2:932009/140‑1 (MQ=255)

                                                                                                          aTTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGc                            >  1:663715/1‑140 (MQ=255)

                                                                                                          aTTTCCCGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGCAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGc                            >  2:1166390/1‑140 (MQ=255)

                                                                                                              cccGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCA‑TTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACcg                       <  2:1600683/140‑1 (MQ=255)

                                                                                                               ccGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACcg                       <  2:732746/140‑1 (MQ=255)

                                                                                                                cGATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACcgc                      <  1:1213209/140‑1 (MQ=255)

                                                                                                                 gATCCTGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACcgcg                     <  1:1359408/140‑1 (MQ=255)

                                                                                                                      tGCGGGTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTg                <  2:975283/140‑1 (MQ=255)

                                                                                                                           gTAGAGAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCgaga           <  1:951398/140‑1 (MQ=255)

                                                                                                                             agagAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGTTGCGTTCAATGGCGTAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAgg         >  1:819824/1‑140 (MQ=255)

                                                                                                                             agagAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCa                           <  2:565372/122‑1 (MQ=255)

                                                                                                                             agagAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCa                           >  1:565372/1‑122 (MQ=255)

                                                                                                                             agagAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAgg         <  2:356254/140‑1 (MQ=255)

                                                                                                                             agagAATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAgg         >  1:1383292/1‑140 (MQ=255)

                                                                                                                               agaATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGc       >  1:549457/1‑140 (MQ=255)

                                                                                                                               agaATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGc       <  2:160955/140‑1 (MQ=255)

                                                                                                                                gaATCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGCa      <  1:861263/140‑1 (MQ=255)

                                                                                                                                 aaTCTTCCGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGCaa     >  1:1259654/1‑140 (MQ=255)

                                                                                                                                    cttccGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGCAaaca  <  2:1154574/140‑1 (MQ=255)

CGTCAGGGTATCCCCAATTAGCTCACGACCACCATTGCCACGTAATATCAGCGCACGTTTGCCCGCAATATTTTGTAATTCAGGTAATTGTAGCAAGACTTCGCTGATTTCCCGATCCTGCGGGTAGAGAATCTTCTGTCCACTTACGGTATGTAGTGCCAGTGCGGTGGTGCGTCCAATGGCGAAATAATCAGGTAGTCGGGGCCATTTACGATCTTGCTGATGCAGCTGTGATTGGGCAAAAGCAACCGCGTGTTGCGAGAGGGCAAACA  >  NC_000913/3989391‑3989662

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