BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,986,916	C→T	A176T (GCC→ACC)	hemD ←	uroporphyrinogen III synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,986,916	0	C	T	100.0%	919.5 / NA	252	A176T (GCC→ACC)	hemD	uroporphyrinogen III synthase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (135/117); total (135/117)

Only 100 of 262 total aligned reads displayed.

ACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACCGCGCGCCG  >  NC_000913/3986826‑3987005

aCCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGttt                                                                                         <  1:1234460/93‑1 (MQ=255)

aCCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGttt                                                                                         <  1:6165971/93‑1 (MQ=255)

aCCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGttt                                                                                         >  2:4595355/1‑93 (MQ=255)

   aGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGcc                                                                                      <  1:3559691/93‑1 (MQ=255)

   aGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGcc                                                                                      <  1:3944642/93‑1 (MQ=255)

    gTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCa                                                                                     <  1:4709060/93‑1 (MQ=255)

     tGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAg                                                                                    <  1:3783459/93‑1 (MQ=255)

     tGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAg                                                                                    >  1:6044064/1‑93 (MQ=255)

       cTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAgcg                                                                                  >  2:2140715/1‑93 (MQ=255)

        tCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAgcgc                                                                                 >  2:1242039/1‑93 (MQ=255)

        tCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAgcgc                                                                                 <  1:5901401/93‑1 (MQ=255)

        tCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAgcgc                                                                                 >  1:4865728/1‑93 (MQ=255)

         cACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCa                                                                                >  1:5422487/1‑93 (MQ=255)

         cACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCa                                                                                >  2:105143/1‑93 (MQ=255)

             aTACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGc                                                                            >  2:906376/1‑93 (MQ=255)

             aTACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGc                                                                            <  2:1677587/93‑1 (MQ=255)

             aTACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGc                                                                            >  1:2898801/1‑93 (MQ=255)

             aTACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGc                                                                            >  1:2699849/1‑93 (MQ=255)

              tACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGct                                                                           >  1:1599871/1‑93 (MQ=255)

               aCCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGctt                                                                          <  2:3568773/93‑1 (MQ=255)

                ccATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGcttc                                                                         <  2:1994977/93‑1 (MQ=255)

                ccATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGcttc                                                                         <  2:6148528/93‑1 (MQ=255)

                  aTTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGcttctt                                                                       >  2:4889922/1‑93 (MQ=255)

                  aTTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGcttctt                                                                       >  1:4113731/1‑93 (MQ=255)

                   ttGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCttcttc                                                                      <  2:1584928/93‑1 (MQ=255)

                   ttGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCttcttc                                                                      <  1:3710687/93‑1 (MQ=255)

                    tgtgGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTtcttct                                                                     <  1:3139466/93‑1 (MQ=255)

                        ggATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCAc                                                                 >  2:5689649/1‑93 (MQ=255)

                        ggATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCAc                                                                 <  1:192365/93‑1 (MQ=255)

                          aTCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCa                                                               <  2:5576369/93‑1 (MQ=255)

                          aTCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCa                                                               <  1:1197049/93‑1 (MQ=255)

                             aGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCg                                                            >  2:1366668/1‑93 (MQ=255)

                             aGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCg                                                            >  2:5223340/1‑93 (MQ=255)

                              gCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGt                                                           <  2:4989424/93‑1 (MQ=255)

                                  ccAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATg                                                       <  2:4821998/93‑1 (MQ=255)

                                  ccAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATg                                                       <  1:1803917/93‑1 (MQ=255)

                                  ccAGAGTAGCTGCAACATTTCACCGCTGATACCATCGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATg                                                       <  2:4147216/93‑1 (MQ=255)

                                   cAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATgg                                                      >  2:384969/1‑93 (MQ=255)

                                    agagTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGa                                                     <  1:856715/93‑1 (MQ=255)

                                    agagTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGa                                                     <  1:2218275/93‑1 (MQ=255)

                                    agagTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGAGCATCGCTTCTTCTGCACCATCGTAATGGa                                                     >  1:1428416/1‑93 (MQ=255)

                                     gagTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGAt                                                    >  1:713237/1‑93 (MQ=255)

                                     gagTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGAt                                                    >  1:1311661/1‑93 (MQ=255)

                                      agTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGAtt                                                   >  2:176974/1‑93 (MQ=255)

                                      agTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGAtt                                                   >  1:1454759/1‑93 (MQ=255)

                                        ttGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTgc                                                 <  2:1289715/93‑1 (MQ=255)

                                        ttGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTgc                                                 >  2:3061757/1‑93 (MQ=255)

                                            tgcAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCAt                                             >  2:729915/1‑93 (MQ=255)

                                             gcAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATc                                            >  1:6002833/1‑93 (MQ=255)

                                             gcAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATc                                            >  2:6173577/1‑93 (MQ=255)

                                             gcAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATc                                            <  1:2109610/93‑1 (MQ=255)

                                             gcAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATc                                            >  1:3736137/1‑93 (MQ=255)

                                               aaCATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGt                                          >  2:1730724/1‑93 (MQ=255)

                                               aaCATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGt                                          <  1:4533766/93‑1 (MQ=255)

                                                aCATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGtt                                         >  2:646098/1‑93 (MQ=255)

                                                aCATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGtt                                         >  1:68184/1‑93 (MQ=255)

                                                aCATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGtt                                         >  2:5419442/1‑93 (MQ=255)

                                                 cATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTg                                        <  2:2901176/93‑1 (MQ=255)

                                                  aTTTCCCCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGa                                       <  2:796903/93‑1 (MQ=255)

                                                   tttCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGAt                                      <  2:6380019/93‑1 (MQ=255)

                                                   tttCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGAt                                      <  1:322589/93‑1 (MQ=255)

                                                   tttCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGAt                                      >  2:3946221/1‑93 (MQ=255)

                                                      cACCGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAAc                                   <  1:1703001/93‑1 (MQ=255)

                                                       acagCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCTCATCGTTGATAACa                                  >  1:6156949/4‑93 (MQ=255)

                                                         cGCTGGTAACATCTACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCCCCATCGTAATGGATTGCGCATCGTTGATAACAtt                                <  1:4468185/93‑1 (MQ=255)

                                                         cGCTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACAtt                                >  2:5108723/1‑93 (MQ=255)

                                                           cTGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTca                              <  2:628727/93‑1 (MQ=255)

                                                            tGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGCTCGCGCATCGTTGATAACATTcac                             >  1:4058022/1‑93 (MQ=255)

                                                            tGGTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTcac                             <  2:5896296/93‑1 (MQ=255)

                                                             ggTAACAACGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTcaca                            >  1:833667/1‑93 (MQ=255)

                                                                aacaacGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACaaaa                         >  1:2388911/1‑93 (MQ=255)

                                                                 acaacGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACaaaaa                        >  1:1078765/1‑93 (MQ=255)

                                                                 acaacGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACaaaaa                        >  2:2795123/1‑93 (MQ=255)

                                                                  caacGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAg                       <  2:3914623/93‑1 (MQ=255)

                                                                  caacGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAg                       <  1:4320714/93‑1 (MQ=255)

                                                                    acGACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTg                     <  1:1637755/93‑1 (MQ=255)

                                                                     cgcccACCGTCCCCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGa                    <  1:3859034/90‑1 (MQ=255)

                                                                      gACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGAc                   <  1:2512846/93‑1 (MQ=255)

                                                                      gACCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGAc                   <  1:1724982/93‑1 (MQ=255)

                                                                       aCCATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAAAATTCACAAAAAGTGAcc                  >  2:5362754/1‑93 (MQ=255)

                                                                        ccATCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCt                 <  1:3675240/93‑1 (MQ=255)

                                                                         cATCGTCACCTCGCTGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTc                <  1:2760994/93‑1 (MQ=255)

                                                                         cATCGTCACCTCGCGGTTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTc                <  2:3691634/93‑1 (MQ=255)

                                                                          aTCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAAGATTCACAAAAAGTGACCTCa               <  1:6325364/93‑1 (MQ=255)

                                                                          aTCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCa               >  1:5565862/1‑93 (MQ=255)

                                                                          aTCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCa               >  1:5759898/1‑93 (MQ=255)

                                                                          aTCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCa               >  2:1932266/1‑93 (MQ=255)

                                                                          aTCGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCa               >  1:2451355/1‑93 (MQ=255)

                                                                            cGTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGc             <  1:3526099/93‑1 (MQ=255)

                                                                             gTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCa            <  2:2260137/93‑1 (MQ=255)

                                                                             gTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCa            >  1:1155618/1‑93 (MQ=255)

                                                                             gTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCa            >  2:5488301/1‑93 (MQ=255)

                                                                             gTCACCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCa            >  2:2448181/1‑93 (MQ=255)

                                                                                aCCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACcg         >  1:4434431/1‑93 (MQ=255)

                                                                                aCCTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACcg         >  2:3375480/1‑93 (MQ=255)

                                                                                 ccTCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACcgc        >  2:3422108/1‑93 (MQ=255)

                                                                                   tCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACcgcgc      >  1:5126334/1‑93 (MQ=255)

                                                                                   tCGCGGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACcgcgc      <  2:3615056/93‑1 (MQ=255)

                                                                                      cgGGTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACCGCGCGcc   >  1:6089749/1‑93 (MQ=255)

                                                                                       gggTTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACCGCGCGCCg  >  2:2641118/1‑93 (MQ=255)

ACCAGTGCTCACGATACCATTGTGGGATCAGCGACCAGAGTTGCTGCAACATTTCACCGCTGGTAACAACGACCATCGTCACCTCGCGGGCTTGCCAGCGCATCGCTTCTTCTGCACCATCGTAATGGATTGCGCATCGTTGATAACATTCACAAAAAGTGACCTCAGCACCGCGCGCCG  >  NC_000913/3986826‑3987005

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