BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,779,837	T→C	H118R (CAC→CGC)	cysE ←	serine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,779,837	0	T	C	100.0%	853.7 / NA	266	H118R (CAC→CGC)	cysE	serine acetyltransferase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (131/135); total (131/135)

Only 100 of 271 total aligned reads displayed.

GCTGCCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGTGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTACGC  >  NC_000913/3779747‑3779926

 ctgcCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGAc                                                                                        <  2:1156111/93‑1 (MQ=255)

   gcCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCg                                                                                      >  2:2805620/1‑93 (MQ=255)

   gcCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCg                                                                                      >  2:2330807/1‑93 (MQ=255)

   gcCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCg                                                                                      >  2:5455875/1‑93 (MQ=255)

     cGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGAt                                                                                    <  1:1682522/93‑1 (MQ=255)

       ggTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGc                                                                                  >  1:4604262/1‑93 (MQ=255)

       ggTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGc                                                                                  <  2:4676086/93‑1 (MQ=255)

        gTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCg                                                                                 >  1:5105378/1‑93 (MQ=255)

         tGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGa                                                                                <  1:3781660/93‑1 (MQ=255)

           aaTATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATa                                                                              <  2:5374406/93‑1 (MQ=255)

              atCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGc                                                                           <  2:2653077/93‑1 (MQ=255)

               tCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGcc                                                                          >  1:5760135/1‑93 (MQ=255)

                cGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCt                                                                         >  1:677766/1‑93 (MQ=255)

                cGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCt                                                                         <  2:3508875/93‑1 (MQ=255)

                  aCCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGc                                                                       >  1:5443182/1‑93 (MQ=255)

                  aCCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGc                                                                       <  2:4422111/93‑1 (MQ=255)

                  aCCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGc                                                                       >  1:1378392/1‑93 (MQ=255)

                  aCCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGc                                                                       <  2:1296285/93‑1 (MQ=255)

                    cTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCaa                                                                     >  2:5971445/1‑93 (MQ=255)

                    cTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCaa                                                                     >  2:4751104/1‑93 (MQ=255)

                    cTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCaa                                                                     <  1:4161757/93‑1 (MQ=255)

                    cTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCaa                                                                     >  2:2948538/1‑93 (MQ=255)

                     tGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGTCTGCAAg                                                                    >  2:2398534/1‑93 (MQ=255)

                     tGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAg                                                                    >  1:2215305/1‑93 (MQ=255)

                      ggAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAgg                                                                   >  2:861393/1‑93 (MQ=255)

                       gAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGc                                                                  <  2:3591560/93‑1 (MQ=255)

                       gAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGc                                                                  >  1:793222/1‑93 (MQ=255)

                       gAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGc                                                                  >  2:1530589/1‑93 (MQ=255)

                        aaCGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCa                                                                 <  2:3438349/93‑1 (MQ=255)

                        aaCGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCa                                                                 >  2:6193661/1‑93 (MQ=255)

                         aCGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCGCTGATTCCACAACCACCGACCGATGCGATAGGCCTGCAAGGCAt                                                                >  1:538000/1‑93 (MQ=255)

                         aCGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCAt                                                                <  2:1351667/93‑1 (MQ=255)

                         aCGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCAt                                                                <  2:4048096/93‑1 (MQ=255)

                           gTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGa                                                              >  2:2629374/1‑93 (MQ=255)

                           gTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGa                                                              >  2:4285407/1‑93 (MQ=255)

                              acaGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGaaaa                                                           >  2:1801078/1‑93 (MQ=255)

                               caGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAAc                                                          >  2:702071/1‑93 (MQ=255)

                                 gAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAAccc                                                        >  2:1129655/1‑93 (MQ=255)

                                 gAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAAccc                                                        <  1:4789552/93‑1 (MQ=255)

                                  aaaCCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCt                                                       >  1:4221946/1‑93 (MQ=255)

                                  aaaCCTGGTTTTGCAGAAAGATTGCCAGGGCCCGACCCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCt                                                       <  1:1611239/93‑1 (MQ=255)

                                    aCCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTc                                                     <  2:1980777/93‑1 (MQ=255)

                                    aCCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTc                                                     <  2:4902260/93‑1 (MQ=255)

                                     ccTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCa                                                    <  1:3713902/93‑1 (MQ=255)

                                      cTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAg                                                   <  1:6151586/93‑1 (MQ=255)

                                      cTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAg                                                   >  1:3358787/1‑93 (MQ=255)

                                       tGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAgg                                                  <  1:1518714/93‑1 (MQ=255)

                                       tGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAgg                                                  >  1:5871897/1‑93 (MQ=255)

                                       tGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAgg                                                  >  1:6213335/1‑93 (MQ=255)

                                       tGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAgg                                                  >  1:1288183/1‑93 (MQ=255)

                                        ggTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGt                                                 <  2:4313715/93‑1 (MQ=255)

                                        ggTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGt                                                 >  2:1095393/1‑93 (MQ=255)

                                         gTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGta                                                >  1:2040329/1‑93 (MQ=255)

                                           tttGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGtata                                              >  2:6282874/1‑93 (MQ=255)

                                            ttGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATaa                                             <  2:2909897/93‑1 (MQ=255)

                                             tGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATaac                                            >  1:4723811/1‑93 (MQ=255)

                                               cAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAGCCCTTCAGGTATaacaa                                          <  2:5068299/93‑1 (MQ=255)

                                               cAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATaacaa                                          <  2:3822405/93‑1 (MQ=255)

                                                aGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGACACCCCTTCAGGTATaacaac                                         >  2:4521585/1‑93 (MQ=255)

                                                 gAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACg                                        >  2:3068729/1‑93 (MQ=255)

                                                 gAAAGATTGCCAGTCCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACg                                        <  1:2358960/93‑1 (MQ=255)

                                                   aaGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACggg                                      >  2:5284544/1‑93 (MQ=255)

                                                   aaGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACggg                                      <  1:1033657/93‑1 (MQ=255)

                                                    aGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACgggg                                     >  1:6170241/1‑93 (MQ=255)

                                                     gATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGt                                    <  2:5783550/93‑1 (MQ=255)

                                                     gATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGt                                    <  2:4833766/93‑1 (MQ=255)

                                                     gATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGt                                    >  2:1629858/1‑93 (MQ=255)

                                                       ttGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTg                                  >  1:71030/1‑93 (MQ=255)

                                                       ttGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTg                                  >  1:4079093/1‑93 (MQ=255)

                                                       ttGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTg                                  >  1:808380/1‑93 (MQ=255)

                                                        tGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGa                                 >  2:1500043/1‑93 (MQ=255)

                                                        tGCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGa                                 <  2:5906386/93‑1 (MQ=255)

                                                         gCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAg                                >  1:6277368/1‑93 (MQ=255)

                                                         gCCAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAg                                <  1:2545956/93‑1 (MQ=255)

                                                           cAGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTa                              >  1:2941043/1‑93 (MQ=255)

                                                            aGTGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTAt                             <  1:327223/93‑1 (MQ=255)

                                                             atGCGCGACGCCCCTAATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACACTTCAGGTATAACAACGGGGTTGAGTAtt                            <  2:2070017/92‑1 (MQ=255)

                                                              tGCGCGACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTAttt                           >  2:78567/1‑93 (MQ=255)

                                                               gcggGACGCCCCCGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAACCCCTTCAGGTATAACAACGGGGTTGAGTATTTa                          <  2:4366864/93‑1 (MQ=255)

                                                               gcgcGACTCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTa                          <  1:3682863/93‑1 (MQ=255)

                                                                cgcgACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTAt                         <  1:2846464/93‑1 (MQ=255)

                                                                cgcgACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTAt                         >  1:2932970/1‑93 (MQ=255)

                                                                   gACGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGa                      <  1:3528389/93‑1 (MQ=255)

                                                                    aCGCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGAc                     <  2:2208997/93‑1 (MQ=255)

                                                                      gCCCCTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATGTATCGACTg                   >  1:384550/1‑93 (MQ=255)

                                                                        cccTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGcc                 <  2:5224758/93‑1 (MQ=255)

                                                                          cTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCgg               >  2:3261676/1‑93 (MQ=255)

                                                                          cTGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCgg               <  2:1825946/93‑1 (MQ=255)

                                                                           tGATTCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCggg              <  1:2467120/93‑1 (MQ=255)

                                                                              ttCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGATGAGTATTTATCGACTGCCGGGTcg           >  1:1486491/1‑93 (MQ=255)

                                                                               tCCACAACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTcgc          <  2:581797/93‑1 (MQ=255)

                                                                                 cacaACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCgg        <  1:2586980/93‑1 (MQ=255)

                                                                                 cacaACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCgg        >  1:4390122/1‑93 (MQ=255)

                                                                                 cacaACCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCgg        <  1:889974/93‑1 (MQ=255)

                                                                                    aaCCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTa     <  1:3973292/93‑1 (MQ=255)

                                                                                    aaCCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTa     <  1:4341830/93‑1 (MQ=255)

                                                                                     aCCAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTAc    <  1:2752601/93‑1 (MQ=255)

                                                                                      ccAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTACg   <  1:3072002/93‑1 (MQ=255)

                                                                                      ccAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTACg   >  1:1864412/1‑93 (MQ=255)

                                                                                       cAGCGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTACGc  >  1:1133331/1‑93 (MQ=255)

GCTGCCGGGTGAATATCGACCTGGAACGTCACAGAAACCTGGTTTTGCAGAAAGATTGCCAGTGCGCGACGCCCCTGATTCCACAACCAGTGACCGATGCGATAGGCCTGCAAGGCATGAAAACCCTTCAGGTATAACAACGGGGTTGAGTATTTATCGACTGCCGGGTCGCGGGTACGC  >  NC_000913/3779747‑3779926

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