BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,985,741	T→C	V313V (GTA→GTG)	hemX ←	uroporphyrinogen III methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,985,741	0	T	C	100.0%	690.2 / NA	222	V313V (GTA→GTG)	hemX	uroporphyrinogen III methylase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (110/112); total (110/112)

Only 100 of 234 total aligned reads displayed.

AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGC  >  NC_000913/3985651‑3985829

aaGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAgcacgcac                                                                                        >  2:2137027/1‑93 (MQ=255)

 aGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAcc                                                                                       >  1:4563198/1‑93 (MQ=255)

  gATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCAccc                                                                                      >  2:408176/1‑93 (MQ=255)

     ccATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAgg                                                                                   <  1:3162181/93‑1 (MQ=255)

      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                  <  2:4345470/93‑1 (MQ=255)

      cATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGt                                                                                  <  1:679713/93‑1 (MQ=255)

         cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                               <  1:3462859/93‑1 (MQ=255)

         cGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGga                                                                               <  2:1012832/93‑1 (MQ=255)

           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                             >  1:2824111/1‑93 (MQ=255)

           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                             >  2:3230430/1‑93 (MQ=255)

           agaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGgaga                                                                             >  1:569193/1‑93 (MQ=255)

            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                            >  1:1084162/1‑93 (MQ=255)

            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                            >  1:3237851/1‑93 (MQ=255)

            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                            <  2:2387165/93‑1 (MQ=255)

            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                            >  1:1426890/1‑93 (MQ=255)

            gaTATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGAc                                                                            <  1:1229591/93‑1 (MQ=255)

             atatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACg                                                                           >  1:1847037/1‑93 (MQ=255)

              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                          >  2:1479769/1‑93 (MQ=255)

              tatTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGt                                                                          <  1:3879338/93‑1 (MQ=255)

                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                        <  2:1346046/93‑1 (MQ=255)

                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                        <  1:3738907/93‑1 (MQ=255)

                ttttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtc                                                                        >  2:4286464/1‑93 (MQ=255)

                 tttGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTtct                                                                       <  1:2745462/93‑1 (MQ=255)

                   tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                     <  2:246257/93‑1 (MQ=255)

                   tGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTcc                                                                     <  1:966753/93‑1 (MQ=255)

                     cTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGGTCTCCAg                                                                   >  2:1717656/1‑93 (MQ=255)

                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                  >  1:2423317/1‑93 (MQ=255)

                      tGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgc                                                                  <  2:2853240/93‑1 (MQ=255)

                       ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAgcg                                                                 >  2:2613519/1‑93 (MQ=255)

                       ggCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGGGGAGACGTTCTCCAgcg                                                                 >  2:4461049/1‑93 (MQ=255)

                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                               <  2:100089/93‑1 (MQ=255)

                         cTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGcc                                                               <  1:2779546/93‑1 (MQ=255)

                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                              <  2:3604381/93‑1 (MQ=255)

                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                              >  1:4183468/1‑93 (MQ=255)

                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                              >  1:4316488/1‑93 (MQ=255)

                          ttAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCt                                                              >  1:2038242/1‑93 (MQ=255)

                             aCTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGc                                                           <  2:594146/93‑1 (MQ=255)

                              cTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCg                                                          >  2:1645982/1‑93 (MQ=255)

                                ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                        <  2:1400518/93‑1 (MQ=255)

                                ggTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGAt                                                        >  2:1250104/1‑93 (MQ=255)

                                 gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                       <  2:2899906/93‑1 (MQ=255)

                                 gTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATa                                                       >  2:4128234/1‑93 (MQ=255)

                                     cccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                   <  1:1635152/92‑1 (MQ=255)

                                     aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                   >  2:813560/1‑93 (MQ=255)

                                     aCCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtc                                                   >  2:912061/1‑93 (MQ=255)

                                      ccTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtct                                                  >  1:1160731/1‑93 (MQ=255)

                                       cTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                 <  2:4234237/93‑1 (MQ=255)

                                       cTCGTCGAGGAACGCTTTGGTGCTGGCATCATCAGTATCGTAGTAAGCGCGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGtctc                                                 <  2:4599086/93‑1 (MQ=255)

                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                >  2:521029/1‑93 (MQ=255)

                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                >  2:1147639/1‑93 (MQ=255)

                                        tcgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTctct                                                <  1:1314162/93‑1 (MQ=255)

                                         cgtcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCtt                                               >  1:4192154/1‑93 (MQ=255)

                                           tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                             <  1:2267905/93‑1 (MQ=255)

                                           tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                             <  2:298453/93‑1 (MQ=255)

                                           tcgAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTcc                                             <  2:669417/93‑1 (MQ=255)

                                             gAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                           <  2:3741581/93‑1 (MQ=255)

                                             gAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTg                                           <  2:1844453/93‑1 (MQ=255)

                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                         >  2:2789238/1‑93 (MQ=255)

                                               ggAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGt                                         <  2:3033130/93‑1 (MQ=255)

                                                gAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTg                                        >  2:1438547/1‑93 (MQ=255)

                                                 aaCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGa                                       <  1:4454896/93‑1 (MQ=255)

                                                  aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                      >  2:3676957/1‑93 (MQ=255)

                                                  aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                      >  1:3378651/1‑93 (MQ=255)

                                                  aCGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAc                                      >  1:3076506/1‑93 (MQ=255)

                                                    gCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGACACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                    <  1:2508210/93‑1 (MQ=255)

                                                    gCTTTGGTGGTGCCATCATCAGCATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCCCCTGGCGATAAGTCTCTTCCTGGTGAcgc                                    <  1:2705829/93‑1 (MQ=255)

                                                      tttGGTGGTGGCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCgg                                  >  1:101293/1‑93 (MQ=255)

                                                              tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                          <  1:4396187/93‑1 (MQ=255)

                                                              tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                          <  2:4310530/93‑1 (MQ=255)

                                                              tggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCtt                          <  2:1108053/93‑1 (MQ=255)

                                                               ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                         <  1:4504494/93‑1 (MQ=255)

                                                               ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                         >  1:581051/1‑93 (MQ=255)

                                                               ggCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCACCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtg                         <  1:3297763/93‑1 (MQ=255)

                                                                gCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                        <  2:2715064/93‑1 (MQ=255)

                                                                gCATCATCAGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTtgt                        >  2:3079258/1‑93 (MQ=255)

                                                                   tcatcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgct                     >  2:3319291/1‑93 (MQ=255)

                                                                     atcaGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGtgctgc                   <  2:1997972/93‑1 (MQ=255)

                                                                        aGTATCGTAGTAAGCACGCCCCCAGGTGGATACGGTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                <  1:2935832/93‑1 (MQ=255)

                                                                        aGTATCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGAc                >  1:297738/1‑93 (MQ=255)

                                                                          tATCGTAGTAAGCATGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACca              <  1:1563060/93‑1 (MQ=255)

                                                                          gaTCGGAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACca              <  1:196538/92‑1 (MQ=255)

                                                                           aTCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcag             >  1:800842/1‑93 (MQ=255)

                                                                           aTCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcag             <  1:2002471/93‑1 (MQ=255)

                                                                            tCGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagc            >  1:1921347/1‑93 (MQ=255)

                                                                             cgaggtaAGCACGCTCCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagca           <  2:3943184/89‑1 (MQ=255)

                                                                             cGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagca           <  2:2770977/93‑1 (MQ=255)

                                                                             cGTAGTAAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagca           <  2:1066402/93‑1 (MQ=255)

                                                                              gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag          <  1:1950324/93‑1 (MQ=255)

                                                                              gtagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACcagcag          >  1:2855576/1‑93 (MQ=255)

                                                                               tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg         >  2:4384993/1‑93 (MQ=255)

                                                                               tagtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAgg         <  1:903121/93‑1 (MQ=255)

                                                                                 gtaAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCg       >  1:4132728/1‑93 (MQ=255)

                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga      <  1:4036957/93‑1 (MQ=255)

                                                                                  taAGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCga      <  2:2210127/93‑1 (MQ=255)

                                                                                   aaGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgag     >  1:4603606/1‑93 (MQ=255)

                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga    <  2:1339390/93‑1 (MQ=255)

                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga    >  2:2363384/1‑93 (MQ=255)

                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga    <  1:2360050/93‑1 (MQ=255)

                                                                                    aGCACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCgaga    <  1:3534024/93‑1 (MQ=255)

                                                                                      cACGCACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGc  <  1:1884651/93‑1 (MQ=255)

AAGATCCATCGAGATATTTTGCTGGCTTAACTGGTCCACCTCGTCGAGGAACGCTTTGGTGGTGGCATCATCAGTATCGTAGTAAGCACGTACCCAGGTGGAGACGTTCTCCAGCGCCTGGCGATAAGTCTCTTCCTGGTGACGCGGTACAGCTTGTGCTGCGACCAGCAGGCGAGAGC  >  NC_000913/3985651‑3985829

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